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THERMO-MECHANICAL  AND  THERMAL 


BEHAVIOR  OF  HIGH-TEMPERATURE  STRUCTURAL  MATERIALS 


PREFACE 


Technical  ceramics  because  of  their  chemical  inertness,  high  melting 
point,  good  wear  resistance,  excellent  mechanical  stability  at  high  tempera¬ 
ture  and  other  unique  properties,  represent  a  class  of  materials  eminently 
suited  for  many  critical  engineering  applications.  Unfortunately,  because 
of  their  brittleness  and  unfavorable  combination  of  pertinent  material 
properties,  technical  ceramics  generally  are  highly  susceptible  to  catas¬ 
trophic  failure  in  non-uniform  thermal  environments,  which  give  rise  to 
thermal  stresses  of  high  magnitude. 

Thermal  stress  failure  analysis  of  structural  materials  represents 
a  multi-disciplinary  problem  which  involves  the  principles  of  heat  trans¬ 
fer,  mechanics  and  materials  engineering.  Over  the  last  few  decades  much 
general  understanding  of  the  nature  of  thermal  stress  failure  of  brittle 
materials  has  been  generated.  However,  due  to  the  multi-disciplinary 
nature  of  the  problem,  the  ability  to  predict  thermal  stress  failure  quan¬ 
titatively  for  design  or  other  purposes  has  lagged  behind  the  progress  made 
in  other  engineering  fields.  The  objective  of  the  present  program  is  to 
improve  the  qualitative  and  quantitative  understanding  of  the  nature  of 
thermal  stress  failure  of  brittle  structural  materials,  including  the  ex¬ 
perimental  as  well  as  theoretical  variables.  In  order  to  achieve  this 
objective,  the  participating  investigators  and  scope  of  the  program  are 
organized  such  that  full  advantage  is  taken  of  the  combined  inputs  from  a 
number  of  engineering  disciplines.  In  a  similar  spirit,  a  number  of  studies 
were  conducted  in  cooperation  with  investigators  at  other  institutions. 

The  effort  of  this  program  consists  of  four  main  themes,  including: 
experimental  thermal  shock  testing  with  supporting  analyses,  measurement 
of  thermophysical  properties  relevant  to  thermal  stress  failure,  the 
analysis  of  mechanisms  of  thermal  stress  failure  and  the  dissemination 
of  information  on  thermal  stresses  in  the  form  of  review  articles,  con¬ 
ferences,  etc. 

Studies  completed  within  the  period  covered  by  this  report  are  pre¬ 
sented  as  individual  chapters  in  the  main  body  of  this  report.  The  title 
of  these  chapters  within  their  main  theme  together  with  a  brief  comment 
are  as  follows: 


A.  Experimental  Thermal  Shock  Testing  and  Suoporting  Analyses 


Chapter  I:  "Observation  on  the  Characteristics  of  a  Fluidized 

Bed  for  the  Thermal  Shock  Testing  of  Brittle  Ceramics" 


This  study  indicates  that  the  heat  transfer  coefficient  of  a  fluidized 
bed  depends  strongly  on  air  flow  rate  and  powder  particle  size.  The  absolute 
values  for  the  heat  transfer  coefficient  for  the  fluidized  bed  are  below 
those  for  silicone  oils  and  aqueous  media.  For  this  reason  fluidized  beds 
extend  the  range  of  values  of  heat  transfer  coefficient  available  to  the 
investigator. 


Chapter  II:  "Effect  of  AT-  and  Spatially  Varying  Heat  Transfer 
Coefficient  on  Thermal  Stress  Resistance  of  Brittle 
Ceramics  Measured  by  the  Quenching  Method" 


Commonly,  analyses  of  experimental  data  for  the  thermal  stress  re¬ 
sistance  of  brittle  ceramics  measured  by  the  quenching  method  assume 
that  the  heat  transfer  coefficient  is  Newtonian  in  nature  and  is  spatially 
uniform  over  the  whole  surface  of  the  specimen.  In  practice,  these  assump¬ 
tions  generally  are  not  realized.  The  results  of  this  study  shows  that 
for  an  accurate  assessment  of  thermal  shock  data,  these  effects  must  be 
taken  into  account. 


Chapter  III:  "Analysis  of  Effect  of  Crack  Interaction  on  Nature  of 
Strength  Loss  in  Thermally  Shocked  Brittle  Ceramics" 


The  results  of  this  study  indicate  that  co-planar  crack  interaction 
during  thermal  stressing  of  brittle  ceramics  can  lead  to  crack  coalescence 
and  complete  disintegration,  not  predicted  by  fracture  theories  which  do 
not  take  crack  interaction  into  account. 


Chapter  IV:  "Analysis  of  Thermal  Fatigue  Behavior  of  Brittle 
Structural  Materials" 


The  purpose  of  this  study  was  to  define  the  proper  role  of  the  various 
fracture  mechanical  parameters  such  as  stress  intensity  factor  exponent  and 
activation  energy  on  the  thermal  fatigue  resistance  of  brittle  structural 
materials  for  a  range  of  defined  thermal  environments. 


B.  Analysis  of  Mechanisms  of  Thermal  Stress  Failure 

Chapter  V:  "Thermal  Stresses  in  a  Partially  Absorbing  Plate 
Asymmetrically  Heated  by  Cyclic  Thermal  Radiation 
and  Cooled  by  Convection” 


This  study  is  a  continuation  of  a  series  of  analyses  of  thermal 
stresses  in  partially  absorbing  ceramics  subjected  to  radiation  heating. 

The  present  study  focussed  on  the  effect  of  cyclic  radiation.  The  results 
obtained  show  that  for  the  same  average  radiative  flux,  cyclic  radiation 
will  result  in  thermal  stresses  of  higher  magnitude  than  steady-state 
radiation. 

C.  Measurement  of  Thermophysical  Properties  Relevant  to  Thermal 
Stress  Failure 

Chapter  VI:  "Anisotropy  Effects  in  the  Thermal  Diffusivity  of  Si^N^-BN 
Composites" 


Experimental  data  for  the  thermal  diffusivity  of  hot-pressed  S^N^-BN 
composites  show  that  the  heat  conduction  properties  of  these  composites 
can  be  highly  anisotropic  due  to  the  preferred  orientation  of  the  BN  inclusions. 

Chapter  VII:  "Thermal  Diffusivity  of  Silicon  Carbide-Silicon  Composites" 


The  results  of  this  study  indicate  that  the  room  temperature  thermal 
diffusivity  of  silicon  carbide-silicon  composites  is  strongly  dependent  on 
the  impurity  content  of  the  silicon  carbide  and  silicon  used  for  their 
manufacture. 


Chapter  VIII:  "Effect  of  Microcracking  on  the  Conduction  of  Heat  in 
Brittle  Composites" 


The  thermal  diffusivity  of  composites  of  silicon  carbide  dispersions 
in  matrices  of  aluminum  oxide,  beryllium  oxide  or  magnesium  oxide  is 
lowered  significantly  by  the  formation  of  microcracks  which  result  from 
the  mismatch  in  the  coefficients  of  thermal  expansion.  The  experimental 
results  also  indicate  that  microcrack  formation  is  a  highly  sensitive 
function  of  pore  content. 


Chapter  IX:  "Measurement  of  the  Thermal  Conductivity  and  Diffusivity 
of  Fine  Silicon  Carbide  Fibers  by  the  Composite 
Technique” 


Values  were  obtained  for  the  thermal  conductivity  and  diffusivity  of 
amorphous  silicon  carbide  fibers  from  corresponding  experimental  data  for 
composites  composed  of  these  fibers  contained  in  a  lithium-alumino-silicate 
glass  ceramic  matrix. 

D.  Reviews 


Chapter  X:  "Thermal  Stress  in  Materials  Heated  Internally  by 
Radiation  Absorption" 


The  results  of  Chapter  V  and  earlier  analysis  of  thermal  stresses  in 
semi-absorbing  ceramics  were  primarily  aimed  at  electro-magnetic  radiation 
such  as  solar  or  laser  radiation.  However,  the  results  obtained  are  identical 
to  the  thermal  stresses  which  would  result  from  the  absorption  of  nuclear 
radiation.  This  chapter  (X)  represents  an  overview  of  all  relevant  results, 
presented  to  the  technical  community  interested  in  fusion  research. 

For  the  convenience  of  the  reader,  Table  1  lists  the  status  of  publi¬ 
cations  presented  in  earlier  reports.  For  most  of  these  studies,  reprints 
are  available  and  will  be  sent  on  request. 
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The  heat  transfer  characteristics  of  a  fluidized  bed  used  as  a 
quenching  medium  for  the  thermal  stress  testing  of  brittle  ceramics, 
were  determined  by  measurements  of  the  thermal  shock  behavior  of  rods  of 
a  soda-lime-silica  glass.  The  heat  transfer  coefficient  was  found  to 
be  strongly  dependent  on  the  mean  particle  size  of  the  powder  and  air 
flow  rate,  and  was  relatively  independent  of  the  position  within  the 
bed.  The  results  indicate  that  the  heat  transfer  coefficient  during 
thermal  shock  fracture  may  have  a  value  lower  than  that  obtained  under 
heat  transfer  conditions  which  more  closely  resemble  steady-state. 

The  heat  transfer  data  inferred  from  the  quenching  experiments  with 
the  glass,  gave  excellent  agreement  between  calculated  and  measured 
values  for  the  thermal  shock  behavior  for  rods  of  a  polycrystalline 
aluminum  oxide.  It  is  concluded  that  fluidized  beds  are  excellent  inert 
quenching  media  with  variable  heat  transfer  coefficient  controlled  by 
particle  size  and  flow  rate. 

1.  INTRODUCTION 

Technical  ceramics,  because  of  their  combination  of  high  melting 
point  and  chemical  stability  are  excellent  candidate  materials  for 
engineering  applications  involving  extreme  thermal  and  chemical 
environments.  However,  due  to  an  unfavorable  combination  of  values 
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for  such  properties  as  thermal  conductivity.  Young's  modulus,  coefficient 
of  thermal  expansion  and  tensile  strength,  many  technical  ceramics 
exhibit  low  resistance  to  failure  by  thermal  stresses,  frequently  en¬ 
countered  during  transient  or  steady-state  heat  flow  [1],  Due  to  the 
inherent  brittleness  of  ceramic  materials,  such  failure  can  be  highly 
catastrophic.  For  reliable  engineering  design  and  optimum  material 
selection,  it  is  imperative  that  the  thermal  stress  behavior  of  technical 
ceramics  be  well  understood.  This  can  be  accomplished  on  the  basis  of 
theoretical  principles  or  by  actual  testing.  One  popular  test  method, 
in  view  of  its  simplicity,  consists  of  quenching  specimens  in  the  form 
of  circular  or  square  cylinders  from  a  pre-selected  higher  temperature 
into  a  fluid  medium  at  lower  temperature.  Such  a  fluid  medium  can  con¬ 
sist  of  either  silicone  oil  [2],  a  eutectic  salt  mixture  [3,4],  water  [5], 
liquid  metals  [6]  or  fluidized  powders  [7].  For  purposes  of  the  quanti¬ 
tative  evaluation  of  the  test  results,  the  heat  transfer  characteristics 
of  these  fluid  media  must  be  well  understood.  In  this  respect,  a  number 
of  experimental  studies  [8,9,10,11]  have  shown  that  water  as  a  quenching 
medium  can  be  highly  variable  in  its  heat  transfer  characteristics, 
primarily  due  to  nucleate  boiling  and/or  steam  film  formation  which 
are  dependent  not  only  on  the  temperature  but  also  on  the  surface  condi¬ 
tion  of  the  material.  The  decrease  of  the  viscosity  of  silicone  oils 
with  increasing  temperature  causes  a  corresponding  increase  in  the  heat 
transfer  coefficient  [11].  Increased  turbulence  of  eutectic  salts  with 
increasing  temperature  tends  to  create  a  corresponding  increase  in  the 
rate  of  heat  transfer  [4],  A  recent  analysis  showed  that  the  results 
of  quenching  experiments  also  depend  on  whether  the  heat  transfer  occurs 
by  natural  or  forced  convection  [12].  Water  has  the  disadvantage 


of  promoting  stress  corrosion  [13].  Oils  are  subject  to  decomposition 
which  can  lead  to  the  formation  of  strongly  adhering  coatings  on  the 
specimen  surface  and  subsequent  changes  in  mechanical  or  other  pro¬ 
perties  [14]. 

Fluidized  beds  which  use  ceramic  powders  as  the  fluid  medium 
offer  the  advantage  of  chemical  inertness  over  a  wide  range  of  tem¬ 
perature.  Literature  data  [7,15]  show  considerable  differences  in  the 
heat  transfer  characteristics  of  fluidized  beds  used  for  thermal  shock 
testing.  These  differences,  as  indicated  by  the  experimental  data  of 
Callahan  [15]  for  the  heat  transfer  coefficient  of  fluidized  beds  obtained 
by  a  direct  method,  may  be  attributable  to  differences  in  the  way  these 
fluidized  bed  were  being  operated.  In  this  respect,  then,  the  primary 
objective  of  this  study  was  to  investigate  the  effect  of  gas  flow  rate 
and  powder  particle  size  on  the  heat  transfer  characteristics  of  fluidized 
beds  used  for  the  thermal  shock  testing  of  brittle  ceramics.  The  secondary 
objective  of  this  study  was  to  compare  the  values  of  the  heat  transfer 
coefficient  inferred  indirectly  from  the  thermal  shock  data  with  those 
measured  by  direct  means. 

2.  EXPERIMENTAL 
2.1  Test  Materials 

A  soda-lime-silica  glass  and  a  polycrystalline  alumina  were 
selected  as  materials  appropriate  for  testing  of  the  fluidized  bed.  Both 
these  materials  were  used  in  a  number  of  previous  studies  [11,13]  which 
facilitated  comparison  of  the  experimental  data.  The  materials  were 

R-6,  Owens-Corning  Fiberglass  Corp.,  Toledo,  Ohio 
**998  Alumina,  McDanel  Refractory  Porcelain  Co.,  Beaver  Falls,  PA. 


supplied  by  the  manufacturers  in  the  form  of  circular  rods.  Their 
property  values  and  diameter  are  given  in  Table  1.  Because  of  the 
ease  of  crack  detection  in  the  glass  by  visual  means,  most  tests  were 
conducted  with  this  material.  The  aluminum  oxide  specimens  were  used 
for  the  purpose  of  illustration  of  the  use  of  data  obtained. 

2.2  Equipment  and  Test  Procedures 

The  fluidized  bed  of  commercial  design5,  suitable  for  operation  at 
elevated  temperatures,  if  desired,  consisted  of  a  vertical  stainless- 
steel  hollow  cylinder  with  inner  diameter  of  approximately  0.1  m  and 
height  of  approximately  0.15  m.  In  order  to  extend  the  range  of  the 
rates  of  air  flow  and  to  minimize  the  loss  of  the  fluidizing  powder,  a 
silica  tube  approximately  0.5  m  high  and  0.09  m  in  inner  diameter  was 
placed  inside  the  steel  cylinder.  Air  was  used  as  the  fluidizing  agent. 

The  rate  of  air  flow  was  controlled  by  a  flow  meter.  To  assure  a  uniform 
flow  of  air  across  the  width  of  the  fluidizing  bed,  the  perforated  bottom 
was  covered  with  a  fiberglass  cloth  pad,  which  also  prevented  the  loss 
of  finer  powders  through  the  bed  bottom  between  experiments.  For  all  tests, 
the  fluidized  bed  was  maintained  at  room  temperature. 

•k 

An  aluminum  oxide  powder  was  selected  as  an  appropriate  fluidizing 
medium.  By  screening,  the  powder  was  divided  into  a  number  of  size  frac¬ 
tions  corresponding  to  screen  sizes  of  60  to  65,  100  to  115,  200  to  270 
and  325-400  mesh,  to  give  average  particle  sizes  of  231,  138,  69  and  38  um. 

§TECAM  Fluidized  Bath  Model  FBS,  Techne  Incorp.,  661  Brunswick  Pike, 
Princeton,  N.J. 

★ 

Calcined  alumina  powder  (-230  mesh),  Aluminum  Corporation  of  America. 
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TABLE  1.  Properties  of  Test  Materials 


★ 

Soda-Lime  Glass 


A1 umi na 


*★ 


Young's  modulus  (GPa)  69 

Coeff.  of  thermal  expansion,  a 

(°C-1)  x  10-6  9.3 

Poisson's  ratio,  v  0.25 

Thermal  conductivity,  k  , 

(cal.cm.-^C-'.s"1)  x  10J  2.5 

Flexural  strength  (MPa)  194 

(Liq.  N2  temp.) 

Diameter,  d(m)  x  10-3  5.3 

Weibull  parameter,  m  8 


393 

9  ( 725°C) 
0.265 

21.9  ( 725°C) 
271 

( Room  temp . ) 
6.35 
16 


See  refs.  11 ,  13. 

**  Young's  modulus  value  was  supplied  by  the  manufacturer,  Poisson's 
ratio  and  coeff.  of  thermal  expansion  were  obtained  from  ref.  16, 
and  other  properties  values  were  obtained  in  the  writer's  laboratory. 
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Prior  to  quenching,  the  specimens  were  preheated  in  an  electrically 
heated,  thermostatically  controlled  (±  0.5°C)  tube  furnace  with  internal 


diameter  of  -  0.05  m,  mounted  vertically  above  the  fluidized  bed.  The 
specimens  in  the  form  of  circular  cylinder  with  length  -  0.06  m  were 
attached  to  the  lower  end  of  a  quartz  rod  mounted  vertically  along  the 
center  line  of  the  furnace  tube.  In  the  upper  position  of  the  rod, 
during  preheating,  the  specimens  were  located  in  the  axial  and  radial 
center  of  the  hot  zone  of  the  furnace.  In  the  lower  position  of  the 
rod,  the  specimen  was  located  within  the  fluidized  bed.  This  lower 
position  was  adjustable  to  permit  measurement  of  the  heat  transfer 
characteristics  as  a  function  of  the  vertical  position  within  the 
bed.  Caps  mounted  over  the  ends  of  the  furnace  tube  assured  temperature 
uniformity  and  eliminated  possible  specimen  cooling  due  to  the  upward 
flow  of  air  through  the  tube  furnace.  Thermal  equilibrium  immediately 
prior  to  quenching  was  obtained  by  holding  the  specimens  for  at  least  15 
min.  at  the  level  of  preselected  temperature.  Quenching  of  the  specimens 
was  accomplished  by  dropping  the  rod  from  the  upper  to  lower  position 
under  the  influence  of  gravity  without  any  mechanical  constraint. 
Following  the  quenching,  the  specimens  were  held  in  the  fluidized  bed 
for  approximately  5  min.  after  which  they  were  removed  and  examined.  All 
specimens  were  tested  singly,  in  order  to  assure  that  the  heat  transfer 
was  not  affected  by  the  presence  of  other  specimens.  However,  in  order 
to  take  account  of  possible  specimen-to-specimen  variation,  a  total  of 
five  (5)  specimens  were  tested  for  any  given  temperature  difference  and 
operating  condition  of  the  fluidized  bed.  To  eliminate  any  fatigue 
effects  on  the  results  obtained  each  specimen  was  tested  only  once. 
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The  critical  temperature  difference  (aTc)  required  to  induce  thermal 
stress  fracture  in  the  specimens  for  a  given  air  flow  rate  and  alumina 
particle  size  was  established  by  varying  the  temperature  difference 
between  the  tube  furnace  and  the  fluidized  bed.  To  reflect  the  variation 
from  specimen  to  specimen,  ATc  for  the  glass  was  reported  in  terms  of 
the  temperature  range  required  to  fracture  one  and  four  of  the  five 
specimens  tested.  For  the  alumina  specimens  in  which  the  cracks  were 
not  easily  detected  by  visual  or  other  means,  aTc  was  determined  by 
measuring  the  fracture  strength  of  the  quenched  specimens  at  room 
temperature  in  four-point  bending  with  an  inner  and  outer  span  of  0.016 
and  0.05  m,  respectively.  For  these  specimens,  ATC  was  defined  as  the 
value  of  quenching  temperature  difference  which  caused  a  50%  or  greater 
decrease  in  strength  compared  to  the  non-quenched  specimens. 

The  heat  transfer  characteristics  of  the  fluidized  bed  were  deter¬ 
mined  by  indirect  as  well  as  direct  means.  Indirectly,  the  heat 
transfer  coefficient  was  determined  from  the  experimental  data  for 
ATc  for  the  glass  specimens  with  the  aid  of  thermoelastic  theory.  Thermal 
stress  failure  in  brittle  materials  generally  occurs  under  the  influence 
of  the  tensile  stresses.  For  the  quenching  conditions  of  the  present 
study,  the  maximum  tensile  thermal  stresses  occur  in  the  specimen  surface. 
For  a  solid  circular  cylinder  and  0  <  Rh/k  <  10,  from  the  solutions  of 
Jaeger  [17]  the  magnitude  of  these  stresses  to  a  very  good  approximation 
can  be  described  by: 


i(r  .  R,  .  L4^  {,♦»!£  }  0, 

where  a,  E,  K  and  v  are  defined  in  Table  1,  R  is  the  radius  of  the  cylinder 
and  h  is  the  heat  transfer  coefficient.  At  aT  =  aTc,  the  maximum  value 


of  tensile  thermal  stress  equals  the  tensile  fracture  stress.  Following 
the  approach  of  Becher  et  al.  [10]  and  Singh,  Tree  and  Hasselman  [11],  the 
heat  transfer  coefficient  was  calcuated  by  substitution  of  the  values 
for  the  material  properties,  R  and  aTc  into  equation  1.  In  the  calculations, 
the  relative  temperature  dependence  of  the  strength  of  the  glass  specimen 
was  assumed  to  be  identical  to  the  strength  data  of  Shand  [18].  Using 
the  strength  of  194  MPa  at  liquid  nitrogen  temperature  [13],  an  estimate 
of  the  strength  of  the  present  glass  specimens  at  different  temperatures 
was  easily  made. 

The  direct  measurement  of  th,*  heat  transfer  coefficient  relied  on 
the  method  of  Callahan  [15]  and  Trantina  [14],  which  consisted  of  measuring 
the  transient  temperature  respor  n  of  a  small,  preheated  specimen  quenched 
in  the  fluidized  bed.  For  tne  present  study  the  specimen  consisted 
of  an  alumina  sphere  of  diameter  1.25  cm,  with  a  thermocouple  in  a  hole 
drilled  to  the  exact  center  of  the  sphere.  For  the  alumina,  the  thermal 
conductivity  and  diffusivity  were  established  independently  by  the  laser- 
flash  method.  A  plot  of  the  non-dimensional  transient  temperature  at 
the  center  of  the  sphere,  on  comparison  with  the  charts  of  Schneider  [19], 
yields  the  appropriate  value  of  Rh/K  from  which  the  value  of  h  can 
be  calculated. 

The  value  of  h  obtained  from  the  data  for  aTc  may  possibly  reflect 
the  highly  transient  nature  (if  any)  of  the  heat  transfer  during  the 
quench.  The  value  of  h  obtained  by  the  immersed-sphere  method  requires 
a  time  period  in  excess  of  40  seconds  well  in  excess  of  the  time  period 
of  <  0.1  seconds  of  insertion  of  the  specimen  in  the  fluidized  bed  or 
the  time  period  of  ^  0.3  sec.  [20]  for  the  tensile  thermal  stress  to 


reach  its  maximum  value.  For  this  reason,  the  value  of  h  obtained  by  the 
latter  method  more  closely  corresponds  to  the  steady-state  heat  transfer 
conditions. 

3.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

Figure  1  shows  the  values  of  &TC  for  the  glass  specimens  as  a  func¬ 
tion  of  air  flow  rate  for  four  values  of  mean  particle  size  of  the  aluminum 
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oxide  powder.  At  approximately  0.8  l  min.  cm  the  air  flow  rate  was 
sufficiently  high  that  some  of  the  powder  with  the  lowest  value  of 
mean  particle  size  was  carried  out  of  the  bed.  The  data  indicate  that 
aT£  depends  strongly  on  the  air  flow  rate  as  well  as  particle  size. 

Such  variation  can  only  be  attributed  to  the  strong  dependence  of  the 
heat  transfer  coefficient  (h)  on  flow  rate  and  particle  size. 

The  validity  of  this  latter  conclusion  is  substantiated  by  the  values 
of  the  heat  transfer  coefficient  shown  in  Fig.  2  evaluated  from  the 
data  for  aTc  of  Fig.  1.  These  values  of  h  vary  by  nearly  as  much  as 
a  factor  of  three.  At  the  lower  range  of  flow  rate,  the  increasing 
value  of  h  with  increasing  flow  rate  most  likely  is  the  result  of  the 
corresponding  increasing  degree  of  fluidization  of  the  bed.  The  decrease 
in  h  with  increasing  flow  rate  at  the  higher  values  of  flow  rate  probably 
results  from  the  decrease  in  net  density  of  the  bed,  which  as  judged  by 
the  observed  increase  in  relative  bed  height  for  the  highest  flow  rate 
could  amount  to  as  much  as  30%.  The  relative  dependence  of  heat  transfer 
coefficient  on  flow  rate  as  shown  in  Fig.  2,  is  similar  to  the  findings 
of  Callahan  [15]. 


Figure  3  shows  the  dependence  of  the  heat  transfer  coefficient 


inferred  from  corresponding  data  for  aTc  for  the  glass  specimens  as  a 

function  of  the  depth  of  the  bed  for  the  specific  value  of  air  flow  rate 
-1  -2 

of  0.47  Ji.min  cm  and  average  alumina  particle  size  of  138  pm.  For 
convenience,  the  depth  of  the  bed  was  defined  as  the  distance  from  the 
bottom  of  the  bed  to  the  center  of  the  vertical  specimens  of  length  ^  6  cms. 
Experimentally,  it  was  observed  that  the  specimen  did  not  undergo  fracture 
at  a  preferred  location  along  their  length.  This  observation  combined 
with  the  relative  independence  of  the  heat  transfer  coefficient  on  depth 
of  the  bed  suggests  that  the  heat  transfer  appears  to  be  relatively 
uniform  along  the  length  of  the  specimen.  Of  practical  significance  is 
that  the  position  of  the  specimen  in  the  bed  is  not  critical  to  the 
results  obtained.  Of  course,  insertion  of  the  specimens  close  to  or  in 
direct  contact  with  the  side-walls  should  be  avoided. 

Figure  4  shows  the  experimental  data  for  the  heat  transfer  coefficients 
obtained  by  the  sphere  immersion  method  for  three  of  the  identical  batches 
of  powders  used  to  establish  the  data  shown  in  Figs.  1  and  2.  Comparison 
of  Figs.  2  and  4  shows  that  the  values  of  the  heat  transfer  coefficient 
obtained  by  both  methods  show  the  same  relative  dependence  on  the  air  flow 
rate.  Quantitatively,  however,  the  maximum  values  of  h  inferred  from 
the  data  for  aTc  significantly  exceed  the  corresponding  values  obtained 
by  the  sphere  immersion  method.  In  explaining  this  discrepancy  at  least  two 
major  effects  must  be  considered.  The  first  effect  results  from  the 
differences  in  the  geometry  of  the  test  specimens.  As  presented  by 
Gebhart  [21], the  forced  convection  heat  transfer  coefficients,  hc  and  h$ 
for  a  cylinder  and  sphere,  respectively  can  be  expressed  by: 


for  40  <  Re  <  4,000 


(2) 


hc  =  0.69(Re)°-466(Pr)1/3Kf/dc 

h$  =  [2.0  +  0.6{Re)0‘5(Pr)1/3]Kf/ds  for  1  <  Re  *  70,000  (3) 

where  Re  =  Vd/n  is  the  Reynolds  number,  V  is  the  fluid  velocity,  n 
is  the  fluid  viscosity  and  d  is  the  specimen  diameter.  Pr  is  the  Prandtl 
number,  is  the  thermal  conductivity  of  the  fluid  and  dc  and  d$  are  the 
diameters  of  the  cylinder  and  the  sphere  respectively. 

Taking  the  case  that  2  «  0.6  (Re)^'^(Pr)^3  and  assuming  no 
-significant  difference  between  the  exponents  of  0.466  and  0.5  for  the 
Reynolds  number  in  eqs.  2  and  3,  yields  the  ratio 

1/2 

hc/hs  =1.16(ds/dc)  (4) 

which  for  the  present  specimens  with  ds/dc  -  2.29  yields: 

hc/hs-1.74  (5) 

This  result  suggests  that  an  upward  adjustment  of  the  data  of  Fig.  4  for 

the  spherical  geometry  by  a  factor  1.74  yields  the  appropriate  values 

for  the  heat  transfer  coefficient  for  the  vertical  glass  circular  cylinders. 

-2  -1  -1 

For  example,  the  value  of  0.0145  cal. cm.  °C  s.  (Fig.  4)  inferred 

from  immersed  sphere  data  for  fluidized  bed  with  particle  size  of  69  urn 

-1  -2 

and  air  flow  rate  of  0.47  a.min.  cm.  must  be  multiplied  by  1.74 
in  order  to  account  for  geometry  effect.  This  brings  the  peak  value  to 
about  0.025  cal .cm.~2°c.~1s.~^  which  is  still  only  about  71%  of  the 
corresponding  value  for  h  determined  from  aTc  (Fig.  2  ). 

The  second  effect  to  which  the  remaining  discrepancy  between  the  data 
of  Fig.  2  and  the  adjusted  data  of  Fig.  4,  can  be  attributed  is  based  on 
the  statistical  nature  of  brittle  fracture,  which  states  that  the  strength 
of  a  brittle  solid  is  a  function  of  the  volume  of  the  stressed  material  as 
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well  as  the  stress  distribution.  Theoretical  analyses  of  this  effect 
have  been  presented  by  Wei bul 1 [22 ]  and  others.  As  a  result  of  this 
statistical  effect,  the  tensile  fracture  stress  of  the  glass  specimens 
under  conditions  of  thermal  shock  with  a  uniform  biaxial  stress  distri¬ 
bution  in  the  specimen  surface  is  expected  to  differ  from  the  strength 
under  conditions  of  4-point  bending  on  which  the  data  of  Fig.  2  were 
based. 

In  a  previous  study  [13]  of  the  thermal  fatigue  behavior  of  the 
identical  glass  specimens  of  the  present  study,  it  was  shown  that  the 
ratio  of  the  tensile  fracture  stress  (ots)  during  thermal  shock  to  the 
corresponding  fracture  stress  in  4-point  bending  (o^)  for  a  population  of 
surface  flaws,  by  means  of  the  Weibull  theory  can  be  calculated  to 
be: 

°ts/ob  “  °-b9 

This  result  suggests  that  multiplication  of  the  data  shown  in 
Fig.  2  by  the  factor  0.59,  will  yield  values  of  the  heat  transfer  coeffi¬ 
cient  which  more  closely  correspond  to  the  magnitude  of  the  tensile 
failure  stress  encountered  during  quenching  in  the  fluidized  bed.  The 

peak  value  for  the  heat  transfer  coefficient  h  determined  from  ATc  for 

-2  -2  -1  -1 

the  glass  specimens  then  becomes  approximately  2.1x10  cal. cm  °C  .s  . 

The  corresponding  h  value  obtained  by  the  correction  for  geometry  of  the 

-p 

results  in  Fig.  4  evaluated  from  the  sphere-immersion  data  is  2.5x10  cal. 
-2  -1  -1 

cm  .°C  .s  ,  which  is  higher  than  the  corrected  value  inferred  from 
aTc>  In  view  of  the  total  uncertainties  in  establishing  both  values, 
possibly  no  major  significance  should  be  attached  to  this  discrepancy. 
Nevertheless,  the  lower  value  of  h  inferred  from  aTc  may  possibly  reflect 


the  highly  transient  nature  of  the  heat  transfer  in  the  very  brief  time 
period  required  for  fracture  to  occur.  Possibly  in  this  brief  period, 
the  transfer  of  heat  from  the  specimen  can  not  benefit  from  the  motion 
of  the  fluidized  powder.  A  detailed  analysis  of  this  effect  is  beyond 
the  scope  of  the  present  discussion. 

Figure  5  shows  the  results  for  the  strength  behavior  of  the  alumina 

specimens  quenched  into  the  fluidized  bed  with  mean  particle  size  of 

-1  -2 

138  ym  and  air  flow  rate  of  0.47  t.min.  cm  .  The  dependence  of  strength 

on  temperature  difference  indicates  a  value  of  .'J  800.  From  the  data 

of  Figs.  2  and  4,  the  values  appropriately  adjusted  for  the  statistical 

and  geometry  effects  as  discussed  earlier  are  approximately  2.0  and  2.5.x 
-2  -2  -1 

10  cal. cm  °C . s  ,  respectively.  Using  these  latter  values,  the  data 
in  Table  1  and  ratio  of  strength  under  thermal  shock  to  the  strength  in 
bending,  o^/a^  ~  0.79  inferred  from  the  Wei  bull  parameter  for  a  dis¬ 

tribution  of  surface  flaws,  the  values  for  ,\T  can  be  calculated  to  be 
830°C  and  675°C.  The  value  of  725CC  for  the  property  values  listed  in 
Table  1,  corresponds  to  the  mean  specimen  temperature  at  the  instant  of 
fracture  for  aT  -  800°C.  The  much  better  agreement  between  this  latter 
value  and  830°C  predicted  from  the  h  value  inferred  from  aT  for  the 
glass  specimen  than  that  predicted  from  the  value  obtained  from  the 
sphere- immersion  method,  provides  confirmation  for  the  earlier  suggestion 
that  the  heat  transfer  in  a  fluidized  bed  during  a  thermal  quench  is  less 
severe  than  under  conditions  which  more  closely  resemble  steady-state 
conditions. 

A  comparison  of  data  shows  that  the  peak  value  of  the  heat  transfer 

-2  -2  -1 

coefficient  of  approximately  2.1  to  2.5  x  10  cal. cm.  °C.s  found  in 
the  present  study,  agrees  quite  well  with  the  value  of  approximately 
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2.0x10  cal. cm.  °C.s  reported  by  Ammann  et  al  [23]  and  Callahar  [15]. 

The  much  lower  value  of  about  8x10  cal. cm.  C. s’  found  by  Trantina  [14] 
probably  can  be  attributed  to  differences  in  experimental  conditions. 

Comparison  of  the  values  found  in  the  present  study  with  those  for 
other  quenching  media,  show  that  even  the  peak  values  for  the  heat  transfer 
coefficient  for  the  fluidized  bed  is  less  than  the  corresponding  values 
for  the  effective  heat  transfer  coefficients  for  silicone  oils  or  water 
obtained  in  a  previous  study  [11].  This  statement  holds  even  if  these 
latter  data  are  adjusted  downward  by  the  factor  0.59  to  account  for  the 
statistical  nature  of  brittle  fracture.  It  appears  then  that  fluidized 
beds  for  purposes  of  thermal  shock  testing  of  brittle  ceramics  extend  the 
range  of  heat  transfer  coefficients  available  to  the  researcher.  As  in¬ 
dicated  by  the  results  of  the  present  study  as  well  as  those  of  Callahan  [15] 
a  particular  advantage  of  fluidized  beds  is  that  they  offer  considerable 
flexibility  in  adjusting  the  heat  transfer  coefficient  by  modification  of 
gas  flow  rate  and  particle  size. 
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AIR  FLOW  RATE.  Gf  (Emin*1  cm'*) 


Effect  of  air  flow  rate  and  alumina  particle 
size  on  the  heat  transfer  coefficient  h  for 
fluidized  bed  as  inferred  from  the  results 
for  &T  in  figure  1 . 


DEPTH  (cm) 


Heat  transfer  coefficient  h  as  a  function  of  specimen 
position  in  the  fluidized  bed  for  air  flor  rate  of  0.47 
£.min.-lcm.~2  and  mean  alumina  particle  size  of  138  pm 
inferred  from  values  of  aTc  for  soda-lime  glass  rods. 


b 
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Fig.  4.  Heat  transfer  coefficient  h  for  the 
fluidized  bed  obtained  by  sphere- 
immersion  method  as  a  function  of 
air  flow  rate  and  alumina  particle 
size. 


QUENCHING  TEMP  DIFF,  AT  CC) 

Fig.  5.  Strength  as  a  function  of  quenching  temperature  dif¬ 
ference  for  alumina  rods  (d  =  0.63  cm)  quenched  in 
fluidized  bed  with  air  flow  rate  of  0.47  t.min.-1  .cm.'2 
and  mean  alumina  particle  size  of  138  pm. 
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ABSTRACT 


A  study  was  conducted  for  the  effect  of  a  spatial  variation  in  heat 
transfer  coefficient  during  forced  convection  and  the  effect  of  temperature 
dependence  of  the  heat  transfer  coefficient  during  free  convection  on  the 
magnitude  of  thermal  stress  encountered  during  the  thermal  shock  testing  of 
brittle  ceramics  by  the  quenching  method.  For  specimens  with  circular  geometry, 
the  results  obtained  by  the  finite  element  method  indicate  that  depending  on 
the  value  of  Biot  number,  a  factor  of  three  spatial  variation  in  heat 
transfer  coefficient  changes  the  tensile  thermal  stresses  on  cooling  by  a 
maximum  of  about  17%  over  the  value  for  spatially  uniform  heat  transfer. 

On  heating,  the  corresponding  tensile  thermal  stresses  are  lower  by  a 
maximum  of  about  6%. 

In  free  convection  for  a  heat  transfer  coefficient,  h,  proportional 
to  AT-1^4,  where  aT1-  is  the  instantaneous  temperature  difference  between 
the  specimen  and  quenching  medium,  the  stresses  are  appreciably  less  than 
those  for  aT^-  independent  heat  transfer  for  the  same  initial  value  of  h 
at  the  onset  of  the  thermal  quench.  The  validity  of  the  dependence  of  h 
on  aT^  was  established  experimentally.  The  relevance  of  the  above  result 
to  the  interpretation  of  thermal  quenching  studies  are  discussed. 
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I.  INTRODUCTION 


Brittle  ceramic  materials  generally  are  highly  susceptible  to  catastrophic 

fracture  under  conditions  of  thermal  shock  which  givesrise  to  thermal  stresses 
1  2 

of  high  magnitude.  ’  For  purpose  of  the  reliable  performance  of  ceramics 
for  structural  applications  their  response  to  thermal  stress  failure  must  be 
well  understood,  both  qualitatively  as  well  as  quantitatively.  For  any 
given  structure  or  component  one  approach  to  this  is  to  calculate  the  magni¬ 
tude  of  thermal  stresses  from  the  known  heat  transfer  conditions  and  pertinent 
material  property  data.  Alternately,  the  thermal  shock  response  can  be 
established  experimentally. 

One  such  method,  very  popular  in  view  of  its  simplicity  consists  of 
quenching  appropriate  specimens  at  some  predetermined  initial  temperature 
into  a  fluid  medium  at  lower  or  higher  temperature.  Such  fluid  media  can  con¬ 
sist  of  water,  liquid  metal,  eutectic  salt  mixtures,  oil,  fluidized  powders 

3-9 

and  stationary  or  flowing  air.  For  quantitatively  meaningful  evaluation  of 
the  results  of  such  quenching  test,  it  is  imperative  that  the  heat  transfer 
characteristics  of  these  fluid  media  are  well  understood.  For  any  given  medium, 
literature  data  suggest  considerable  variation  in  heat  transfer  coefficient. 

In  this  respect,  water  as  a  quenching  medium  shows  a  large  variation  in  heat 
transfer,  due  to  nucleate  boiling  and  film- formation,  which  is  anticipated 
to  depend  on  the  temperatures  of  both  the  specimen  as  well  as  the  water  it¬ 
self."*’  Silicon  oils  as  thermal  quenching  media,  also  show  considerable 

variation  in  heat  transfer  behavior,  attributable  in  large  part  to  the 
temperature  dependence  of  the  viscosity.6  Literature  data  for  the  heat  transfer 
characteristics  of  fluidized  beds  show  an  order  of  magnitude  variation, 
presumably  due  to  differences  in  the  flow  rate  of  the  fluidizing  gas,  the 
particle  size  of  the  fluidized  powder  and  other  variables. 


1 


In  the  analysis  of  quenching  data,  it  is  generally  assumed  that  the 

transfer  of  heat  occurs  in  a  Newtonian  manner  such  that  the  heat  flux  is 

directly  proportional  to  the  instantansous  difference  in  temperature  of 

the  specimen  surface  and  the  fluid  medium.  It  is  also  assumed  that  the 

transfer  of  heat  is  spatially  uniform  over  the  specimen  surface.  Both  these 

assumptions  permit  the  use  of  analytical  solutions  for  the  thermal  stresses 

13  14  5 

in  circular  solid  or  hollow  cylinders,  ’  flat  plates  or  spheres  sub¬ 
jected  to  instantaneous  changes  in  ambient  temperature  with  Newtonian  con¬ 
vective  heat  transfer. 

The  transfer  of  heat,  however,  is  Newtonian  only  in  the  case  of  forced 
convection,  in  which  the  fluid  medium  is  passed  over  the  ceramic  specimen 
or  during  the  movement  of  a  specimen  through  a  stationary  fluid.  This 
latter  condition  occurs  during  the  type  of  quenching  in  which  the  specimen 
is  permitted  to  fall  through  the  fluid  by  gravity.  On  the  other  hand,  if 
a  ceramic  specimen  is  held  stationary  in  a  stationary  fluid,  as  frequently 
is  the  case  in  automated  thermal  fatigue  studies,  the  specimen  is  subjected 
to  natural  convection,  with  the  exception  of  the  brief  period  during  which 
the  specimen  is  inserted  in  the  fluid,  where  forced  convection  is  the 
primary  mode  of  heat  transfer. 

The  heat  transfer  coefficient  under  conditions  of  natural -convection  is 

a  function  of  the  instantaneous  difference  in  temperature  of  the  specimen 
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surface  and  the  fluid  medium.  This  implies  that  the  heat  transfer  coeffic¬ 
ient  is  a  maximum  at  the  onset  of  the  thermal  quench  and  decreases  as  the 
specimen  approaches  thermal  equilibrium.  For  this  reason,  the  magnitude 
of  the  thermal  stresses  for  natural  and  forced  convective  heat  transfer  are 
not  expected  to  be  equal.  The  validity  of  this  latter  conclusion  was 
demonstrated  in  an  earlier  study^  by  approximate  analytical  expressions 
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valid  for  very  low  values  of  the  Biot  number.  Ceramic  specimens  of  large 
size  and  low  values  of  thermal  diffusivity  (such  as  glasses)  can  exhibit 
time  periods  to  maximum  thermal  stress  which  can  be  quite  large  compared 
to  the  time  period  required  for  insertion  of  the  specimen  into  the  fluid 
medium.  In  this  case,  for  a  specimen  held  stationary  in  a  stationary 
fluid, the  transfer  of  heat  will  occur  primarily  by  natural  convection. 

A  variable  which  further  complicates  analyses  of  thermal  quenching 
results  is  that  the  heat  transfer  coefficient  for  both  natural  and  forced 
convection  can  be  spatially  highly  non-uniform.^  This  effect  is  well 
known  in  the  heat  transfer  literature.  However,  its  role  during  the 
thermal  shock  testing  of  ceramic  materials  by  the  quenching  method  does 
not  appear  to  have  been  investigated. 

For  the  proper  analysis  of  the  results  of  quenching  studies,  quanti¬ 
tative  data  for  the  transient  thermal  stresses  for  natural  convection  and  for  forced 
convection  with  a  spatially  variable  heat  transfer  coefficient  are  required.  No 

such  information  appears  to  be  available  in  the  literature,  probably  be¬ 
cause  it  is  not  easily  obtained  by  either  analytical  or  experimental 
means.  The  purpose  of  this  study  was  to  obtain  the  required  information 
by  numerical  means. 

II.  GEOMETRY,  BOUNDARY  CONDITIONS,  NUMERICAL  PROCEDURES  AND 

EXPERIMENT 

A  solid  circular  eyi inder  was  chosen  as  the  geometry  appropriate 
for  this  study,  since  this  specimen  shape  frequently  is  chosen  for  thermal 
quenching  experiments.  For  simplicity,  the  cylinder  was  assumed  to  be  of 
infinite  length  which  eliminates  the  complexity  of  possible  end-effects. 

In  practice,  thermal  fracture  near  the  center  of  finite-sized  specimens 
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All  relevant  material 
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is  promoted  by  thermal  insulation  of  the  end. 
properties  of  the  cylinder  such  as  the  thermal  conductivity,  density, 
specific  heat,  coefficient  of  thermal  expansion.  Young's  modulus  of 
elasticity  and  Poisson's  ratio  were  assumed  to  be  independent  of  tem¬ 
perature  and  position.  In  order  to  obtain  information  on  the  separate 
effects  of  the  aT-  and  position-dependent  h,  these  effects  were  not 
superposed  in  the  calculations  of  the  temperatures  and  thermal  stress. 

This  was  justified  in  retrospect,  as  the  effect  of  the  spatial  variation 
in  heat  transfer  coefficient  on  the  magnitude  of  thermal  stress  was 
found  to  be  relatively  small.  Furthermore,  superposition  of  these  two 
effects  would  have  rendered  the  calculations  quite  lengthy  and  prohibitive¬ 
ly  expensive.  The  natural  and  forced  convection  modes  of  heat  transfer 
during  thermal  quench  are  both  expected  to  be  accompanied  by  transient  effects. 
Equations  in  the  literature  for  these  two  modes  of  heat  transfer  are  generally 
applicable  to  steady  state  conditions.  For  the  purpose  of  present  analysis,  it 
is  assumed  that  during  the  thermal  quench,  the  steady  state  conditions  prevail. 
For  transient  effects,  the  reader  is  referred  to  studies  such  as  those  cited 
by  Lin.19 

For  natural  convection  the  dependence  of  the  average  heat  transfer 
coefficient  on  the  instantaneous  temperature  difference  aT^  between  the 
surface  of  the  cylinder  and  the  fluid  was  chosen  to  be  of  the  form:1*1 


h  =  C(aT1-)1/4 

(1) 

where 

C  =  C'(gYfpfcf/vf)1/4d'1/4Kf'3/4 

(2) 

In  which  C'  is  a  constant  which  depends  on  the  values  of  the  Grashof  and  Prandtl 
numbers,  g  is  the  gravitational  acceleration,  y^  is  the  coefficient  of  volume  thermal 
expansion  of  the  fluid,  p^-  is  the  density  of  the  fluid,  c^  is  the  specific 
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heat  of  the  fluid,  vf  is  the  fluid  viscosity,  d  is  the  diameter  of  the  cylinder 
and  k.p  is  the  thermal  conductivity  of  the  fluid. 

The  spatial  variation  of  the  heat  transfer  coefficient  h  depends  on  the 
experimental  conditions.17  For  this  study  the  spatial  variation  of  h 

was  taken  as  described  by  the  experimental  data  of  Krall  and  Eckert17 

for  a  cylinder  with  its  axis  oriented  perpendicular  to  the  direction 

of  fluid  flow  for  the  value  of  the  Reynolds  number,  R0  =  4640.  These  data  are 

included  in  Fig.  6,  expressed  in  terms  of  the  angular  position  around  the 

surface,  relative  to  the  average  value  of  h.  Since  the  heat  transfer  is 

symmetric  about  the  diameter  parallel  to  the  fluid  flow,  the  data  for  h 

are  given  over  180°  only.  The  spatial  variation  of  the  heat  transfer  coefficient 

as  described  in  Fig.  6,  would  correspond  to  a  quench  specimen  falling  through 

the  quenching  medium.  Axial  variation  of  heat  transfer  was  assumed  to  be 

negligible. 

The  temperatures  and  thermal  stresses  were  calculated  by  the  finite  element 
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method  with  the  aid  of  a  computer  program  used  for  a  number  of  previous  studies.  ’ 
The  program  was  modified  to  incorporate  the  boundary  conditions  for  the  heat 
transfer  described  above.  The  model  consisted  of  150  triangular  elements 
forming  a  half  sol  id  cyl  inder  with  the  smallest  elements  being  near  the  cylinder 
surface  where  the  temperature  and  stress  gradients  are  highest.  Thirty  of 
these  elements  lay  along  the  half  cylinder  circumference  yielding  good  spatial 
resolution  of  the  heat  transfer  coefficient.  The  coefficient's  distribution 
given  in  Fig.  6  was  interpolated  to  give  the  values  on  the  surface  elements 
of  the  model.  The  computed  temperatures  at  the  model  node  points  provided 
the  thermal  loading  on  the  stress  model.  Generalized  plane  strain  conditions 
were  assumed  throughout.  The  rod  used  in  quenching  tests  generally  has  its 
ends  unconstrained.  For  this  reason,  no  net  axial  load  or  flexural  bending 
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moment  due  to  the  axial  stress  distribution  was  permitted.  The  axial 
stresses  were  first  computed  from  the  in-plane  stress  components  assuming  the 
cylinder  ends  to  be  fully  constrained  (plane  strain  conditions).  The  axial 
load  and  bending  moment  due  to  the  constrained  axial  stress  distribution 
were  then  calculated.  The  net  thermal  stresses  in  the  unconstrained  cylinder 
were  then  obtained  by  subtraction  of  the  latter  values  from  the  stress  values 
for  the  fully  constrained  cylinder. 

For  the  calculations,  the  numerical  values  of  the  materials  and  fluid 

properties  were  chosen  to  be  typical  of  a  polycrystalline  alumina  or  a  soda- 
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lime-glass  quenched  into  fluid  media  over  a  range  of  quenching  temperature 
difference  required  to  induce  thermal  stress  fracture.6  Variation  of  the 
radius  and  thermal  conductivity  of  the  cyclinder  permitted  calculations  of  the 
thermal  stresses  over  a  range  of  values  of  the  Biot  number,  6  =  Rh/kc  where 
R  is  the  cylinder  radius  and  kc  is  the  thermal  conductivity  of  the  cylinder. 

In  order  to  facilitate  comparison,  the  stress  results  for  the  &T- 
dependent  h  (eq.  1)  were  reported  in  terms  of  the  Biot  number  (eQ)  calculated 
using  the  value  of  the  heat  transfer  coefficient  (hQ)  at  the  beginning  of  the 
thermal  quench,  i.e.,  the  value  of  h  in  eq.  1  corresponds  to  the  value  of 
AT  equal  to  the  initial  temperature  difference  over  which  the  cylinder  is 
being  quenched.  For  the  spatially  dependent  h,  the  data  were  reported  in 
terms  of  the  Biot  number,  calculated  using  the  average  value  of  the  heat 
transfer  coefficient. 

Furthermore,  for  convenience,  the  transient  temperature  and  stresses  were 
expressed  in  terms  of  the  non-dimensional  stress: 

o  =  o(l-v)/aEAT  (3) 


the  non-dimensional  temperatures: 


(4a) 


T*  =  (T  -  Tf)/(To  -  Tf )  (cooling) 

and  * 

T  =  (T  -  TQ)/(Tf  -  To)  (heating)  (4b) 

and  non-dimensional  time: 

*  ? 

t  =  ict/R 

where  a  is  the  stress,  v  is  Poisson's  ratio,  a  is  the  coefficient  of  thermal 
expansion,  E  is  Young's  modulus,  *  is  the  thermal  di ffusi vi ty,  t  is  the  time, 

Tq  is  the  initial  specimen  temperature,  is  the  temperature  of  the  fluid 
medium  and  aT  =  TQ  -  Tf. 

The  validity  of  the  dependence  of  the  heat  transfer  coefficient  on  (aT.)^4  in 
natural  convection  was  established  experimentally  by  measuring  the  temperature 
response  of  an  appropriate  specimen  consisting  of  a  right  circular  cylinder  5 
cm  in  height  by  5  cm  in  diameter.  The  cylinder  was  composed  of  an  industrial 
polycrystall ine  alumina  for  which  the  thermal  diffusivity  and  conductivity 
were  measured  by  the  laser-flash  method.  The  pre-heated  cylinder  initially 
at  thermal  equilibrium  at  approximately  40°C  was  subjected  to  sudden  cooling 
by  immersion  into  a  silicon  oil  over  a  temperature  range  of  20°C.  The  tempera¬ 
ture  at  the  exact  center  of  the  cylinder  was  measured  by  a  thermocouple  inserted 
into  a  hole  drilled  to  the  center  along  the  cylinder  axis.  The  relatively  low 
thermal  diffusivity,  the  large  specimen  size  and  the  low  heat  transfer  coefficient 
for  aT  =  20°C  assured  that  time  period  of  less  than  1  sec  required  for  insertion 
of  the  cylinder  into  the  oil  was  very  small  compared  to  the  total  time  period 
of  approximately  30  min.  over  which  the  temperature  was  recorded.  The  data 
obtained  were  converted  to  non-dimensional  temperature  and  plotted  against  the 
non-dimensional  time  calculated  with  the  aid  of  the  value  of  thermal  diffusivity 
and  cylinder  dimensions. 

The  transient  temperature  response  at  the  center  of  the  cylinder  was  cal- 

22 

culated  from  the  analytical  solutions  of  Fautz.  For  the  AT-dependent  heat 
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transfer,  the  latter  solutions  were  modified  to  permit  the  calculations  of 
the  transient  temperatures  by  iterative  means.  Calculations  were  made  only 
for  value1'  of  the  Biot  number  near  the  value  corresponding  the  dimensions 
ird  thermal  conductivity  of  the  cylinder  and  the  heat  transfer  characteristics 
of  the  silicon  oil  for  the  particlar  conditions  of  the  experiment. 

III.  RESULTS  AND  DISCUSSION 

A.  Natural  Convection  With  AT-Dependent  h 

figures  la  and  lb  show  the  time-dependence  of  the  temperatures  at  the 

surface  and  center  of  the  cylinder,  respectively,  for  both  natural  (h  =  C(aT.)^) 

and  forced  (h  =  constant)  convective  cooling  for  a  range  of  values  of  the  Biot 

number.  The  transient  temperatures  for  convective  heating  can  be  obtained  by 

substracting  the  data  shown  in  Fig.  1  from  unity.  The  values  for  constant  h 
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agree  with  the  analytical  solution  of  Schneider.  Comparison  of  the  data  for 
the  two  modes  of  heat  transfer  indicates  that  any  given  temperature  for 
natural  convection  lags  behind  the  corresponding  value  of  temperature  for 
forced  convection.  This  is  in  accordance  with  general  expectations,  because 
h  decreases  with  decreasing  aT^  for  natural  convection,  which  is  not 

the  case  for  forced  convection. 

Thermal  stress  failure  in  brittle  ceramic  materials  subjected  to 

quenching  occurs  usually  in  tension.  For  this  reason,  the  maximum  values 

of  the  tensile  thermal  stresses  are  of  primary  interest.  Figures  2a  and  2b 

compare  the  tensile  thermal  stresses  in  the  surface  of  the  cylinder  cooled 

by  natural  and  forced  convection  for  a  range  of  values  of  the  Biot  number  as 

a  function  of  time.  On  the  surface  of  the  cylinder  the  thermal  stress  state 

is  uniform  biaxial,  so  that  the  data  shown  in  Figs.  2a  and  2b  apply  both  to 
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the  axial  and  tangential  stresses.  For  forced  convection  (h  =  constant) 
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the  data  show  excellent  agreement  with  the  solutions  of  Jaeger  .  Comparison 
shows  that  for  any  given  value  of  Biot  number  the  stresses  for  natural 
convection  are  less  than  those  for  forced  convection.  The  time  period  required 
for  the  stresses  to  reach  their  maximum  value  is  also  less  for  natural  convection 
than  for  forced  convection. 

Figures  3a  and  3b  compare  the  transient  tensile  thermal  stresses  at  the  center 
of  the  cylinder  for  natural  and  forced  convective  heating.  These  stresses  are 
a  maximum  in  the  axial  direction  and  are  twice  the  value  of  the  stresses  in  the 
radial  and  tangential  direction.  The  stresses  for  natural  convection  are 
less  than  those  for  forced  convection,  as  was  found  for  cooling.  The 
time  required  for  the  stresses  to  reach  their  maximum  value  for  natural  convec¬ 
tion  are  somewhat  greater  than  those  for  forced  convection. 

Figure  4  summarizes  the  values  of  maximum  tensile  thermal  stress.  The 
relative  difference  in  the  magnitude  of  the  stresses  for  natural  and  forced 
convection  is  greater  for  cooling  than  for  heating. 

Figure  5  compares  the  experimental  data  for  the  temperature  response  at 

the  center  of  the  finite  cylinder  with  the  calculated  values  for  h  =  constant 

and  h  =  C(aT..)^.  For  the  latter  case  the  calculated  temperature  response 

lags  behind  the  corresponding  data  for  the  former  case.  The  experimental 

data  over  most  of  the  temperature  range  are  parallel  to  the  curve  for 

natural  convection.  This  implies  that  the  choice  of  h  =  C(aT)^,  also 
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backed  by  extensive  literature  data,  appears  to  be  appropriate  to  describe 
the  heat  transfer  coefficient  for  quenching  in  which  a  specimen  is  held 
stationary  in  a  stationary  fluid. 

B.  Spatially  Varying  Heat  Transfer  Coefficient 

Figure  6  shows  the  angular  distribution  of  the  relative  heat  transfer 
coefficient,  temperature,  axial  and  tangential  stresses  on  the  surface  of 


9 


the  cylinder  at  the  instant  of  maximum  tensile  thermal  stress  for  Biot  number 


3=1.5  and  an  average  value  of  the  heat  transfer  coefficient,  h  =  0.13  cal.°C'^. 

-1-2  6 
s  .cm  typical  of  values  encountered  in  quenching  studies.  The  data  indicate 

that  the  angular  relative  variation  of  the  temperature  is  less  than  the  rela¬ 
tive  variation  in  heat  transfer  coefficient.  This  effect  arises  because  heat 
flow  within  the  cylinder  can  occur  in  the  tangential  as  well  as  the  radial 
direction.  The  axial  and  tangential  stresses  are  no  longer  equal  as  is  the 
case  for  the  spatially  invariant  heat  transfer  coefficient.  An  additional 
contributing  factor  to  this  observation  is  that  because  of  the  circumferential 
non-uniform  temperature  distribution,  the  cylinder  undergoes  bending  de¬ 
formation  with  a  corresponding  reduction  in  the  thermal  stresses. 

Figure  7  compares  the  time-dependence  of  the  value  of  the  maximum  axial 
tensile  stress  for  spatially  uniform  and  non-uniform  convective  cooling  for  the 
values  of  Biot  number  and  average  heat  transfer  coefficient  for  the  data  shown 
in  Fig.  6.  The  stresses  for  the  non-uniform  heat  transfer  exceeds  the  correspond 
ing  values  for  uniform  h.  Also,  for  the  non-uniform  heat  transfer  the  maximum 
value  of  stress  occurs  at  a  somewhat  shorter  value  of  time  than  for  the  uniform 
h. 

Figure  8  shows  the  values  for  the  maximum  tensile  (axial)  thermal  stress 

as  a  function  of  Biot  number  for  the  spatially  uniform  and  variable  h.  The 

relative  difference  in  the  stresses  for  the  two  heat  transfer  conditions  is 
greater  for  cooling  than  for  heating.  In  fact,  on  heating  the  stresses  for 
the  variable  h  are  less  than  for  the  uniform  h.  Figure  9  gives  the  actual 
percentage  differences  for  the  data  shown  in  Fig.  8. 

C.  Discussion 

The  results  presented  above  indicate  that  both  the  spatial  variation  as 

well  as  the  AT-dependence  of  the  heat  transfer  coefficient  have  an  effect  on  the 
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magnitude  of  the  maximum  tensile  stresses  encountered  during  quenching.  For 
this  reason,  they  should  have  relevance  to  the  interpretation  of  test  results 
as  well  as  the  design  of  quenching  experiments. 

It  is  encouraging  to  note  that  the  spatial  non-uniformity  of  the  heat  trans¬ 
fer  coefficient  for  the  specific  numerical  data  chosen,  has  a  relatively  small 
effect  on  the  magnitude  of  the  maximum  tensile  thermal  stress.  In  the  case 
of  heating,  in  fact,  this  effect  to  a  first  approximation  is  negligible.  Even 
on  cooling  the  maximum  difference  in  stresses  for  the  spatially  uniform  and  non- 
uniform  h  is  only  about  17%,  inspite  of  the  factor  of  three  variation  in  the 
heat  transfer  coefficient  as  pointed  out  earlier.  At  least  qualitatively, 
the  stresses  are  affected  only  slightly,  possibly  because  of  the  combined 
effects  of  non-radial  heat  flow  and  the  accommodation  of  the  differential 
thermal  strains  by  bending. 

The  numerical  differences  in  the  maximum  thermal  stresses  due  to  the  non¬ 
uniformity  of  the  heat  transfer  coefficient  are  relevant  to  prediction  or 
interpretation  of  the  results  of  thermal  quenching  experiments.  Such 
experiments  involving  single  thermal  cycles  concentrate  on  determining  the 
critical  temperature  difference  (aTc)  required  to  initiate  thermal  stress 
failure.  Calculations  of  aT  from  data  for  physical  properties  and  spatially 
indepedent  heat  transfer  coefficient  will  result  in  an  over-estimate  of 
aTc  on  cooling  and  an  under-estimate,  albeit  small,  for  heating.  On  the 
other  hand,  if  experimental  values  of  AT  are  used  to  evaluate  the  tensile 
failure  stress,  ignoring  the  spatial  dependence  of  the  heat  transfer 
coefficient  will  result  in  an  underestimate  on  cooling  and  an  overestimate 
for  heating  at  least  for  the  data  of  the  present  study. 

However,  the  statistical  nature  of  brittle  fracture  is  expected  to 
reduce  the  above  over-  and  underestimates.  In  the  tensile  failure  of 


brittle  materials  each  stressed  element  makes  a  contribution  to  the  failure 

probability.  For  this  reason,  estimates  of  the  effect  of  non-uniform  h  on 

thermal  stress  failure  should  be  based  on  statistical  theories  such  as  the 
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Weibull  theory. 

It  should  be  noted  that  the  stress  values  shown  in  Fig.  7  and  8  are 
strictly  valid  only  for  the  specific  numerical  example  chosen.  Different  re¬ 
lative  effects  may  be  encountered  for  other  geometries  and/or  heat  transfer 
conditions.  For  this  reason,  each  case  must  be  examined  individually. 

Regardless  of  the  numerical  differences  involved,  the  existence  of  an 
effect  Of  non-uniform  h  on  magnitude  of  thermal  stress  suggests  care  in  the 
performance  of  quenching  studies.  In  the  example  considered  above,  the  tangen¬ 
tial  non-uniformity  of  h  occurs  because  the  relative  motion  of  the  fluid  flow 
Occurs  perpendicularly  to  the  axis  of  the  cylinder.  Such  tangential  non-uniformity 
does  not  exist  if  the  fluid  flow  occurs  parallel  to  the  cylinder  axis.  For 
this  reason,  specimen  orientation  relative  to  the  fluid  is  expected  to  influence 
the  results  of  quenching  studies.  Unspecified  or  random  orientation  can  handicap 
data  comparison  and  contribute  to  data  scatter.  Identical  specimen  orientation, 
either  perpendicular  or  parallel  to  the  fluid  flow  is  preferred.  In  the  latter 
case  axial  variation  in  heat  transfer  coefficient  may  occur  which  may  move  the 
position  of  maximum  tensile  stress  and  resulting  fracture  (if  any)  towards  the 
leading  end  of  a  cylinder  of  finite  length.  Such  an  effect,  however,  was  not 

observed  during  a  recent  study  of  the  heat  transfer  characteristics  of  a  fluidized 
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bed  for  thermal  quenching  purposes. 

The  large  relative  difference  between  the  maximum  tensile  stresses  for 
forced  and  natural  convection  for  the  same  initial  value  of  the  Biot  number  also 
has  significant  implications  for  the  analysis  of  thermal  quenching  experiments. 

The  range  of  the  applicability  of  the  thermal  stresses  calculated  for  these 
two  nodes  of  heat  transfer  requires  examination  in  some  detail.  Clearly,  in 
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the  case  of  a  specimen  moving  relative  to  the  fluid,  the  stress  data  for 
forced  convection  (h  =  constant)  will  apply.  For  specimens  held  stationary 
in  a  situationary  fluid,  the  question  of  which  set  of  stress  data  apply 
is  more  complex.  Quenching  into  a  stationary  fluid  requires  insertion  of 
the  specimen  into  the  fluid.  During  this  period,  because  of  the  relative 
motion  between  the  specimen  and  fluid,  forced  convective  conditions  will 
prevail.  Once  the  specimen  has  reached  its  stationary  position,  after 
a  transient  during  which  any  fluid  turbulence  is  damped  out,  heat  transfer 
will  occur  by  natural  convection.  For  these  reasons  the  dominant  mode  of 
heat  transfer  is  expected  to  critically  depend  on  the  time  period  required  for 
specimen  insertion  relative  to  the  time  required  for  the  stresses  to  reach 
their  maximum  value.  Two  limiting  cases  can  be  defined.  In  the  first  case, 
the  time  of  maximum  stress  is  less  than  the  time  period  required  for  specimen 
insertion.  In  this  case,  the  specimen  is  subjected  to  forced  convective 
heat  transfer  with  constant  h  only.  In  the  second  case,  the  time  to 
maximum  stress  greatly  exceeds  the  time  period  required  for  specimen  insertion, 
such  that  the  heat  transferred  during  insertion  becomes  negligible  compared  to 
the  total  heat  transfer  while  the  specimen  is  held  stationary.  In  this  case, 
the  specimen  is  primarily  subjected  to  natural  convection  mode  of  heat  transfer. 

These  two  limiting  cases  can  be  illustrated  by  a  simple  numerical 
example.  Let  it  be  assumed  that  it  takes  0.1  second  for  a  quench  specimen  to 
move  through  the  fluid  to  reach  its  stationary  position.  For  a  Biot  number, 

8=2,  the  non-dimensional  time,  *ct/R  required  for  the  maximum  value  of  stress 
is  of  the  order  of  0.08.  For  thermal  stress  failure  to  be  caused  by  forced 
convective  heat  transfer  during  insertion  of  the  specimen  into  the  fluid, 
this  value  of  t*  should  be  such  that  the  real  time  of  maximum  stress  is  equal 
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or  less  than  0.1  sec.  Taking  a  value  for  the  thermal  diffusivity,  k=0.03  cm.  s 


(typical  for  a  polycrystalline  alumina),  forced  convection  heat  transfer  will 
prevail  for  a  cylinder  radius  R  <  0.2  cm.  Similar  values  of  specimen  sizes 
can  be  calculated  for  other  materials  with  ranges  of  thermal  diffusivity 
and  Biot  number  encountered  for  typical  ceramic  materials.  A  specimen  radius 
of  R  <  0.2  cm  is  of  the  order  common  to  the  vast  majority  of  laboratory  studies 
of  the  quenching  behavior  of  many  ceramic  materials.  Obviously,  then  the 
analysis  of  the  results  of  such  studies  should  be  based  on  the  thermal  stresses 
calculated  for  a  constant  heat  transfer  coefficient. 

Thermal  stress  failure  by  naturally  convective  heat  transfer  will 
occur  when  the  time  of  maximum  stress  greatly  exceeds  the  time  period  required 
for  specimen  insertion,  so  that  to  a  first  approximation  the  heat  transferred 
during  this  latter  period  is  negligible  relative  to  the  total  heat  transferred 
by  natural  convection.  For  an  estimate  of  the  minimum  specimen  size  required 
for  this  latter  condition  to  prevail  it  may  be  assumed  that  the  time  of  maximum 
stress  should  be  greater  than  10  sec.  For  the  same  values  for  the  Biot  number, 
time  to  maximum  stress  and  thermal  diffusivity  considered  above,  this  condition 
is  met  for  a  cylinder  radius  R  1!  2  cm.  This  value  is  equal  to  the  radius  of 
the  alumina  specimen,  which  as  shown  by  the  data  in  Fig.  5,  demonstrated  that 
natural  convection  was  the  primary  mode  of  heat  transfer.  Radii  larger  than 
2  cms  are  encountered  in  actual  ceramic  components  for  industrial  or  other 
applications  subjected  to  thermal  shock  during  quenching  or  during  removal 
from  a  kiln,  furnace  etc.,  for  cooling  in  room  air.  The  results 
of  thermal  fracture  under  such  conditions  should  be  based  on  thermal  stresses 
calculated  from  heat  transfer  coefficients  for  natural  convection.  Analysis 
of  such  results  on  the  basis  of  forced  convection  will  result  in  an  over¬ 
estimate  of  the  thermal  stresses  at  which  failure  occurred.  On  the  other  hand, 
of  course,  the  assumption  of  heat  transfer  by  forced  convection  when  in  fact, 


natural  convection  heat  transfer  prevails,  will  yield  a  conservative 
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estimate  of  thermal  stress  resistance.  As  pointed  out  elsewhere  ,  major 
differences  in  thermal  fatigue  life  can  be  found  between  forced  and  natural  con¬ 
vective  heat  transfer,  not  only  due  to  the  differences  in  stress  values 
but  also  due  to  the  difference  in  the  duration  of  each  thermal  stress  cycle.  For 
specimen  sizes  intermediate  to  values  below  and  above  which  heat  transfer 
occurs  primarily  by  forced  and  natural  convection,  resp. ,  heat  transfer 
will  occur  by  both  mechanisms.  Thermal  stress  analysis,  in  that  case,  should 
be  based  on  the  superposition  of  the  thermal  stress  by  an  initial  pulse  of 
forced  convection,  followed  by  natural  convection  thereafter. 

In  general,  the  results  of  this  study  illustrate  that  a  quantitatively 
meaningful  analysis  of  experimental  results  of  thermal  stress  fracture  or 
their  prediction,  requires  a  quantitative  understanding  of  the  heat  transfer 
environment. 
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Time  dependence  of  temperature  at  surface  of  solid  cir¬ 
cular  cylinder  subjected  to  change  in  ambient  temperature 
for  natural  (h  =  C(aT)~)  and  forced  (h  -  constant)  con¬ 
vection  for  a  range  of  values  of  Biot  number. 


Fig.  lb.  Time  dependence  of  temperature  at  the  center  of  .  cir 
cular  cylinder  subjected  to  change  in  ambient  temperature 
for  natural  convection  (h  =  C(,'T)'4)  and  forced  convection 
(h  =  constant)  for  a  range  of  values  of  the  Riot  number. 
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Fig.  6.  Angular  distribution  of  relative  heat  transfer  coefficient, 
temperature,  axial  and  tangential  stresses  on  the  surface 
of  a  solid  circular  cylinder  subjected  to  change  in 
ambient  temperature  at  the  instant  of  time  of  the  maximum 
axial  tensile  stress:  h_..  =  0. 1  3  cal .  ‘  c . -1  s . -1  cm"2  and 
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Fig.  7.  Time  dependence  of  the  axial  tensile  stress  on 
the  surface  of  a  solid  circular  cylinder  sub¬ 
jected  to  sudden  cooling  for  spatially  uniform 
and  variable  heat  transfer  coefficient,  h. 
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Relative  difference  in  the  maximum  axial  tensile  stresses 
as  a  function  of  Biot  number  B  in  a  solid  circular  cylinder 
subjected  to  change  in  ambient  temperature  for  spatially 
uniform  and  variable  heat  transfer  coefficient,  h:  sur¬ 
face  stresses  for  cooling  and  center  stresses  for  heating. 
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ABSTRACT 


A  fracture  mechanics  analysis  is  presented  to  study  the  effect 
of  collinear  crack  interaction  on  the  crack  stability  and  propagation 
behavior  in  a  brittle  material  subjected  to  thermal  stress  failure. 

It  is  shown  that  collinear  crack  interaction  may  result  in  crack 
coalescence  and  complete  failure. 


I 


l 


I. 


INTRODUCTION 


Brittle  ceramic  materials  for  use  at  high  temperature  are  highly 

susceptible  to  failure  by  transient  or  steady-state  thermal  stresses.^ 

The  extensive  crack  propagation  which  frequently  accompanies  such  failure 

can  have  a  profound  adverse  effect  on  the  structural  integrity  of  the 

2  3  A 

ceramic  and  corresponding  load-bearing  ability.  ’  ’  For  the  purpose 
of  engineering  design,  material  selection  and  failure  analysis,  it  is 
critical  that  the  variables  which  control  the  nature  and  extent  of 
crack  propagation  in  thermal  stress  fields  is  qualitatively  and  quantitative¬ 
ly  well  understood.  Previous  analyses^’ ^  of  crack  propagation  in  and 
retained  load-bearing  ability  of  brittle  materials  subjected  to  thermal 
stresses,  generally  considered  the  simultaneous  propagation  of  a  number 
of  cracks,  without  considering  possible  crack-interaction  effects.  As 
a  result  of  this  assumption,  these  analyses  predict  that  for  stable  crack 
propagation,  strength  will  decrease  monotonically  with  increasing  severity 
of  thermal  shock,  above  the  level  required  to  initiate  crack  propagation. 

The  same  effect  is  expected  at  levels  of  thermal  shock  for  materials  with 
large  crack  sizes  which  resulted  from  unstable  crack  propagation.  This 
conclusion  implies  that  the  ceramic  always  should  exhibit  some  residual 
load-bearing  ability  regardless  of  the  severity  of  thermal  shock  and 
corresponding  extent  of  crack  propagation. 

Experimentally,  however,  complete  disintegration  of  a  ceramic  structure 
or  laboratory  specimen  on  occasion  is  observed.  One  explanation  for  such 
observation  is  that  due  to  the  combination  of  the  pertinent  material, 
geometric  and  environmental  variables,  the  size  of  the  cracks,  which  re¬ 
sulted  from  thermal  stress  failure,  exceeded  the  dimensions  of  the  ceramic 


structure,  component  or  specimen.  Crack  interaction,  however,  should  also 
be  considered  as  an  explanation  for  complete  fracture.  In  particular 
this  should  be  the  case  for  two  or  more  cracks  located  on  the  same  plane 
of  propagation.  Coalescence  of  such  cracks  may  result  in  a  macro-crack 
exceeding  the  dimensions  of  the  ceramic,  even  when  the  extent  of  propaga¬ 
tion  of  each  crack  can  be  relatively  small.  It  is  the  purpose  of  this 
study  to  analyze  the  effect  of  co-planar  crack  propagation  and  coalescence 
on  the  nature  of  the  corresponding  retained  load-bearing  ability  of  brittle 
ceramics  subjected  to  thermal  stress. 

II.  ANALYSIS 

The  mechanical  model  and  assumptions  for  the  present  study  will  be 
similar  to  those  chosen  for  the  two  previous  studies^’ ^  which  did  not  consider 
crack  interaction.  Specifically,  the  mechanical  model  consists  of  a  thin 
plate  uniaxially  constrained  from  thermal  expansion  on  cooling  over  a 
temperature  range  AT.  A  plate  of  unit  thickness  contains  a  rectangular 
array  of  cracks  of  equal  size  oriented  perpendicularly  to  the  direction  of  con¬ 
straint,  with  crack  dimension  and  spacings  shown  in  Fig.  1.  Interaction  between 
cracks  within  a  column  is  assumed  absent  in  order  to  permit  an  analytical  assess¬ 
ment  of  the  nature  of  crack  propagation.  In  a  previous  study  the  effect  of 
columnar  interaction  on  the  temperature  difference  (AT^)  required  for  the 
onset  of  crack  propagation  was  assessed  on  the  basis  of  the  numerical  results 

of  Delameter  et  al.  It  was  found  that  due  to  a  decrease  in  stress-intensity 

t 

factor  with  decreasing  columnar  crack  spacing  b,  the  value  of  AT^  increased 
over  the  value  for  a  single  crack.  An  estimate  of  the  extent  of  crack  pro¬ 
pagation  including  crack  coalescence  would  require  extrapolation  of  the 


numerical  data  of  Delameter  et  al . ,  with  the  corresponding  uncertainties 
inherent  in  such  a  procedure.  For  this  reason,  an  analytical  approach  is 
preferred.  The  assumption  of  the  absence  of  columnar  crack  interaction 
implies  that  the  results  to  be  obtained  in  this  study  by  analytical  means 
for  AT  as  well  as  the  extent  of  unstable  crack  propagation  represent 
overestimates  and  therefore  are  conservative. 

The  basic  approach  used  for  the  present  analysis  is  identical  to  the 

one  followed  in  previous  studies. From  the  expression  for  the 

potential  energy  in  the  plate  expressed  in  terms  of  its  effective  elastic 

behavior  and  the  fracture  surface  energy  of  the  cracks,  the  condition 

for  the  onset  of  crack  instability  can  be  derived,  following  the  energy 
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balance  of  Griffith.  From  the  total  change  in  potential  energy,  the 
total  extent  of  unstable  crack  propagation  can  be  obtained. 

An  expression  can  be  derived  for  the  effective  Young's  modulus  in  the 
direction  perpendicular  to  the  plane  of  the  cracks,  from  the  solution  of 
Yokobori  and  Ichikawa10  for  the  increase  in  potential  energy  (AW)  per 
crack  for  a  single  row  of  co-planar  interacting  cracks  (for  plane  strain) 
in  a  plate  of  unit  thickness  expressed  by: 


AW 


2(l-v2)d2o2 

ttE 

o 


In  cos  ( 


d 


) 


(1) 


where  v  is  Poisson's  ratio,  a  is  the  stress  and  Eq  is  Young's  modulus  of 
the  plate  in  the  absence  of  cracks. 

For  N  cracks  per  unit  area  of  the  plate  with  N  =  (bd)  1 ,  the  change 
in  potential  energy  AW^  becomes: 

2(l-v2)d2o2N 

AW  - - - -  In  cos  (■— )  (2) 

N  irE  d 
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The  total  elastic  energy  (W^)  in  the  plate  is  the  elastic  energy 
in  the  plate  without  cracks  plus  the  change  in  energy  due  to  the  cracks; 


,,  °2  r-i  2(l-v2)d2N  , 

WP  =  W~  [1 - 7 - ln  cos  (T)] 


from  which  the  effective  Young's  modulus,  E  ^  can  be  obtained: 

„  „  r,  4(l-v2)d2N  ,  /iTaNT_l  ,, 

Eeff  =  Eo[l - 7 -  ln  cos  (T)] 

For  a  temperature  difference  AT  over  which  the  plate  is  cooled,  the 
sum  of  the  strain  energy  and  the  fracture  surface  energy  per  unit  area  of 


plate  is: 


a2 (AT)2E  ...  2\ j2m 

o  r,  4 1 l-v/ ) d^N  ,  ,iras-|-l  ,  „ 

W  =  - x -  Ll  - - r -  cos  (-r-)J  +  4GaN  (5) 


where  a  is  the  coefficient  of  thermal  expansion  (with  a  >  0)  and  G  is 
the  energy  required  to  create  unit  area  of  new  fracture  surface,  which  includes 
all  forms  of  energy  dissiptation  including  plastic  flow,  viscous  deformation, 
acoustic  energy  losses  and  others. 

The  cracks  are  unstable  at  a  value  of  AT  =  AT  when: 

c 

d(W  ) 

sr- 1 0  (6) 

which  for  simultaneous  crack  propagation,  yields: 


AT  =  { [2G/ (l-v2)a2Ed]  cot  (na/d)  }^[l  -  — -  In  cos  (-^r)  ]  (7) 

c  i  D  a 

In  eq.  7  when  a/d  +  0,  ATc  approaches  the  solution  for  the  absence  of 
crack  interaction. 

The  extent  of  unstable  crack  propagation  for  crack  instability  at  AT  =  ATf 
can  be  derived  by  equating  the  total  loss  in  elastic  energy  on  crack 
propagation  from  an  initial  value,  to  a  final  value,  a ^  to  the 


) 


required  change  in  surface  fracture  energy,  which  yields) 


a2 (AT  )2E 
c  o 


,  4(l-v2)d“N1  /  ,  4(l-u2)d-N,  /  f-, 

1  ~  — ln  cos  (— )  1  -  - —  In  cos  (— -) 


=  4GN(ac  -  a  ) 
f  o 


For  extensive  crack  propagation  such  that  >>  aQ  the  fina]  crack 
length  (if  arrested,  to  be  discussed  later)  to  a  good  approximation  can 
be  expressed  as: 


a2 (AT  )2E 
c  o 


4 ( 1-v2) d2N 


In  cos  (-—•-) 


III.  RESULTS  AND  DISCUSSION 


In  Fig.  2,  the  solid  curves  show  the  normalized  values  of  ATc  as  a 
function  of  the  ratio  a/d  for  a  range  of  values  of  b/d.  With  the  ex¬ 
ception  for  the  highest  values  of  a/d,  above  b/d  =  1.0,  AT^  shows  a  U- 
shaped  behavior  qualitatively  similar  to  the  results  found  for  the  absence 
of  crack  interaction.  At  the  low  values  of  a/d  and  high  values  of  b/d,  AT 
is  not  affected  by  crack  interaction  due  to  the  large  inter-planar  and  inter- 
columnar  crack  spacings  compared  to  the  crack  size.  In  this  range  of  values 
of  a/d  and  b/d,  AT^  decreases  with  increasing  crack  length.  The  increase  in 
AT  for  values  of  a/d  larger  than  the  value  for  the  minimum  in  AT  is  because 
the  effective  Young's  modulus  of  the  plate  decreases  faster  with  increasing 
a/d  than  the  corresponding  decrease  in  stress  required  for  fracture  of  the 
plate.  This  latter  quantity  is  controlled  by  the  collinear  crack  inter¬ 
action  only,  whereas  Young's  modulus  is  governed  by  the  total  crack  arrav. 

The  value  of  crack  length  corresponding  to  the  minimum  in  ATc  will  be  denoted 


An  interesting  effect  occurs  at  the  higher  values  of  a/d.  For  a 
given  value  of  b/d,  AT_  exhibits  a  maximum  at  a  crack  length  denoted  as 

3max‘  F°r  2  >  3max’  ATc  decreases  rapidly  to  approach  zero  as  a/d  '  0.5. 

This  latter  effect  corresponds  to  crack  coalescence  and  the  formation  of 

a  macrocrack  and  complete  disintegration  of  the  material.  The  maximum  in 

ATc  at  the  higher  values  of  a/d,  contrasts  with  the  findings  of  the  previous 

studies  which  indicated  that  without  crack  interaction  AT^  would  increase 

monotonically  with  increasing  crack  length.  (Of  course,  as  long  as  the 

crack  size  did  not  exceed  the  dimensions  of  the  material.) 

As  discussed  earlier,  for  values  of  a/d  smaller  than  a  .  /d ,  the 

min 

strain  energy  release  rate  for  crack  instability  at  AT  ^  AT^  exceeds 

the  fracture  surface  energy  required  to  create  the  new  fracture  surfaces. 

This  fact  will  cause  the  crack  to  propagate  in  a  so-called  unstable  or 

dynamic  mode.  In  contrast,  cracks  with  a  .  /d  <  a/d$  a  /d  will  propagate  in 
1  nun  max  r  r 

a  stable  manner  for  AT  i  ATC>  with  the  relationship  between  ATC  and  a/d 
at  any  stage  of  crack  propagation  process  described  by  eq.  7. 

The  dotted  curves  in  Fig.  2  denote  the  maximum  values  of  the  final  crack 
length,  a ^  which  result  from  the  unstable  crack  propagation  of  cracks 
with  an  initial  crack  length,  a  <  a  .  calculated  by  means  of  eq .  9.  For 
purpose  of  clarity.  Fig.  3  shows  the  value  of  AT^  and  a ^  on  an  expanded 
scale  for  0.4  1  a/d  H  0.5.  As  long  as  the  value  of  a^  for  a  given  value  of 
b/d  lies  below  the  corresponding  curve  for  ATc ,  the  crack  will  be  arrested. 
Furthermore,  it  will  be  subcritical  and  for  continued  propagation  will  require 
a  finite  increase  in  AT  over  the  value  of  AT£  originally  required  to 
initiate  crack  propagation  at  a  =  aQ.  Qualitatively,  a  similar  result  was 
predicted  for  the  absence  of  crack  interaction.  However,  due  to  the  co-planar 
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CRITICAL  TEMPERATURE  DIFFERENCE 


crack  interaction  considered  in  the  present  study,  an  interesting  effect 

will  occur  for  values  of  a^  of  sufficient  size  that  the  corresponding  value 

for  ATc  for  continued  propagation  is  smaller  than  the  value  of  AT^  which 

initiated  the  crack  propagation  at  a  =  aQ.  In  Figs.  2  and  3,  this  corresponds 

to  values  of  a/d  equal  to  or  greater  than  the  critical  value  of  a^/d 

(i.e.  a^/d)  at  the  points  of  intersection  between  the  solid  curves  for 

AT  and  the  dotted  curves  for  ac.  If  unstable  crack  extension  is  sufficiently 
c  r 

extensive  that  this  effect  occurs,  crack  arrest  will  not  take  place.  Instead, 
the  crack  will  continue  to  propagate  in  a  stable  manner  until  a/d  =  0.5  when 
all  cracks  coalesce,  to  result  in  complete  fracture.  It  is  interesting  to 
note  that  this  result  differs  from  that  for  a  single  crack  for  which  under 
"fixed  grip"  conditions,  the  stress  intensity  factor,  Kj.  is  controlled  by  com¬ 
pliance  effects  and  decreases  due  to  stress  relaxation  with  increasing  crack 
length.  In  the  present  case  for  multiple  collinear  cracks,  when  the  cracks 
increase  in  length  and  approach  each  other,  the  stress  intensity  factor 
becomes  very  large  due  to  the  existence  of  a  double  singularity.  This  effect 
dominates  over  the  effect  of  increasing  compliance  which  tends  to  reduce  K^. 
The  combined  result  of  the  two  effects  leads  to  crack  coalescence  and  complete 
fracture.  Clearly,  to  avoid  this  mode  of  failure,  the  initial  crack  size 

in  the  material  prior  to  thermal  stress  failure  should  be  larger  than  a 
critical  initial  size,  aQ  which  would  result  in  the  critical  final  size,  a^. 

The  lengths  of  cracks,  a  and  al  as  well  as  the  length  a  .  and  a 

s  c  f  jjin  max 

corresponding  to  minimum  and  maximum  values  of  ATc  (i.e.,  at  d(ATc)/da  =  0) 

are  shown  schematically  in  Fig.  4  for  given  values  of  planar  and  columnar 

crack  spacings.  Figure  4  also  indicates  the  crack  propagation  behavior  after 

AT  has  been  reached  for  different  values  of  initial  crack  length  a  . 
c  o 


U 
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CRITICAL  TEMPERATURE  DIFFERENCE  (ATC) 


3 


Fig.  4 


Diagram  for  crack  stability  showing  ac,  amin,  amax  and  af-  and 
the  crack  propagation  behavior  for  cracks  of  different  initial 
lengths. 
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For  a  <  a  and  a  >  a  ,  Fig.  4  shows  that  the  crack  length  remains 
o  c  o  max  ° 

constant  until  AT  equals  AT^  when  the  crack  becomes  unstable  and  results  in 

a  complete  fracture  for  reasons  explained  before.  The  paths  for  crack 

propagation  are  shown  in  chain  lines.  For  a  <  a  <  a  .  ,  the  crack  length 

remains  constant  until  the  critical  ATc  is  reached  when  the  crack  propagates 

instantaneously  to  a  new  length  and  becomes  subcritical.  At  this  stage,  further 

increase  in  AT  to  some  value  AT^  is  required  to  re-propagate  the  crack.  For 

AT  >  AT^,  the  crack  propagates  stably  until  the  crack  length  equals  amax 

when  the  material  fractures  completely.  For  a  .  <  a  <  a  ,  the  crack  length 

remains  constant  until  AT  equals  ATc.  For  AT  >  ATc,  the  crack  propagates  in 

a  stable  manner  until  a  is  reached  and  a  complete  fracture  occurs. 

max  r 

Figure  5  indicates  schematically  the  corresponding  strength  loss  behavior 

as  a  function  of  AT.  For  a  <  a  and  a  >  a  ,  the  ceramic  will  exhibit 

o  c  o  max’ 

a  total  loss  in  strength.  For  ac  <  aQ  <  amin  strength  loss  is  governed  by 
both  unstable  and  stable  crack  propagation.  For  a  .  <  a  <  a  ,  the  strength 
less  is  governed  by  stable  crack  propagation  only  and  does  not  exhibit  the 
discontinuity  in  strength  encountered  for  unstable  crack  propagation.  Com¬ 
parison  of  the  four  diagrams  in  Fig.  5  with  the  nature  of  strength  loss  be¬ 
havior  predicted  for  the  lack  of  crack  interaction  shows  that  the  effect  of 
such  interaction  can  have  pronounced  effect  on  the  strength  loss  behavior  of 
a  brittle  ceramic  subjected  to  severe  thermal  shock.  It  should  be  noted,  how¬ 
ever,  that  for  a  batch  of  specimens,  flaw  sizes  and  the  degree  of  possible  flaw 
interaction  is  subject  to  statistical  variation.  For  this  reason,  the  crack 
propagation  behavior  predicted  by  the  present  analysis  and  corresponding 
strength  loss  characteristics  may  be  observed  only  for  a  few  specimens  within 
a  given  batch.  In  fact,  an  unusual  data  scatter  in  retained  strength  values 
at  least  in  part  could  be  attributed  to  varying  degrees  of  crack  interaction. 
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STRENGTH 


For  this  reason,  the  experimenter  may  wish  to  base  his  data  interpretation 
on  the  examination  of  each  specimen  tested  rather  than  on  comparison  of 
averaged  data. 


It  is  suggested  here  that  similar  crack  propagation  behavior  is  ex¬ 
pected  also  in  the  formation  of  microcracks  due  to  the  thermal  expansion 
anisotropy  in  polycrystalline  materials  or  composites,  because  the  variables 
which  affect  microcrack  formation  and  thermal  stress  failure  in  the  absence 
of  external  body  forces  are  quantitatively  identical. 


ACKNOWLEDGMENT 

The  present  study  was  conducted  as  part  of  a  research  program  on  the 
thermal  and  thermo-mechanical  beahvior  of  structural  ceramics  supported 
by  the  Office  of  Naval  Research  under  Contract:  N00014-78-C-0A31 . 


REFERENCES 


1.  W.  D.  Kingery,  "Factors  Affecting  Thermal  Stress  Resistance  of 
Ceramic  Materials,"  J.  Am,  Ceram.  Soc . ,  38  (1)  3-15  (1955). 

2.  W.  R.  Morgan,  "Thermal  Shock  Effect  on  the  Transverse  Strength  of 
Clay  Bodies,"  J.  Am.  Ceram.  Soc.,  14  (2)  913-23  (1931). 

3.  D.  P.  H.  Hasselman,  "Strength  Behavior  of  Polycrystalline  Alumina 
Subjected  to  Thermal  Shock,"  J.  Am.  Ceram.  Soc.,  53  (9),  490-95  (1970). 

4.  T.  K.  Gupta,  "Crack  Healing  and  Strengthening  of  Thermally  Shocked 
Alumina,"  J.  Am.  Ceram.  Soc.,  59  (5-6)  259-62  (1976). 

5.  D.  P.  H.  Hasselman,  "Unified  Theory  of  Thermal  Shock  Fracture  Initia¬ 
tion  and  Crack  Propagation  of  Brittle  Ceramics,"  J.  Am.  Ceram.  Soc., 

52  (11)  600-04  (1969). 

6.  D.  P.  H.  Hasselman,  "Thermal  Stress  Crack  Stability  and  Propagation  in 
Severe  Thermal  Environments,"  pp.  89-103  in  Materials  Science  Research, 
Vol.  V,  Ceramics  in  Severe  Environments,  Ed.  by  W.  W.  Kriegel  and  Hayne 
Palmour  III,  Plenum  Press,  New  York  (1971). 

7.  C.  Shih,  J.  P.  Singh  and  D.  P.  H.  Hasselman,  "Effect  of  Crack  Inter¬ 
action  on  the  Fracture  Initiation  and  Crack  Propagation  in  Brittle 
Ceramics  Subjected  to  Severe  Thermal  Shock,"  High  Temperatures- 
High  Pressures  (in  press). 

8.  W.  R.  Delameter,  G.  Herrmann  and  D.  M.  Barnett,  J.  Appl.  Mech.,  Trans. 
ASME  42,  74-80  (1975). 

9.  A.  A.  Griffith,  "Theory  of  Rupture,"  pp.  55-63  in  Proceedings  First 
International  Congress  of  Applied  Mechanics.  Ed.  by  J.  Waltman,  Jr., 
Delft,  Holland  (1924). 

10.  T.  Yokobori  and  M.  Ichikawa,  "Elastic  Solid  with  Infinite  Row  of  Co- 
linear  Cracks  and  the  Fracture  Criterion,"  J.  Phys.  Soc.  Japan,  19, 
2341-42  (1964). 


CHAPTER  IV 


ANALYSIS  OF  THERMAL  FATIGUE  BEHAVIOR 
OF 

BRITTLE  STRUCTURAL  MATERIALS 


by 


J.  P.  Singh,  K.  Niihara  and  D.  P.  H.  Hasselman 

Department  of  Materials  Engineering 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia  24061  USA 


ANALYSIS  OF  THERMAL  FATIGUE  BEHAVIOR  OF  BRITTLE  STRUCTURAL  MATERIALS 


J.  P.  Singh,  K.  Niihara  and  D.  P.  H.  Hasselman 
Department  of  Materials  Engineering,  Virginia 
Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia  24061,  USA 

ABSTRACT 

The  effect  of  the  level  of  maximum  temperature  (Tm  ),  the  tempera¬ 
ture  range  (aT)  and  the  mode  of  convective  heat  transfer  on  the  thermal 
fatigue  resistance  of  brittle  structural  materials  is  analyzed.  Expres¬ 
sions  are  derived  for  the  number  of  thermal  cycles  to  failure  in  terms  of 
the  appropriate  mechanical  and  thermal  properties,  crack  growth  parameter, 
aT  and  T  For  simultaneous  changes  in  T  and  aT  commonly  used  in 
practice,  the  change  in  thermal  fatigue  life  is  governed  by  both  the  thermal 
stress  intensity  exponent  (n)  and  the  activation  energy  (Q)  for  subcritical 
crack  growth,  in  contrast  to  the  results  of  other  studies.  For  constant 
T  but  variable  aT,  thermal  fatigue  life  is  affected  by  n  only,  whereas, 

fTIa  X 

for  constant  aT  but  variable  T  ,  the  value  of  Q  alone  governs  changes  in 

m3  x 

fatigue-life.  Heat  transfer  by  natural  or  forced  convection  will  result 
in  differences  in  thermal  fatigue  resistance. 

Recommendations  are  made  for  the  design  and  analysis  of  thermal  fatigue 
experiments.  Figures-of-merit  for  the  selection  of  materials  with  high 
thermal  fatigue  resistance  are  presented. 


1.  INTRODUCTION 

Brittle  structural  materials  can  exhibit  static  [1,2],  dynamic  [3,4], 
cyclic  [5,6]  as  well  as  thermal  fatigue  [7,8].  These  mechanisms  of  fatigue 
result  from  sub-critical  crack  growth  due  to  stress-corrosion  at  low  tem¬ 
peratures  or  diffusional  mechanisms  at  high  temperatures.  For  the  purpose 
of  reliable  design  of  structures  or  components  made  of  brittle  materials, 
it  is  imperative  that  the  mechanisms  responsible  for  fatigue  are  qualitatively 
as  well  as  quantitatively  well  understood.  An  extensive  literature  indicates 
that  the  "failure-prediction"  of  brittle  materials  subjected  to  static  [9,10], 
dynamic  or  cycl ic  fatigue  [11]  under  isothermal  conditions  appears  well-under¬ 
stood.  The  statistical  nature  of  brittle  fracture  can  also  be  incorporated 
into  predictions  of  fatigue-behavior. 

The  thermal  fatigue  behavior  of  brittle  materials  is  more  complex  for 
the  principal  reasons  that  both  the  stresses  and  temperatures  change  simul¬ 
taneously,  requiring  numerical  techniques  for  the  calculations  of  thermal 
fatigue  behavior.  In  a  number  of  such  studies  [12-14],  reasonable  to 
excellent  agreement  was  found  between  laboratory  data  and  the  predicted 
number  of  thermal  cycles-to-failure.  More  recently,  by  analyzing  experimental 
data  for  thermal  fatigue  behavior  of  a  number  of  brittle  structural  material, 
Kamiya  and  Kamigaito  [15]  obtained  quantitative  information  on  the  behavior 
of  sub-critical  crack  growth  during  thermal  fatigue. 

In  order  to  provide  a  basis  for  the  objectives  of  the  present  study, 
it  should  be  noted  that  for  convective  heat  transfer  used  in  most  thermal 
fatigue  studies  the  magnitude  of  the  transient  thermal  stresses  is  a  func¬ 
tion  of  the  temperature  differences  encountered.  The  rate  of  sub-critical 
crack  growth,  however,  is  a  function  of  the  thermal  stresses  as  well  as  the 
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absolute  temperatures  involved.  These  two  separate  effects  are  critical 
to  the  design  of  a  thermal  fatigue  experiment  as  well  as  to  the  analysis  of 
the  data  obtained. 

In  general,  the  thermal  fatigue  behavior  of  a  material  is  established 
by  measuring  the  number  of  thermal  cycles  (N)  required  to  cause  failure  of 
appropriate  specimens  cycled  over  a  temperature  range  AT,  between  an 
upper  temperature  T^ax  and  a  lower  temperature,  T  Figure  1  schematical ly 
show  the  time-dependence  of  the  temperatures  and  stresses  in  the  surface  of 
a  thermal  fatigue  specimen. 

P  complete  thermal  fatigue  curve  is  established  by  determining  the  number 
of  thermal  cycles  N  for  a  range  of  values  of  AT.  However,  critical  to 
the  results  obtained  is  the  specific  way  in  which  AT  is  varied.  Commonly, 
aT  is  varied  by  lowering  the  upper  temperature  Tmax  while  the  lower  tempera¬ 
ture  T  .  is  kept  constant.  In  this  case  the  magnitude  of  the  stresses 
min 

as  well  as  the  absolute  values  of  temperature  are  varied.  If  AT  were 
changed  by  raising  (or  lowering)  Tmin,  the  stresses  are  varied,  but  the 
absolute  temperatures  at  which  most  of  the  crack  growth  takes  place  are 
kept  constant.  Similarly,  AT  can  be  kept  constant,  by  changing  Tmax  and 
Tmin  by  equal  amount.  Depending  on  the  kinetics  of  crack  growth,  these 
three  different  conditions  are  expected  to  lead  to  differences  in  thermal 
fatigue  behavior  as  measured  by  the  number  of  cycles-to-failure.  A  further 
complexity  is  introduced  by  noting  that  heat  transfer  can  occur  in  a  forced 
or  natural  convective  mode.  In  forced  convection,  heat  transfer  is  Newtonian, 
i.e.,  the  heat  transfer  coefficient  is  independent  of  the  difference  in 
instantaneous  temperature  of  the  surface  of  the  specimen  and  the  tempera¬ 
ture  of  the  fluid  medium  used  for  heating  or  cooling  the  specimen.  In 


contrast,  for  natural  convection,  the  heat  transfer  coefficient  is  a 
non-linear  function  of  the  difference  in  instantaneous  surface  tempera¬ 
ture  and  the  temperature  of  the  fluid  medium  [16].  Under  these  conditions 
the  heat  transfer  coefficient  changes  as  the  specimen  changes  in  tem¬ 
perature  (see  Fig.  1)  towards  thermal  equilibrium.  Detailed  calculation 
of  this  effect  by  the  finite  element  method,  to  be  reported  elsewhere  [17], 
showed  that  for  a  given  material  and  specimen  size  and  identical  values  of 
the  heat  transfer  coefficient  at  the  beginning  of  the  thermal  cycle,  the 
magnitude  of  the  maximum  thermal  stress  under  conditions  of  natural  convec¬ 
tion,  are  significantly  less  than  the  corresponding  values  for  forced  con¬ 
vection. 

It  is  the  purpose  of  this  study  to  present  an  analytical  treatment 
of  the  temperature  and  heat  transfer  variables  in  the  thermal  fatigue  re¬ 
sistance  of  brittle  materials,  in  order  to  establish  guidelines  for  the 
design  of  thermal  fatigue  experiments  and  analysis  of  the  data  in  terms  of 
crack-growth  behavior. 

2.  ANALYSIS 
2.1.  General 

The  analysis  will  concentrate  on  thermal  fatigue  specimens  in  the 
shape  of  long  cylinders  with  circular  cross-section.  For  this  geometry, 
solutions  for  the  thermal  stresses  are  available  [17,18].  The  ends  of 
the  specimens  are  assumed  to  be  thermally  insulated  so  that  fracture  at  the 
ends  is  avoided.  This  assumption  will  permit  the  use  of  the  thermal 
stresses  calculated  for  circular  cylinders  of  infinite  length.  It  is 
further  assumed  that  the  thermal  fatigue  mechanism  operates  only  during 
the  cooling  part  of  the  cycle.  This  latter  assumption  is  reasonable  since 
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under  those  conditions  the  stresses  exceed  those  during  heating  and  are  a 
maximum  in  the  specimen  surface  most  susceptible  to  stress-corrosion  and 
the  introduction  of  flaws  during  specimen  preparation  and  handling.  The 
results  of  the  analysis  can  be  easily  extended  to  other  specimen  geometries 
or  thermal  fatigue  during  heating  or  combinations  of  heating  and  cooling. 
Finally  it  is  assumed  that  between  thermal  cycles  the  specimen  comes  to 
thermal  equilibrium  so  that  overlapping  of  the  transient  thermal  stresses 
of  successive  thermal  cycles  need  not  be  considered. 


2.2.  Theory 

The  rate  of  sub-critical  crack  growth  (a)  under  conditions  of  thermal 
fatigue  will  be  assumed  to  be  identical  to  the  expression  for  isothermal 
conditions  given  by: 

a  =  AK"  exp  (-Q/RT)  (1) 


where  A  is  a  constant,  n  is  the  stress  intensity  exponent,  Kj  is  the  mode 
I  stress-intensity  factor,  Q  is  the  activation  energy,  R  is  the  gas  constan'. 
and  T  is  the  absolute  temperature. 

In  terms  of  the  thermal  stress  (a)  and  crack  size  (a),  Kj  is  defined 
by: 


Kj  =  Ya/T 


(2) 


where  Y  is  a  constant  related  to  the  stress  distribution  and  the  geometry 
of  the  crack.  Implicit  in  the  use  of  Equation  2,  is  the  assumption  that  over 
the  total  range  of  crack  growth,  the  crack  size  remains  sufficiently  small 
compared  to  the  size  of  the  specimen  that  the  compliance  of  the  specimens 
remains  unaffected.  This  will  permit  calculation  based  ori  thermo-elastic 
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theory.  In  this  manner,  the  complexity  of  the  effect  of  the  presence 
of  the  crack  on  the  thermal  stress  field  and  thermal  stress  intensity 
factor,  which  would  require  numerical  methods,  is  avoided. 

For  a  given  value  of  thermal  stress,  failure  will  occur  at  a  critical 
crack  size,  ac  such  that  KI  =  KIcs  the  critical  stress  intensity  factor. 

The  thermal  fatigue  life  is  defined  by  the  number  of  thermal  cycles  required 
to  propagate  the  crack  from  an  initial  crack  size,  aQ  to  the  critical 
crack  size,  ac. 

In  direct  analogy  to  the  time-to-failure  under  isothermal  conditions, 
the  number  of  cycl es-to-fai 1 ure  (N)  in  thermal  fatigue  can  be  expressed: 

K. 

N  =  B  /  1C  KTdKT/ao2(t)V2  (3) 

o  1  1 

where  B  is  a  constant  and  a(t)  is  the  transient  thermal  stress.  Because 
of  the  relative  complexity  of  the  analytical  expressions  for  the  transient 
temperatures  and  stresses  even  for  the  simple  geometry  such  as  a  solid 
circular  cylinder,  the  integration  of  Equation  3  is  most  conveniently  car¬ 
ried  out  by  numerical  methods,  such  as  those  employed  in  a  number  of  pre¬ 
vious  investigations  [12,19,20]. 

However,  in  order  to  obtain  a  convenient  analytical  expression  for  the 
number  of  cycles-to-fail ure,  a  number  of  simplifying,  but  entirely  reason¬ 
able  assumptions  are  made.  Firstly,  as  indicated  by  Fig.  1  any  temperature 
during  the  thermal  cycle  can  be  expressed  as  a  fractional  value  of  the  maxi¬ 
mum  temperature,  T  v  at  the  beginning  of  the  thermal  cycle.  For  this  reason, 
the  effect  of  the  thermally  activated  nature  of  the  crack  growth  on  thermal 
fatigue  life  can  be  expressed  uniquely  in  terms  of  a  factor  exp(Q/RT  ) . 
Secondly  and  similarly,  any  value  of  the  transient  thermal  stress  during 


the  thermal  cycle  can  be  expressed  as  a  function  of  the  maximum  value  of  ther¬ 
mal  stress,  omax  encountered  during  the  thermal  cycle.  Furthermore,  the 
analytical  results  for  the  time-to-failure  [10]  under  isothermal,  constant 
stress  conditions,  indicate  that  unless  the  total  extent  of  crack  growth 
is  small,  the  time-to-failure  is  an  inverse  function  of  the  value  of  the 
stress  intensity  factor  at  the  time  (t=0)  of  the  application  of  the  load. 

For  these  latter  two  reasons,  it  will  be  assumed  that  the  number  of  cycles- 
to-failure,  as  affected  by  the  stress-intensity  factor,  can  be  uniquely 
defined  in  terms  of  an  inverse  function  of  the  maximum  value  of  the  thermal 
stress  intensity  factor  (K^. )  encountered  during  the  first  thermal  cycle. 

In  effect  this  changes  the  lower  limit  of  the  integral  of  Equation  3  from 
zero  to  Kj ^  Furthermore,  it  will  be  assumed  that  the  number  of  thermal 
cycles-to-failure  is  an  inverse  function  of  the  time-duration  over  which 
the  transient  thermal  stress  acts  for  each  cycle.  Analytical  results  [18] 
for  the  thermal  stresses  indicate  that  the  magnitude  of  stress  as  a  function 
of  time  can  be  expressed  in  terms  of  a  non-dimensional  time,  t  =  xt/R2, 
where  k  is  the  thermal  diffusivity,  t  is  the  real  time  and  R  is  the  radius 
of  the  cylinder.  For  any  material  and  specimen  size,  the  time  period  of  each 
thermal  stress  pulse  is  proportional  to  R2/x  making  the  number  of  cycles- 
to-failure  proportional  to  x/R2. 

Substitution  of  Equation  1  into  Equation  3  and  incorporation  of  the 
above  conclusions,  the  number  of  thermal  cycles-to-failure,  in  analogy  to 
the  time-to-failure  under  iso-stress  and  iso-thermal  conditions  can  be 
written  in  the  general  form: 

N  =  C«  «p(Q/RTmax)/ACT^ax^Y2(n-2)KI1(n-2)  (4) 
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where  C  is  a  constant  which,  in  principle,  can  be  obtained  by  numerical 
integration  of  Equation  3  and  K^.  is  defined  by: 


KT .  =  Vo  v'ST 
Ii  max  i 


(5) 


Substitution  of  Equation  5  into  Equation  4,  results  in  the  value 
of  N  expressed  directly  in  terms  of  the  maximum  value  of  the  thermal 
stress 

n-2 

N  -  Ck  «p<Q/KTmx)/B2£xA(n-2>Y"af  2  (6) 


2.3.  Derivation  of  thermal  fatigue-life 
2.3.1.  Forced  convection 

Under  forced  convection,  the  heat  transfer  coefficient  is  independent 
of  the  range  of  temperature  difference  over  which  the  specimen  is  being 
cycled.  Over  the  range  of  the  Biot  number,  1  <  &  <  20,  the  maximum  value 
of  tensile  thermal  stress  in  the  surface  of  a  circular  rod  to  a  very  good 
approximation  can  be  expressed  [18,21] 


-1 

a 

max 


(1-v) 

ctEAT 


(1 .45  +  4.95/6) 


(7) 


where  a  is  the  coefficient  of  thermal  expansion,  E  is  Young's  modulus 
of  elasticity,  v  is  Poisson's  ratio,  aT  is  the  temperature  range  over  which 
the  specimen  is  being  cycled  (aT  =  Tmax  -  Tm-jn*  see  Fig.  1)  and  6  =  Rh/k 
where  R  is  the  cylinder  radius,  h  is  the  heat  transfer  coefficient  and  k 
is  the  thermal  conductivity. 

Substitution  of  Equation  7  into  Equation  6  yields: 


N 


Ck 


Yn(n-2)AR2a.(n'2^/2 


1-v 

nEAT 


1 .45  +  4.95/6 


exp(Q/RTmax) 

(8) 
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In  order  to  establish  the  nature  (slope,  etc.)  of  a  thermal  fatigue 
curve  for  a  given  material  it  is  most  convenient  to  derive  expressions  for 
the  ratio  of  the  values  of  N  for  two  conditions  of  thermal  fatigue.  As 
discussed  earlier,  there  are  three  basic  ways  in  which  a  thermal  fatigue 
experiment  can  be  carried  out:  1,  AT  can  be  varied  by  raising  or  lowering 
T  ,  keeping  T  .  constant;  2,  aT  can  be  varied  by  keeping  T  constant, 
and  raising  or  lowering  T  ..  and  3,  aT  can  be  kept  constant  by  raising 
or  lowerinq  T  and  T  .  by  equal  amounts. 

The  effect  of  these  three  changes  in  temperature  on  thermal  fatigue- 
life  will  be  different. 

For  analytical  purposes,  it  is  most  convenient  to  examine  these 

quantitatively  by  deriving  expressions  for  the  ratio  of  the  number  of 

cycles-to-failure  for  small  differences  in  T  ,  T  .  or  aT,  which  yield 

max  min 

the  following  results. 


2. 3. 1.1.  AT  varied;  T  constant 

FTId  X 

For  this  condition,  the  ratio  of  the  number-of-cycl es  to  failure  Ni 
and  for  values  of  AT j  and  aT2,  respectively  is: 

-n 

=  (AT-j/A^)  (9) 

2. 3. 1.2.  aT  constant;  varied 

ilia  X 

For  two  values  of  T,  and  T,  the  ratio  of  the  cycles-to-failure 

I  max  c.  niaX 

and  N2,  respectively  is: 


N1/N2  =  exr» 


-L _ .1.  ] 

r  t 

1  max  2  max) 


(10) 


i  i  i» 
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2. 3. 1.3.  aT  varied;  T  .  constant 

min 

For  this  condition  the  ratio  of  the  cycles-to-failure  N-j  and  N£  for 
values  of  aT-|  and  aT£,  respectively,  becomes: 


N1 

( 

AT1 

-n 

exp 

Q 

r 

1  1 

n2- 

AT2  1 

R 

_ 

Jl  max  "*"2  max. 

2.3.2'.  Natural  convection 

Under  conditions  of  natural  convection  the  heat  transfer  coefficient 
is  a  function  of  the  instantaneous  temperature  difference  (aT^ )  between 
the  specimen  surface  and  the  fluid  medium.  For  values  of  the  Grasshof, 

Prandtl  and  Nusselt  numbers  appropriate  for  laboratory  studies  of  thermal 
fatigue  behavior,  the  heat  transfer  coefficient  can  be  expressed  [16] 

h  =  C1  ( AT  -j ) 1  /4  (12) 

where  C1  is  a  constant  which  depends  on  the  properties  of  the  fluid  medium 
such  as  the  viscosity,  density,  thermal  expansion,  specific  heat  and  ther¬ 
mal  conductivity  as  well  as  on  the  size  of  the  structure  being  heated  or 
cooled.  Equation  12  indicates  that  the  heat  transfer  coefficient  decreases  as 
the  difference  in  temperature  between  the  specimen  surface  and  medium  (aT^) 
decreases  during  the  transient  heat  transfer.  In  accepting  the  validity  of 
Equation  12,  it  is  implicitly  assumed  that  during  the  transient  heat  transfer, 
conditions  for  natural  free  convection  pertain,  which  essentially  is  a  steady- 
state  phenomenon. 

The  transient  thermal  stresses  in  a  solid  cylinder  under  natural  heat 
conditions  described  by  Equation,  12  were  calculated  by  finite  element 
methods  to  be  reported  in  detail  elsewhere.  It  was  found  that  the 


maximum  value  of  the  tensile  thermal  stresses  in  the  surface  of  the  cylinder 
over  the  range  0.5  <  6  <  20,  to  an  excellent  approximation  can  be  written  [17] 


-1 

a 

max 


1-v 

aEAT 


(1  .90  +  6.0/3) 


(13) 


where  the  Biot  number,  3  now  is  defined  in  terms  of  the  maximum  value  of  the 
heat  transfer  coefficient  encountered,  i.e.,  at  the  value  of  aT..  in  Equation 
12  equal  to  aT  of  the  thermal  fatigue  environment,  such  that 

6  =  RC’(AT)1/4/k  (14) 

With  3  defined  in  this  manner,  a  direct  comparison  of  Equations  7  and 
13,  shows  that  for  the  same  Biot  number,  the  magnitude  of  maximum  thermal 
stress  in  natural  convection  is  considerably  less  than  the  corresponding 
value  for  forced  convection. 

Substitution  of  Equation  13  into  Equation  6  results  in  the  number  of 
thermal  cycles-to-failure  under  conditions  of  natural  convection: 


N  = 


_ Ck 

Y2(n^2)AR2a.^n"2)/2 


1-v 

aEAT 


n  f 

'1.9 


_ 6.0k 

RC'(AT)1/4 


'  exP(Q/RTmax 


) 


(15) 


In  direct  analogy  to  the  expressions  derived  earlier  in  case  of  forced 
convection  for  the  ratio  of  the  thermal  cycles  to  failure,  corresponding 
Equations  for  condition  of  natural  convection  can  also  be  derived  as  fol¬ 
lows: 


2. 3. 2.1.  aT  varied;  T  constant 

max 

From  Equation  15,  the  ratio  of  the  thermal  cycles-to-failure  and 
Ng,  corresponding  to  values  of  aT-j  and  aT^,  becomes: 
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1.g0  +  __2^-__  1,90  +  — 

RC ' ( AT-j ) ' j  l  RC'(aT2)1/4 


Equation  16  can  be  presented  in  simplified  form  for  the  two  limiting 
values  of  the  Biot  number: 


ATi  P 


aT-j  "j  -5n/4 


e  » i 


3  <<  1 


(17a) 


(17b) 


2. 3. 2. 2.  aT  constant;  Tmax  varied 


For  two  values  of  T-j  max  and  T2  max  the  ratio  of  the  thermal  cycles- 
to-failure  N,  and  N2,  respectively,  with  the  aid  of  Equation  15  can  be 


derived  to  be: 


n;  =  exP  f  t: - 

2  1  max 


2. 3. 2. 3.  aT  varied;  T  ^  constant 


2  max 


For  convenience,  the  ratio  of  the  thermal  cycl es-to-failure  N-j  and 
N2  corresponding  to  values  of  aT^  and  aT^  and  values  of  max  and  T2max 
will  be  written  directly  for  the  high  (s  >>  1)  and  low  (b  <<  1)  Biot 
number  approximations: 


exp  f 


max  2  max 


B  >>  1 


(19a) 


11 
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|  -5/4n 


Q  _1 
R  T,  . 


«  1  (19b) 


3.  DISCUSSION 

The  analytical  results  indicate  that  the  role  of  the  individual 
parameters  specifically  n  and  Q,  which  affect  the  kinetics  of  thermal 
fatigue  crack  growth,  depends  on  the  method  by  which  a  thermal  fatigue 
curve  is  established.  For  simplicity,  relative  thermal  fatigue-life 
will  be  considered  only,  such  that  the  discussion  can  focus  in  detail 
on  the  crack  growth  variables  which  affect  the  ratios  of  cycles-to- 
failure  for  a  given  thermal  environment,  material,  specimen  and  crack 
geometry. 

As  indicated  by  Equations  9  and  17,  if  thermal  fatigue  life  is 

established  by  keeping  T  constant  and  varying  aT  by  varying  T  .  , 

ma  x  mi n 

relative  changes  in  thermal  fatigue  1 ife  are  governed  only  by  the  stress 
intensity  factor  exponent,  n.  In  contrast,  as  indicated  by  Equations 
10  and  18,  if  a  thermal  fatigue  curve  is  established  by  keeping  aT  con¬ 
stant  and  varying  Tmax  and  T^n  by  equal  amounts,  changes  in  thermal 
fatigue  life  are  influenced  only  by  the  activation  energy  for  sub- 
critical  crack  growth. 

If,  however,  as  is  common  in  practice,  thermal  fatique  behavior  is 

measured  by  keeping  T  .  constant  and  varying  T  and  aT  by  the  same 

min  m3 x 

amount,  changes  in  thermal  fatigue  life  as  indicated  by  Equations  11  and  19, 
are  affected  by  both  the  stress  intensity  exponent  as  well  as  the  activation 
energy  for  crack  growth.  Comparing  the  analytical  results  for  forced  and 
natural  convection  shows  that  for  the  latter  (with  b  <<  1)  the  role  of  the 


stress  intensity  factor  exponent  (n)  is  greater  than  in  the  former.  In  the 
high  Biot  number  approximation,  the  role  of  n  is  identical  for  both  forced 
and  natural  convection. 

The  above  conclusions  are  critical  to  the  design  of  thermal  fatigue 
experiments  and  the  analysis  of  experimental  data  obtained.  One  significant 
conclusion  is  that  thermal  fatigue  data  obtained  by  varying  aT  and  T 
simultaneously  (as  is  the  usual  case),  cannot  be  used  to  obtain  a  quantitative 
value  for  the  stress  intensity  exponent,  by  obtaining  the  value  of  the  slope 
of  log  N  vs  log  AT.  This  latter  conclusion  is  at  variance  with  the  recent 
results  of  Kamiya  and  Kamigaito  [15],  which  could  lead  to  misleading  results, 
unless  due  care  is  exercised. 

For  the  purposes  of  illustration,  a  numerical  example  will  be  considered 
on  thermal  fatigue  for  hot-pressed  polycrystall ine  silicon  nitride  subjected 
to  forced  convection  heat  transfer  by  Ammann  et  al .  [13]  by  cycling  appro¬ 
priate  specimens  at  an  initial  temperature  of  near  1600°K  into  a  fluidized 
bed  at  approximately  313°K.  Appropriate  values  for  the  stress  intensity 
exponent,  and  activation  energy  for  sub-critical  crack  growth  are  n  =  6 
and  Q  =  170  Kcal/mole,  respectively.  The  change  in  thermal  fatigue  life 
caused  by  a  simultaneous  decrease  in  T  a-nd  aT  of  20°C  will  be  considered. 

Substitution  of  the  above  values  of  n  and  Q  into  Equation  11  yields  the  re¬ 
sult  that  due  to  the  decrease  in  magnitude  of  thermal  stress  (i.e,  aT) 
thermal  fatigue  life  is  increased  by  15«  (i.e.  by  a  factor  of  1.15).  In 
contrast,  the  decrease  in  the  value  of  T  increases  thermal  fatigue-life 
by  a  factor  equal  to  1.97,  for  a  total  increase  in  fatigue-life  due  to  both 
effects  equal  to  2.26.  These  numerical  results  indicate  that  in  silicon 
nitride  the  activation  energy  for  slow  crack  growth  plays  a  far  more  im- 


portant  role  in  establishing  thermal  fatigue  behavior  than  the  stress 
intensity  exponent.  For  this  reason,  a  plot  of  log  N  vs  log  aT  of  the 
data  of  Ammann  et  al  [13]  cannot  yield  a  reliable  value  of  n  in  contrast  with 
the  findings  of  Kamiya  and  Kamigaito  [15].  Because  of  the  relatively 
small  effect  of  the  value  of  n  on  thermal  fatigue  life,  it  is  more  appro¬ 
priate  to  plot  log  N  vs  1/T  to  yield  a  value  of  Q.  Doing  this  with  the 
data  of  Ammann  et  al  [13]  results  in  a  value  of  Q  »  130  Kcal/mole  which 
is  less  than  the  value  found  by  experiment.  In  fact,  because  the  role 
of  n  is  ignored  in  this  approach,  the  activation  energy  determined  from 
the  slope  of  the  In  N  vs  1 /Tmax  P^ot  should  be  higher  than  the  value  of 
Q  determined  by  actual  measurement.  It  is  the  view  of  these  writers  that  the 
total  number  (seven)  of  data  points  in  the  set  of  data  of  Ammann  et  al  [13] 
is  too  small  to  reliably  establish  the  thermal  fatigue  behavior  of  silicon 
nitride.  Of  the  seven,  four  show  significant  data  scatter  for  nearly 
identical  values  of  T  ,  leaving  only  three  data  points  to  establish 

IDa  X 

the  effect  of  aT  or  T  on  fatigue-life.  Three  data  points  are  not  con- 
sidered  adequate  in  view  of  the  statistical  nature  of  brittle  fractjre  of 
such  materials  as  silicon  nitride  resulting  from  variations  in  crack  size, 
geometry,  orientation  and  other  variables. 

Such  statistical  effects  were  largely  eliminated  in  the  study  of 
Hasselman  et  al  [30]  of  the  thermal  fatigue  behavior  of  circular  rods  of 
soda-lime-silica  glass,  thermally  cycled  from  higher  temperature  into  a 
water  bath  at  constant  lower  temperature.  The  excessive  scatter  in  ther¬ 
mal  fatigue  data  was  reduced  by  promoting  failure  at  artificial  surface 
flaws  introduced  by  diamond  indentation.  The  data  obtained  and  the 
property  value  for  the  glass  can  be  used  to  numerically  illustrate  the 
relative  influence  of  n  and  Q  on  thermal  fatigue  life. 
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For  the  soda-1 ime-sil ica  glass  rods  investigated,  n  -  16  and  0-25 

Kcal/mole.  In  the  fatigue  experiments  T  =  433°K  and  aT  =  130°C.  The 

r  max 

relative  change  in  thermal  fatigue  life  which  results  from  a  decrease 

and  aT  of  10°C  will  be  calculated.  Assuminq  that  heat  transfer  oc- 
curred  by  forced  convection  as  the  result  of  the  motion  of  the  specimen 
through  the  water,  substitution  of  the  above  values  of  n  and  0  in  Equation 
11  yields  an  increase  in  thermal  fatigue-life  due  to  the  decrease  in  aT 
by  a  factor  of  3.6  whereas  the  decrease  in  Tmax  increases  thermal  fatigue- 
life  by  a  factor  approximately  equal  to  1.92.  In  the  case  of  soda-lime- 
silica  glass,  then,  the  stress  intensity  exponent  plays  a  more  important 
role  in  establishing  fatigue-life  than  the  activation  energy.  This  contrasts 
with  the  previous  findings  for  silicon  nitride.  Of  course,  such  an  effect 
is  expected  since  for  the  silicon  nitride  the  values  of  n  and  Q  are  lower  and 
higher  respectively  than  the  corresponding  values  for  the  soda-lime-silica- 
glass. 

The  data  for  the  soda-lime-silica  glass  rods  of  Hasselman  et  al .  [20] 
can  be  used  to  illustrate  the  erroneous  values  of  n  and  Q,  which  can  be 
obtained  unless  care  is  taken.  Figures  2a  and  2b  show  the  experimental 
data  for  the  number-of-cycles  to  failure  for  the  fifth  specimen  out  of  a 
total  set  of  nine.  Figure  2a  shows  In  N  plotted  as  a  function  of  In  aT, 
whereas  Fig.  2b  shows  the  identical  data  plotted  as  a  function  of  1/T 

Pla  X 

Except  for  the  data  for  1  cycle-to-failure  (log  N  =  0),  which  may  be 
governed  by  failure  well  before  the  end  of  the  first  cycle,  both  sets  of 
data  show  reasonable  linear  behavior.  The  slope  of  the  data  in  Fig.  2a 
yields  an  apparent  value  for  the  stress  intensity  exponent  of  22,7,  This 
is  significantly  higher  than  the  literature  value  [10]  of  n  for  this 
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material,  as  expected  since  a  plot  of  In  N  vs  In  ,\T  does  not  reflect 


changes  in  the  value  of  T  .  Figure  2b  results  in  an  apparent  value  of 

Q  =  69  Kcal/mole.  This  is  far  in  excess  of  the  experimental  value,  because 

a  plot  of  In  N  vs  1/T  does  not  reflect  the  role  of  the  stress  intensity 
r  ma  x 

exponent  in  thermal  fatigue-life.  Assuming  the  existence  of  natural  con¬ 
vection  during  the  thermal  quench  of  the  soda-lime-silica  glass  rods  results 
in  an  only  slightly  lower  value  of  n,  because  the  value  of  the  Biot  number  due 
to  the  thermal  conductivity  has  a  value  of  near  5  to  10. 

Since  the  above  numerical  examples  indicate  that  care  needs  to  be 
exercised  in  the  analysis  of  thermal  fatique  data,  the  following  recommenda¬ 
tions  are  made:  Values  of  thermal  fatigue-life  should  be  established  for 
two  basic  conditions;  1.  T  v  should  be  kept  constant  with  AT  being 

IMd  X 

varied  by  changing  T  .  ;  2.  AT  should  be  kept  constant  by  varying  T 

ini  n  n  id  x 

and  T  .  by  equal  amount.  The  data  obtained  for  the  first  condition 
should  result  in  an  unambiqious  value  of  n  whereas  the  results  for  the 
second  condition  should  aive  an  unbiased  value  for  Q.  The  data  obtained  in 
this  manner  should  be  sufficient  to  make  predictions  of  thermal  fatigue  • 
life  for  any  combination  of  Tmax,  T^  and  AT.  If  desired,  additional 
experimental  data  can  be  obtained  for  verification  of  such  predictions. 

Equation  8  and  15  indicate  that  even  for  identical  specimen  tempera¬ 
tures,  the  thermal  fatigue  life  for  natural  and  forced  convention  modes  of 
heat  transfer  will  be  different.  Therefore,  for  quantitative  evaluation  of 
thermal  fatigue-life,  the  mode  of  heat  transfer  must  be  established.  Typical¬ 
ly,  in  a  thermal  fatigue  experiment,  the  hot  specimen  is  inserted  into  a 
quenching  bath  at  lower  temperature  and  held  there  for  a  specified  period 
of  time.  There,  for  a  small  specimen  in  a  quenching  medium  (such  as  water) 
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with  a  high  heat  transfer  coefficient  where  the  maximum  stress  could  develop 
during  the  time  period  when  the  specimen  is  moving  through  the  medium, 
the  thermal  fatigue-life  will  be  controlled  by  the  forced  convection  mode  of 
heat  transfer.  On  the  other  hand  for  a  large  specimen  in  a  quenching  medium 
(such  as  oil)  with  a  low  heat  transfer  coefficient  where  the  maximum  stress 
will  develop  during  the  hold  period  after  insertion  into  the  fluid,  the 
fatigue  life  will  be  controlled  by  the  natural  convection.  This  suggests 
the  need  for  a-priori  rough  estimation  of  the  mode  of  heat  transfer  before 
the  quantitative  data  analysis  can  be  made. 

Some  materials  may  show  the  existence  of  a  pronounced  fatigue-limit, 
i.e.,  a  value  of  stress  intensity  factor  below  which  no  slow  crack  growth 
occurs.  This  appears  to  be  the  case  for  polycrystalline  mullite  investigated 
by  Kamiya  and  Kamigaito[l 5] .  In  case  of  such  a  Dronounced  fatigue-limit, 
the  expressions  presented  earlier  will  need  to  be  modified.  A  further  com¬ 
plexity  in  data  analysis  arises  for  heat  transfer  coefficient  which  may  in¬ 
dicate  a  strong  temperature  dependence.  For  many  fluid  media,  this  can  re¬ 
sult  from  the  effects  of  nucleate  boiling  and  film-formation.  The  correspond 
ing  changes  in  h  with  changes  in  T  ,  T  .  and  \T  could  well  play  a  major 
role  in  governing  fatigue-life.  In  particular,  this  could  be  the  case  with 
water,  commonly  used  as  a  quenching  medium.  Extra  caution  is  recommended 
in  analyzing  data  of  thermal  fatigue  experiments  obtained  with  water  baths. 
For  this  reason,  fluidized  beds  may  be  preferred  over  water  bath.  However, 
in  the  use  of  a  fluidized  bed,  caution  must  be  exercised  with  specimens  such 
as  glass  whose  thermal  fatigue  life  may  be  greatly  influenced  by  surface 
conditions.  The  surface  of  a  glass  specimen  may  be  damaqed  during  its  motion 
through  the  fluidized  bed  due  to  particle  impact  and  the  fatigue  life  may  be 
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greatly  reduced.  Therefore,  the  effect  of  surface  damage  due  to  Darticle 
contact  in  the  fluidized  bed  must  be  independently  evaluated  and  incorpo¬ 
rated  in  the  quantitative  evaluation  and  interpretation  of  the  experimental 
data. 

The  above  analyses  also  shed  light  on  the  effect  of  specimen  size  and 
material  properties  on  thermal  fatigue  life.  For  high  values  of  Biot  number. 
Equations  8  and  15  indicate  that  thermal  fatigue-life  is  inversely  proportional 
to  the  square  of  the  dimensions  of  the  specimen.  On  dimensional  grounds,  this 
conclusion  should  be  generally  valid  regardless  of  geometry.  For  low  values 
of  Biot  number,  freqnently  encountered  in  laboratory  studies  and  such  com¬ 
ponents  as  silicon  nitride  and  silicon  carbide  turbine  blades  and  valves,  a 
very  pronounced  size  effect  exists.  For  both  forced  and  natural  convection. 
Equations  8  and  15  suggest  that 

N  -  1/Rn+2  8  «  1  (20) 

This  result  indicates  that  major  improvements  in  thermal -fatigue  life 
can  be  achieved  by  minor  design  changes  in  the  form  of  even  slight  reduction 
in  component  size. 

In  practice,  the  materials  technologist  may  wish  to  select  the 
material  with  highest  thermal  fatigue  resistance.  For  a  specimen  of  given 
size  and  geometry  with  a  given  crack  size  and  configuration  and  given  heat 
transfer  environment,  highest  thermal  fatigue  resistance  can  be  obtained  by 
optimizing  the  relevant  material  properties  which  govern  thermal  fatigue 
failure.  Materials  with  optimum  thermal  stress  resistance  can  be  selected 
on  the  basis  of  thermal  stress  resistance  parameters  or  figures-of-merit 
which  are  available  for  a  wide  variety  of  heat  transfer  conditions  and 


modes  of  thermal  stress  failure.  Similar  figures-of-merit  can  be  defined 
for  optimizing  thermal  fatigue  resistance.  For  both  forced  and  natural 
convection  these  can  be  obtained  from  Equations  8  and  15. 

'  n 

■  exp(Q/RTmax)  e  «  1  (21a) 

n 

exp(Q/RTmax)  8  >:>  1  (21b) 

These  figures-of-merit  indicate  that  thermal  fatigue-life  is  governed 
by  as  many  as  seven  material  properties.  Without  quantitative  information 
on  these  properties,  estimates  of  thermal  fa Jgue-life  are  not  feasible. 

Of  interest  to  note  is  that  high  thermal  fatigue  resistance  requires  high 
values  of  the  thermal  diffusivity  «,  thermal  conductivity  k,  stress  intensity 
factor  exponent  n  and  activation  energy  0  with  low  values  of  constant  A, 
Young's  modulus  E  and  coefficient  of  thermal  exnansion  a. 

In  this  respect,  for  basic  studies  of  thermal  fatigue  life,  materials 
with  low  thermal  diffusivity  and  conductivity  may  be  most  useful  since  thermal 
fatigue  curves  can  be  established  with  minimum  number  of  cycles.  Glassy 
materials,  for  this  reason,  appear  most  useful. 
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Fig.  2b.  Thermal  fatigue  life  of  soda-lime-silica  glass  specimen 
subjected  to  repeated  water  quench  (after  ref.  [20]). 
a;  as  a  function  of  the  initial  temperature  difference, 
b;  as  a  function  of  the  maximum  temperature. 
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ABSTRACT 


An  analysis  is  presented  for  the  temperatures  and  thermal  stresses 
in  a  partially  absorbing  flat  plate  subjected  to  normally  incident  cyclic 
thermal  radiation  on  the  front  face,  and  cooled  by  convection  on  the  rear 
face. 

The  resulting  temperature  and  stress  responses  are  cyclic  in  nature 
with  the  temperatures  and  stresses  in  the  center  and  the  back  face  lagging 
behind  those  in  the  front  face.  The  amplitude  of  the  temperature 
fluctuation  is  found  to  be  maximum  in  the  front  face.  During  the  initial 
thermal  cycles  the  maximum  stress  in  the  plate  is  compressive.  In  con¬ 
trast,  when  the  plate  approaches  thermal  equilibrium  after  many  cycles, 
the  maximum  stress  is  tensile.  The  values  of  the  maximum  tensile  and 
compressive  stresses  were  found  to  depend  on  the  frequency  of  the  incident 
cyclic  radiation. 
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INTRODUCTION 


Radiation  heat  transfer  can  result  in  transient  and  steady-state  thermal 
stresses.  This  problem  is  particularly  severe  in  components  and  structures 
such  as  lenses,  windows,  heat  exchangers,  etc.,  subjected  to  intense  radiation 
from  lasers,  solar  energy  concentrators  and  other  sources. 

Solutions  for  the  thermal  stresses  due  to  suddenly  applied  constant 
flux  black-body  radiation  were  presented  earlier  for  materials  totally 
opaque  to  the  incoming  radiation  [1] .  This  solution  was  extended  to  in¬ 
clude  materials  opaque  above  a  given  wavelength  and  totally  transparent 
below  this  wavelength  of  black-body  radiation  [2],  More  recently,  for 
a  variety  of  boundary  conditions,  solutions  were  obtained  for  the  thermal 
stresses  in  materials  suddenly  subjected  to  a  constant  radiative  heat  flux, 
absorbed  internally  as  described  by  the  absorption  coefficient  [3-6]. 

In  practice,  radiative  heat  fluxes  can  be  constant  or  vary  with  time. 

This  latter  condition  is  encountered  for  pulsed  lasers  or  solar-energy  con¬ 
centrators  under  conditions  of  intermittent  cloud  cover.  Such  rapid 
changes  in  radiation  intensity  are  expected  to  give  rise  to  transient 
thermal  stresses.  The  purpose  of  this  paper  is  to  present  an  analysis 
of  the  thermal  stresses  in  a  partially  transparent  flat  plate  subjected  to 
thermal  radiation  with  an  intensity  which  varies  sinusoidally  in  time. 

ANALYSIS 

Boundary  and  Initial  Conditions 

The  flat  plate,  infinite  in  extent,  is  located  in  the  yz  plane,  and 
-a  <  x  <  a.  At  the  front  face  (x  *  -a)  the  plate  is  subjected  at  time 


y. 


..  .v  .*■ 
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t  =  0  to  normally  incident,  spatially  uniform  radiation  which  oscillates 
in  time  with  angular  frequency  w. 

At  the  back  face  (x  =  a) ,  the  plate  is  cooled  by  Newtonian  convt  ction 
with  a  heat  transfer  coefficient  h.  As  for  the  previous  analysis  [3-6] 
the  optical  properties  of  the  plate  are  assumed  to  be  "grey,"  i.e.,  independent 
of  wavelength.  The  reflectivity  of  the  plate  is  assumed  sufficiently  low 
so  that  the  effect  of  multiple  reflections  within  the  plate  can  be  neglected. 

The  material  properties  that  affect  the  temperature  distribution  and  thermal 
stresses,  such  as  the  emissivity  (absorptivity),  absorption  coefficient, 
thermal  conductivity  and  thermal  diffusivity,  coefficient  of  thermal  expansion. 
Young’s  modulus  and  Poisson's  ratio  are  assumed  to  be  independent  of  tempera¬ 
ture.  It  is  also  assumed  that  the  temperatures  in  the  plate  remain  sufficiently 
low  so  that  the  intensity  of  any  re-emitted  radiation,  to  a  first  approximation, 
is  negligible  in  comparison  to  the  intensity  of  the  incident  radiation. 

The  intensity  q(t)  of  the  radiation  as  it  approaches  the  plate  is  assumed 
to  be  of  the  form 

q(t)  =  qQ[l  +  «  sin(uit)  ]  (1) 

where  0  *  6  S  1  and  t  is  the  time.  For  this  study,  6  will  be  taken  as  1,  so 
that  the  magnitude  of  incident  radiation  heat  flux  varies  between  the  values 
2q^  and  0. 

The  intensity  of  the  heat  flux  q  at  x  within  the  plate  due  to  the  thermal 
radiation  entering  at  x  ■  -a  is 

q(x,t)  *  eqoe  ^a+x) [i  +  $  sin(mt)]  (2) 

where  e  ■  1-r  with  r  being  the  reflectivity  and  y  is  the  absorption  coefficient. 
The  rate  of  internal  heat  absorption  per  unit  volume  g"'  within  the  plate  at 


x  is  thus 
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g"'(x,t)  =  yeqQe  [1  +  6  sin(wt)]  (3) 

Solutions  for  the  transient  temperature  T  (x,t)  are  obtained  by  solving 
the  heat  conduction  equation  [7] 


2 

rj 

3x" 


9tl  ,  g"'(x,t)  .  1  3T 
2  k  k  3x 


(4) 


where  k  is  the  thermal  conductivity  and  k  is  the  thermal  diffusivity. 

For  the  given  heating  and  cooling  conditions  the  initial  and  boundary 
conditions  are 


T(x, 0)  =  T 


3T 

Dx 

3T 


(-a, t)  =  0 


(5) 

(6) 

(7) 


ax  ‘  -  k'T(a’t>-T°! 

where  h  is  the  heat  transfer  coefficient. 

Solutions  for  the  transient  temperatures  were  obtained  using  standard 
techniques  [7] 


T(x, t)  =  Tq  +  E  Bn(t)  cos(Xnx) 
n=l 

where  X  are  the  roots  of  the  transcendental  equation 
n 

X  tan(2X  a)  =  -r 
n  n  k 
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an<^  N  =  — —  [Sin  (2X  a)  cos  (2X  a)  +  2X  a] 
n  2X  n  n  n 

n 


(12) 


THERMAL  STRESSES 

Expressions  for  the  thermal  stresses  o  (x,t)  are  obtained  [8] 

y  > 2 

by  substitution  of  eq.(8)for  T  into 

a  a 

o  =  [-T  +  •—  j  Tdx  +  /  Txdx] 

y’Z  ^  23  -a  2a3  -a 

which  yields 


(13) 


o  =  -  £  ..B  (t)cos[X  (x+a)  ] 

y , z  1-v  n=l  n  1  n  J 
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™  For  convenience,  the  numerical  results  are  reported  in  terms  of  the 

non-dimensional  temperature 

*  (T-T  )k 

T  =  - — 

cqoa 

the  non-dimensional  stress 


*  _  g(l-v)k 

aEqQca 

*  2 

the  non-dimensional  time,  t  =  Kt/a  and  the  non-dimensional  distance, 
x/a . 


RESULTS  AND  DISCUSSION 


The  analytical  results  for  the  temperatures  and  stresses  obtained 
above  are  most  easily  illustrated  and  discussed  by  means  of  a  numerical 
example.  In  order  to  keep  the  amount  of  data  to  be  presented  within 
reasonable  bounds,  emphasis  will  be  placed  on  the  effect  of  the  cyclic 
nature  of  the  incident  radiation  for  single  values  of  the  other  parameters 
(k,h,K,a, pa) .  These  latter  quantities  were  taken  as  k  =  0.3  watts. cm  ^°C  \ 
h  =  6  x  10  watts. cm  °C-^,  k  =  0.1  cm.  S  and  a=  1  cm.  The  value  of 
optical  thickness  pa  was  chosen  as  3,  which  is  near  the  value  which  re¬ 
sulted  in  the  maximum  magnitude  of  stresses  under  steady-state  conditions 
for  a  flat  plate  asymmetrically  heated  by  a  constant  radiative  heat 
flux  and  cooled  by  convection  on  either  the  back  or  front  face  [5,6]. 

Figure  1  for  an  arbitrarily  selected  value  of  w  =  0.21  s  1  illustrates 

* 

the  magnitude  of  the  temperature  as  a  function  of  t  for  three  different 

time  intervals  and  three  positions  within  the  plate,  x  =  -a,  0,  and  a. 

Figure  la  shows  the  temperature  response  during  the  first  few  cycles. 

Figure  lc  corresponds  to  the  magnitude  of  temperature  at  large  values 

of  time  for  which  the  temperatures  in  the  plate  approach  the  values 

encountered  at  t  =  °°.  This  condition  corresponds  to  a  time-averaged 

* 

back-face  temperature  T  =  49.9  at  which  the  average  convective  heat 
flux  equals  the  time-averaged  radiative  heat  flux  absorbed  within  the 
plate.  For  pa  =  3  for  this  numerical  example,  this  corresponds  to  99.8% 
of  the  radiant  flux  which  enters  the  plate  at  x  =  -a.  As  shown  in  Fig.  1 
for  all  three  time  intervals,  as  expected,  the  cyclic  temperature  response 
at  the  center  and  back-face  lags  behind  the  corresponding  value  for  the 
front  face.  Also  the  rate  of  rise  of  average  temperature  decreases  with 
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increasing  time  as  the  rate  of  heat  removal  by  convection  increases  with 


increasing  back-face  temperature. 

Figure  2  shows  the  spatial  distribution  and  amplitude  of  the  temperature 
oscillation  for  the  time- interval  in  Fig.  lc.  The  amplitude  of  the  temperature 
oscillation  is  greatest  at  the  front-face  and  decreases  with  increasing  dis¬ 
tance  towards  the  back-face.  The  temperature  distribution  exhibits  a  pro¬ 
nounced  downward  concavity.  This  results  from  the  asymmetric  distribution 
of  the  rate  of  internal  heat  generation. 

Figure  3  shows  the  thermal  stresses  at  x  =  -a,  0  and  a  for  the  time 
intervals  shown  in  Fig.  1.  In  Fig.  3c  the  cycle- to-cycle  variation  in 
the  magnitude  of  thermal  stress  is  less  than  0.02%.  Of  interest  is  the 
general  reversal  of  stress  as  time  proceeds.  The  stresses  in  the  center 
of  the  plate,  which  are  primarily  tensile  during  the  first  few  cycles  of 
the  radiation  become  primarily  compressive  at  the  higher  values  of  time  as  the 
plate  approaches  its  final  maximum  temperature.  The  opposite  effect  occurs 
for  the  stresses  at  the  back  face.  At  the  front  of  the  plate  the  stresses 
initially  cycle  between  tension  and  compression,  to  become  entirely  tensile 
at  large  values  of  time.  Such  tensile  stresses  occur  in  the  front  face, 
which  is  at  the  highest  temperature,  because  of  the  downward  concavity  of 
the  temperature  distribution  [9].  The  stress  values  shown  in  Fig.  3  were 

_2 

calculated  for  a  given  value  of  heat  transfer  coefficient  (h=  0.006  watts. cm. 

°C  1).  However,  as  indicated  by  the  results  of  earlier  studies  of  time- invariant 
radiation  [5,6],  the  magnitude  of  the  stresses  are  not  expected  to  depend 
strongly  on  the  heat  transfer  coefficient  values  encountered  in  practice. 

During  the  first  few  cycles,  the  back  face  temperature  is  not  sufficiently 
high  that  the  temperature  distribution  is  affected  strongly  by  the  convective 
heat  losses.  Similarly,  at  the  large  values  of  time,  although  the  absolute 
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values  of  temperature  encountered  are  an  inverse  function  of  the  heat 
transfer  coefficient,  the  temperature  profile,  which  controls  the  magnitude 
of  thermal  stresses,  is  governed  by  the  distribution  of  internal  heat 
generation  rather  than  the  heat  transfer  coefficient. 

Figure  4  shows  the  maximum  values  of  the  tensile  and  compressive 
stresses  encountered  during  each  cycle  with  increasing  total  number  of 
cycles  of  incident  radiation.  During  approximately  the  first  25  cycles, 
plate  failure  most  likely  will  occur  in  compression  unless  the  plate  is 
composed  of  a  highly  brittle  material  with  high  ratio  of  compressive  to 
tensile  strength.  After  more  than  about  twenty-five  cycles,  however,  failure 
is  expected  in  a  tensile  mode  regardless  of  whether  the  plate  is  composed 
of  brittle  or  ductile  material. 

! 

The  magnitude  of  the  thermal  stresses  are  expected  to  be  a  function  of  the 
frequency  of  incident  radiation  in  view  of  the  expected  coupling  between  the 
time  period  of  a  thermal  cycle  and  the  transit  time  of  a  thermal  wave  through 
the  plate.  This  latter  quantity  is  a  function  of  the  plate  dimensions  and  the 
thermal  diffusivity.  Evidence  for  the  validity  of  this  hypothesis  is  shown  in 
Figures  5a  and  5b.  Figure  5a  shows  the  dependence  of  the  maximum  tensile  and 
compressive  stresses  during  the  first  cycle  in  the  center  and  the  back  face  of 
the  plate,  respectively  as  a  function  of  angular  frequency,  w.  Figure  5b  shows 
the  corresponding  stresses  at  steady  state  conditions  in  the  front  face  and 
the  center  of  the  plate  respectively.  It  can  be  noted  that  for  this  numerical 
example  at  steady  state,  a  maximum  in  tensile  stress  value  is  obtained  at  an 
angular  frequency  of  w  «  0.6  (a  period  of  approx.  10  sec.)  whereas  a  peak  in 

i 

compressive  stress  value  occurs  at  an  angular  frequency  of  a  =  0.1.  At  u  -  0 
and  <*>,  the  plate  in  effect  is  subjected  to  a  constant  heat  flux,  qQ,  with  a 
corresponding  value  of  a*  =  0.167  in  agreement  with  earlier  results  [5], 
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The  results  of  Fig.  5  clearly  indicate  that  for  the  same  time  averaged 
radiative  heat  flux,  the  stresses  for  cyclic  conditions  can  exceed  those 
for  a  constant  heat  flux  by  a  considerable  margin. 
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Fig.  2.  Envelope  of  the  temperature 
oscillation  in  a  partially 
absorbing  flat  plate  heated 
in  front  by  normally  incident 
cyclic  radiation  and  cooled  at 
the  rear  surface  by  convection 
corresponding  to  the  time  interval 
in  figure  1c:  h  =  0.006  W.cm.~^°C 
k  =  0.3  W.cm.-loC_1,  k  =  0.1  cm.  s 
id  =  0.21  s.~l.  pa  =  3  and  a  =  1  cm 
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Fig.  3.  Transient  thermal  stresses  in  a 
partially  absorbing  flat  plate 
heated  in  front  by  normally  in¬ 
cident  cyclic  thermal  radiation 
and  cooled  at  the  rear  surface 
by  convection  for  three  different 
time  intervals:  h  =  0.006  W.cm.~2 
°C._1.  k  =  0.3  W.cm.-loc.-1,  k  = 
0.1  cm.^s.-!,  a:  =  0.21  s.“l,  ya  *= 

3  and  a  *  1  cm. 
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y.  4.  Maximum  tensile  and  compressive 
stresses  in  a  partially  absorb¬ 
ing  flat  plate  heated  in  front 
by  normally  incident  cyclic 
thermal  radiation  and  cooled  at 
the  rear  surface  by  convection: 
h  =  0.006  W.cm.-2°C.-1,  k  =  0.3 
W.cm._l°C.~l,  k  -  0.1  cm.2s. 
w  =  0.21s.--*,  pa  =  3  and  a  =  1 
cm. 


ANGULAR  FREQUENCY  Irod 


Variation  of  the  maximum  tensile  and 
compressive  stresses  with  the  fre¬ 
quency  of  the  incident  cyclic  thermal 
radiation  in  a  partially  absorbing 
flat  plate  heated  in  front  and  cooled 
at  the  rear  surface  by  convection: 
h  =  0.006  W.cm.-2°C.-1.  k  =  0.3  W.cm. 
°C-1,  k  =  0.1  cm.^s.“l,  ya  =  3  and  a 
1  cm.  (a)  for  first  cycle,  (b)  for 
steady  state. 
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ANISOTROPY  OF  THE  THERMAL  DIFFUSIVITY 
OF  Si-jN^-BN  COMPOSITES 

K.  Niihara,  L.  D.  Bentsen,  D.  P.  H.  Hasselraan 
Department  of  Materials  Engineering 
Virginia  Polytechnic  Institute  and  State  University 

K.  Mazdiyasni 

Air  Force-Wright  Aeronautical  Laboratory 
Wright-Pat terson  Air  Force  Base,  Ohio  45433 

The  thermal  diffusivity  of  silicon  nitride-boron  nitride  composites 
is  shown  to  exhibit  considerable  anisotropy  with  respect  to  the  hot- 
pressing  direction  due  to  the  preferred  orientation  of  the  boron  nitride 
inclusions. 

During  the  powder -metallurgical  processing  of  composites  involving 
mechanical  compaction,  the  particles  of  one  or  both  components  can  take 
on  a  preferred  crystallographic  orientation  with  respect  to  the  direction 
of  compaction.  Materials  such  as  mica  [l],  graphite  [2],  or  boron  nitride, 
frequently  exhibit  a  pronounced  anisotropy  of  many  of  their  properties 
such  as  elastic  moduli  [3],  thermal  expansion  [4],  thermal  conductivity 
and  diffusivity  [5,6,7].  For  this  reason,  it  is  expected  that  the  bulk 
properties  of  composites,  in  which  such  particles  take  on  a  preferred 
orientation,  also  will  be  anisotropic.  Evidence  for  this  effect  in  the 
form  of  experimental  data  for  the  thermal  diffusivity  of  the  silicon 
nitride  containing  a  dispersed  phase  of  boron  nitride  is  presented. 
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The  composites  were  made  by  hot-pressing  mixtures  of  silicon  nitride 
powder  of  99.9%  purity  and  average  particle  size  of  0.5  pm  and  boron 
nitride  of  1  pm  particle  size.  The  silicon  nitride  powder  contained 
6  wt.%  of  an  alkoxy-derived  Ce02  as  hot-pressing  aid.  The  powder  mix¬ 
tures  were  dry-milled  in  polyethylene  jars  with  alumina  balls  for  a 
period  of  4  hours.  Discs  5.1  to  7.6  cm  in  diameter  by  0.6  cm  thick 

were  uniaxially  hot-pressed  in  graphite  dies  coated  with  BN  and  lined 
§ 

with  graphite  foil  in  a  nitrogen  atmosphere.  The  hot-pressing  schedule 
consisted  of  slow  heating  to  1200°C  at  a  pressure  of  about  3  MPa  fol¬ 
lowed  by  rapid  heating  to  1700-1750°C  and  holding  for  a  period  of  30  to 
60  minutes  at  about  30  MPa.  The  hot-pressed  discs  were  allowed  to 
cool  slowly  under  pressure  within  the  hot-press. 

The  preferred  orientations  of  the  boron  nitride  and  silicon  nitride 
were  established  by  x-ray  analysis  of  surfaces  cut  parallel  and  perpendi¬ 
cular  to  the  hot-pressing  direction.  For  all  five  compositions  studied. 
Figs,  la  and  lb  show  the  ratios  of  the  relative  summed  intensities 
normal  and  parallel  to  the  hot-pressing  direction  for  the  boron  nitride 
and  silicon  nitride,  respectively.  The  boron  nitride  inclusions  are 
highly  oriented  with  the  c-axis  parallel  to  the  direction  of  hot-press¬ 
ing.  The  silicon  nitride  also  has  a  preferred  orientation  which  has 
been  observed  by  other  investigators  [8,9].  However,  the  degree  of 
preferred  orientation  is  reduced  by  the  presence  of  the  boron  nitride 
particles,  probably  because  they  interfere  with  the  direction  of  growth 

*SN  402,  GTE  Laboratories,  Towanda,  PA. 

**Cerac  Inc.,  Menonomee  Falls,  WI. 

*Grafoil,  Union  Carbide  Corp. ,  Chicago,  IL. 
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of  the  g-Si^N^  grains  perpendicular  to  the  hot-pressing  direction.  Trans¬ 
mission  electron  microscopy  combined  with  electron-beam  diffraction  of  ion- 
thinned  sections  confirmed  the  preferred  orientation  of  the  BN  particles. 

Also,  the  aspect  ratio  of  the  BN  particles  was  found  to  have  a  value  near 
three.  It  is  postulated  that  this  particle  shape  contributed  to  the  preferred 
orientation  achieved  due  to  particle  rearrangement  during  hot-pressing. 

Specimens  for  the  measurement  of  the  thermal  diffusivity  in  the  form 
of  disks  approximately  10-12  mm  in  diameter  by  2  mm  thick  were  cut  from 
each  of  the  hot-pressed  blanks.  For  all  compositions,  three  specimens  were 
cut  for  heat  flow  parallel  and  perpendicular  to  the  hot-pressing  direction. 

The  thermal  diffusivity  at  room  temperature  was  measured  by  the  laser-flash 
technique  [10]  using  equipment  described  elsewhere  [ll,12]. 

Figure  2  shows  the  experimental  data,  which  represent  the  average  of 
the  three  specimens  for  each  composition  and  orientation.  Specimen-to- 
specimen  variation  was  less  than  5%.  In  agreement  with  previous  findings  [12], 
the  silicon  nitride  without  BN  exhibits  anisotropic  thermal  diffusivity.  This  is 
due  to  the  preferred  orientation  of  the  8-grains,  confirmed  by  the  x-ray  data 
given  in  Fig.  1.  Upon  addition  of  BN  the  anisotropy  of  the  thermal  diffusivity 
increases  to  be  most  pronounced  at  30  wt.%  BN.  At  low  BN  content  the  thermal 
diffusivity  perpendicular  as  well  as  parallel  to  the  hot-pressing  direction 
shows  a  decrease  with  increasing  BN  content.  This  is  thought  to  be  due  to  the 
combination  of  two  separate  effects.  The  first  effect  results  from  the  reduc¬ 
tion  in  the  degree  of  preferred  orientation  of  the  silicon  nitride.  The  second 
effect  results  from  the  absolute  values  and  anisotropy  of  the  thermal  diffusi¬ 
vity  of  the  BN  particles.  As  judged  by  data  [6]  for  pyrolytic  BN,  the  thermal 
diffusivity  in  the  basal  plane  exceeds  the  value  for  silicon  nitride  by  about  a 
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factor  of  5.  On  the  other  hand,  the  thermal  diffusivity  of  BN  perpendicular 
to  the  basal  plane  is  only  about  one  twentieth  of  that  of  silicon  nitride. 

For  the  present  samples  then  it  would  be  expected  that  parallel  to  the  hot- 
pressing  direc tion  the  BN  particles  would  cause  a  decrease  in  thermal  dif¬ 
fusivity  of  the  composite,  as  observed. 

Perpendicular  to  the  hot-pressing  direction,  however,  an  increase 
in  the  thermal  diffusivity  would  be  expected  as  observed  at  the  higher 
value  of  BN  content.  The  observed  decrease  at  the  lower  volume  fraction 
of  BN,  is  due  to  the  combined  effects  of  the  decrease  in  degree  of  pre¬ 
ferred  orientation  of  both  the  silicon  nitride  and  boron  nitride. 

Attempts  to  calculate  the  thermal  diffusivity  from  the  thermal  con¬ 
ductivity,  specific  heat  and  density  for  silicon  nitride  and  pyrolytic  BN, 
by  using  composite  theory  did  not  give  good  quantitative  agreement  with 
the  observed  data.  This  suggests  that  a  third  effect  may  also  be  operative 
such  as  an  interfacial  effect  resulting  from  less  than  perfect  adhesion 
between  the  silicon  nitride  and  boron  nitride.  Most  likely  the  lack  of  agree¬ 
ment  between  theory  and  experiment  occurs  because  the  BN  particles  in  the 
present  composites  and  pyrolytic  BN  may  have  different  degrees  of  structural 
perfection. 

The  anisotropy  in  thermal  diffusivity  observed  for  the  present  composites 
possibly  could  be  used  to  advantage  for  structural  components  which  should 
show  good  thermal  insulating  behavior  in  one  direction  and  high  thermal 
conductivity  in  a  perpendicular  direction,  for  the  purpose  of  improving 
thermal  stress  resistance  depending  on  the  direction  of  heat  flow. 
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Fig. 


2.  Room  temperature  thermal  diffusivity 
of  SijN^-BN  composites  for  two  orienta 
tion  relative  to  the  hot-pressing 
direction . 
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ABSTRACT 

The  room  temperature  thermal  d i . . i v i l v  o!  silicon  carbide- 
silicon  composites  was  measured  bv  the  laser-! lush  method.  The 
experimental  results  indicate  that  the  thermal  diliusivitv  and  con¬ 
ductivity  of  these  composites  decreases  strongly  with  increasing 
level  of  impurity  content  in  both  the  silicon  and  the  silicon 
carbide.  This  suggests  that  the  heat  transfer  propertio  s  ol  these 
composites  can  be  tailored  by  control  over  impurity  levels. 


INTRODUCTION 

Increased  efficiency  of  energy-conversion  systems  can  be 
achieved  by  increasing  their  operating  temperature  levels.  For 
this  reason  a  demand  has  been  created  for  engineering  materials 
with  mechanical  properties  and  corrosion  and  oxidation  resistance 
well  in  excess  of  those  exhibited  by  super-alloys.  Materials  such 
as  the  refractory  nitrides  and  carbides,  in  view  of  their  chemical 
and  structural  stability  at  high  temperature,  appear  to  be  excel¬ 
lent  candidate  materials  for  such  applications  as  the  all-ceramic 
turbine  engine  or  internal  combustion  engine  with  retro-fitted 
ceramic  components. 

The  in-service  performance  of  these  materials  depends  not  only 
on  their  chemical  and  mechanical  properties,  but  also  on  their 
thermal  properties  as  well.  These  latter  include  the  coefficient 
of  thermal  expansion  and  the  thermal  conductivity  and  diffusivity. 
Low  values  of  the  thermal  conductivity  lead  to  improved  energy- 
conversion  efficiencies  since  heat  losses  are  kept  to  a  minimum. 


Due  to  their  brittle  nature,  however,  the  above  materials  are  high¬ 
ly  susceptible  to  catastrophic  failure  due  to  thermal  stresses  of 
high  magnitude  which  result  from  non-linear  transient  or  steady- 
state  temperature  distributions.  In  order  to  minimize  the  possi¬ 
bility  of  this  mode  of  failure,  the  optimum  candidate  materials 
should  have  values  of  the  coefficient  of  thermal  expansion  as  low 
as  possible  in  combination  with  values  of  the  thermal  conductivity 
and  thermal  dlffusivity  as  high  as  possible.  This  latter  require¬ 
ment  is  incompatible  with  the  requirement  of  low  thermal  conducti¬ 
vity  for  high  operating  efficiency.  Clearly,  appropriate  trade¬ 
offs  must  be  made.  The  ability  to  make  such  trade-offs,  in  turn, 
requires  quantitative  information  for  the  values  of  the  thermal 
conductivity  and  thermal  diffusivity  and  a  detailed  understanding 
of  the  materials  and  intrinsic  variables  which  affect  their  values. 
In  particular,  such  effects  at  higher  levels  of  temperature  are  of 
considerable  technical  significance. 

Heat  transport  through  a  material  occurs  primarily  by  phonon, 
photon  or  electron  transport.  The  heat  flux  which  results  from 
these  mechanisms  depends  strongly  on  the  associated  specific  heat, 
the  temperature,  and  the  existence  of  structural  and  chemical  im¬ 
perfections  such  as  vacancies,  dislocations,  grain  boundaries, 
foreign  atoms,  optical  discontinuities  and  any  other  factors  which 
contribute  to  phonon,  photon  and  electron  scattering  [1,2].  At  the 
microstructural  level,  the  conduction  of  heat  is  affected  by  the 
presence  of  pores,  Inclusions  and  cracks,  grain  boundary  phases  and 
preferred  crystallographic  orientation  [2,3-9]. 

Many  of  the  variables  which  affect  the  conduction  of  heat 
depend  on  processing  history,  which  suggests  that  a  measure  of 
control  can  be  exerted  over  the  value  of  the  thermal  conductivity 
or  thermal  diffusivity  of  a  material  required  for  a  given  applica¬ 
tion.  The  purpose  of  this  paper  is  to  present  experimental  data 
for  some  of  the  variables  that  affect  the  heat  conduction  properties 
of  silicon  carbide-silicon  composites  which  have  demonstrated  po¬ 
tential  for  ceramic  turbine  applications. 


MATERIALS  AND  EXPERIMENTAL  PROCEDURES 

The  SiC-Si  composites  were  made  by  a  warm  molding  process 
using  an  industrial  grade  silicon  and  silicon  carbide  powder  made 
by  the  Acheson  process.  For  two  different  series  of  samples,  400 
and  1000  mesh  size  silicon  carbide  powders  were  used.  A  single 
sample  was  also  prepared  from  1200  mesh  SiC  powder.  The  silicon 
content  of  the  composites  was  calculated  from  the  measured  density 
of  the  composites  and  values  of  the  densities  of  silicon  and  sili¬ 
con  carbide  of  2.33  and  3.21  g/cm^,  respectively.  For  the  compo¬ 
sites  made  with  the  silicon  carbide  particle  size  of  400,  1000  and 
1200  mesh,  the  final  size  of  the  silicon  inclusions  within  the  SiC 
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Fig.  1.  Optical  micrographs  of  silicon  carbide-silicon  composites 
with  average  Si  particle  size  of  a:  15  pm;  b:  5  pm  and 
c:  A  pm. 

matrix  was  approximately  15,  5  and  A  pm,  respectively.  Typical 
photo-micrographs  are  shown  in  Fig.  1. 

For  purposes  of  interpretation  of  the  data  additional  samples 
included  in  the  study  consisted  of  dense  sintered  a-SiC,  and  a 
representative  sample  of  the  silicon  metal  used  for  the  prepara¬ 
tion  of  the  SiC-Si  composites.  The  principal  impurities  in  this 
silicon  metal  as  determined  by  spec t ro-cliemical  means  are  listed 
in  Table  1.  An  additional  sample  consisted  of  zone-refined  sili¬ 
con  of  at  least  99.99%  purity. 

TABLE  1.  Principal  Impurities  in  Industrial  Grade  Silicon  Used 
For  Preparation  of  Silicon  Carbide-Silicon  Composites 


ELEMENT 

PERCENT 

ELEMENT 

PERCENT 

Boron 

0.6 

Ti tanium 

0.00A 

Magnesium 

0.001 

Vanadium 

0.001 

Manganese 

0.01 

Copper 

0.1 

Iron 

0.1 

Chromi urn 

0.002 

Nickel 

0.00A 

Calcium 

0.06 

Aluminum 

0.01 

The  heat  conduction  properties  of  the  SiC-Si  composites  were 
determined  by  measurements  of  the  thermal  diffuslvlty  by  the  laser- 
flash  method  [10],  using  a  glass-Nd  laser.  The  specimens  of 
appropriate  size  and  geometry  were  carbon-coated  to  prevent  direct 


r 


transmission  of  the  laser-beam.  For  measurements  of  the  thermal 
diffusivity  above  room  temperature,  the  specimens  were  held  in  a 
graphite  resistance  furnace  containing  a  nitrogen  atmosphere.  At 
room  temperature  and  up  to  about  500°C  the  transient  temperature 
of  the  specimen  was  monitored  with  a  liquid  N2  cooled  InSb  infra¬ 
red  detector.  Above  500°C  the  specimen  temperature  was  monitored 
with  a  Si-photodiode. 


EXPERIMENTAL  RESULTS  AND  DISCUSSION 

Figures  2  and  3  show  the  dependence  of  the  thermal  diffusivity 
at  room  temperature  on  silicon  content  for  the  two  series  of  sili¬ 
con  carbide-silicon  composites  with  Si  inclusion  sizes  of  15  and 
5  pm,  respectively.  Table  2  lists  the  data  for  the  thermal  dif¬ 
fusivity  at  room  temperature  for  the  single  specimens  of  single¬ 
phase  a-SiC,  the  silicon  used  for  the  preparation  of  the  specimens, 
the  zone-refined  silicon  and  the  SiC-Si  specimen  with  the  4  ym  Si 
inclusions, which  contained  20.82  vol.%  Si.  For  the  latter  mate¬ 
rial,  the  thermal  diffusivity  was  found  to  be  independent  of 
orientation. 
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Fig.  2.  Thermal  diffusivity  at  room  temperature  of  silicon  carbide 
with  15  ym  Si  inclusions  as  a  function  of  silicon  content. 
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Fig.  3.  Thermal  diffusivity  at  room  temperature  of  silicon  carbide 
with  5  pm  Si  inclusions  as  a  function  of  silicon  content. 


TABLE  2.  Values  for  the  Thermal  Diffusivity  at  Room  Temperature 
for  Individual  Samples  of  Silicon,  SiC  and  SiC-Si  Composites 


MATERIALS 


2  -1 

THERMAL  DIFFUSIVITY  (cm  s  ) 


Single  Phase  Silicon  Carbide 


Industrial  Silicon 


0.694 


0.181 


Zone-refined  Silicon 


0.668 


SiC-Si  composites  (8  specimens) 


0.774  ±  8% 


Coefficient  of  variation. 

Upon  comparison  with  the  value  for  the  thermal  diffusivity  of 
the  single-phase  a-SiC  and  the  inverse  dependence  of  the  thermal 
diffusivity  on  silicon  content,  the  data  shown  in  Figs.  2  and  3 
indicate  that  the  addition  of  the  silicon  dispersed  phase  to  the 


silicon  carbide  matrix  results  in  a  significant  decrease  in  the 
thermal  diffusivity.  This  is  expected  from  the  general  theory  of 
the  thermal  conductivity  of  composites,  which  shows  that  a  dispersed 
phase  with  low  thermal  conductivity  within  a  matrix  of  higher  ther¬ 
mal  conductivity  will  lead  to  a  decrease  in  overall  thermal  diffusi¬ 
vity  of  the  composite,  regardless  of  the  nature  of  the  distribution 
of  the  dispersed  phase  [3-6].  The  same  conclusion  applies  to  the 
thermal  diffusivity  of  composites  for  components  with  comparable 
values  of  the  specific  heat  per  unit  volume.  This  latter  condi¬ 
tion  holds  for  the  present  samples  witli  values  for  the  thermal 
diffusivity  of  the  silicon  and  silicon  carbide  of  0.181  and  0.694 
cm^s-!,  respectively.  For  this  reason,  the  effect  of  silicon  con¬ 
tent  on  the  thermal  diffusivity  shown  in  Figs.  2  and  3  is  in  general 
agreement  with  composite  theory. 

Comparison  of  the  data  shown  in  Figs.  2  and  3  shows  that  for 
a  given  silicon  content,  the  thermal  diffusivity  is  a  function  of 
the  silicon  particle  size.  Such  an  effect  of  inclusion  size  is  not 
predicted  by  the  theory  of  the  thermal  conductivity  of  composites, 
unless  such  a  change  in  inclusion  size  is  accompanied  also  'by  a 
corresponding  change  in  the  geometry  or  orientation  of  the  inclu¬ 
sions.  The  photo-micrographs  in  Fig.  1  show  no  conclusive  evidence 
for  such  a  change  in  geometry  or  orientation.  Also,  such  an  effect 
of  orientation  should  result  in  an  anisotropy  of  the  thermal  dif¬ 
fusivity,  which  was  not  observed.  If,  however,  the  hypothesis  that 
a  change  in  Si  inclusion  size  is  accompanied  by  a  change  in  particle 
geometry  or  orientation  is  correct,  such  differences  at  least  quali¬ 
tatively  should  be  observed  in  the  elastic  properties  of  the  compo¬ 
sites  as  well.  As  shown  in  Fig.  4,  no  such  effect  is  indicated  by 
the  data  for  Young's  modulus  [11]  for  all  the  specimens  with  values 
of  thermal  diffusivity  shown  in  Figs.  2  and  3.  This  represents 
positive  evidence  that  the  apparent  effect  of  the  Si  inclusion 
size  on  thermal  diffusivity  cannot  be  attributed  to  corresponding 
changes  in  inclusion  geometry  and/or  orientation. 

For  an  alternative  explanation  for  this  apparent  effect  of  in¬ 
clusion  size  it  should  be  noted  that  generally,  theories  for  the 
thermal  conductivity  of  composites  do  not  consider  the  effect  of 
an  interfacial  resistance  to  heat  flow.  One  possible  source  of 
such  an  interfacial  barrier  to  heat  flow  is  the  difference  in 
Young's  modulus  of  the  silicon  and  silicon  carbide  of  a  factor  of 
three.  The  interface  between  the  Si  and  SiC  in  these  composites 
represents  an  elastic  discontinuity  which  could  contribute  to 
increased  phonon  scattering.  At  room  temperature  such  a  discon¬ 
tinuity  may  affect  the  contribution  to  heat  flow  by  low  frequency 
phonons  primarily.  Such  an  effect  is  expected  to  be  more  pronounced 
for  the  smaller  Si  inclusions  than  for  the  larger  ones,  and  may 
decrease  the  thermal  diffusivity  at  room  temperature  by  perhaps  a 
fraction  of  a  percent,  but  does  not  explain  the  average  of  5  to  6 
percent  difference  between  the  two  sets  of  data  shown  in  Figs.  2 
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Fig.  A.  Young's  modulus  of  elasticity  at  room  temperature  for 
specimens  of  Figs.  2  and  3. 


and  3. 

The  most  likely  explanation,  albeit  qualitative  only  at  this 
time,  is  that  the  changes  in  Si  inclusion  size  are  accompanied  by 
corresponding  changes  in  the  thermal  diffusivity  of  the  silicon  and 
possibly  the  silicon  carbide  as  well.  Such  changes  in  the  thermal 
diffusivity  can  result  from  changes  in  impurity  content.  Proof 
for  the  validity  of  this  hypothesis  is  offered  by  the  large  dif¬ 
ference  in  the  value  of  the  thermal  diffusivity  of  0.668  cm^s-^ 
for  the  zone-refined  Si  and  the  corresponding  value  of  0.181 
cm^s-^  for  the  silicon  used  for  preparation  of  the  composite  speci¬ 
mens  . 

2  -1 

The  value  of  0.77A  cm  s  for  the  thermal  diffusivity  of  the 
sample  with  20.82  vol.%  Si  with  a  particle  size  of  A  pm  also  is 
relevant  to  the  above  discussion.  This  value  of  thermal  diffusi¬ 
vity  is  higher  than  the  corresponding  values  of  the  sintered  single¬ 
phase  a-SiC  as  well  as  the  zone-refined  and  industrial  grade  Si  in 
the  composites.  This  result  contradicts  composite  theory,  which 


states  that  Che  highest  value  of  the  thermal  conductivity  of  a  com¬ 
posite  is  governed  by  the  law  of  mixtures.  The  same  holds  for  the 
thermal  diffusivity  of  composites  with  components  of  comparable 
values  for  the  specific  heat  per  unit  volume.  If  it  is  assumed 
that  the  thermal  diffusivity  of  the  silicon  carbide  in  this  sample 
is  identical  to  the  value  for  the  single-phase  silicon  carbide,  then 
the  value  for  the  thermal  diffusivity  of  the  Si,  inferred  from  com¬ 
posite  theory,  would  be  far  in  excess  of  the  value  for  the  zone- 
refined  Si.  This  does  not  seem  likely.  A  more  plausible  ex¬ 
planation  is  that  the  thermal  diffusivity  of  the  SiC  in  this  SiC- 
Si  sample  exceeds  the  value  for  the  single  phase  silicon  carbide. 
This  latter  material  contains  a  higher  level  of  foreign  elements 
introduced  as  a  sintering  aid,  which  are  absent  in  SiC-Si  compo¬ 
sites.  It  is  speculated  here  that  for  high-purity  silicon  carbide, 
the  value  of  the  thermal  diffusivity  at  room  temperature  may  well 
approach  1  cm^s"^.  Regardless  of  the  actual  value,  it  seems  likely 
that  impurities  or  added  elements  in  both  the  silicon  and  silicon 
carbide  can  affect  the  heat  conduction  behavior  of  SiC-Si  composites 
significantly . 
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Fig.  5.  Temperature  dependence  of  the  thermal  diffusivity  of  two 
silicon  carbide-silicon  composites  for  two  values  of  si¬ 
licon  content. 


Figure  5  shows  the  temperature  dependence  of  the  thermal 
diffusivity  of  the  sample  with  20.82  vol.%  of  4  pm  Si  inclusions 
and  for  one  of  the  samples  of  Fig.  3.  These  data  indicate  that 
the  relative  difference  in  the  thermal  diffusivity  decreases  with 
increasing  temperature. 

Figure  6  compares  the  temperature  dependence  of  the  thermal 
diffusivity  of  two  of  the  samples  of  Figs.  2  and  3  with  comparable 
silicon  content.  These  data  are  of  interest  in  that  they  indicate 
that  above  700°C  the  values  for  the  thermal  diffusivity  are  almost 
identical.  It  is  suggested  here  that  this  effect  may  constitute 
evidence  for  an  electronic  contribution  to  the  heat  transfer  by 
thermally  activated  holes  and  electrons,  intrinsic  as  well  as  ex¬ 
trinsic  due  to  the  presence  of  the  various  impurities.  In  parti- 
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Fig.  6.  Temperature  dependence  of  thermal  diffusivity  of  silicon 
carbide-silicon  composites  with  comparable  silicon  con¬ 
tents  and  two  different  particle  sizes. 


cular,  this  would  be  the  case  for  the  sample  with  the  lower  value  of 
the  thermal  diffusivity  at  room  temperature,  presumably  because  of 
its  higher  impurity  content.  Indeed  if  this  hypothesis  is  correct, 
impurities  in  silicon  or  silicon  carbide  will  lower  the  thermal 
diffusivity  at  low  temperature,  but  may  increase  it  again  at  higher 
temperatures  if  these  impurities  are  of  the  type  which  affect  the 
electrical  conductivity.  The  higher  boron  content  in  the  silicon 
in  the  present  samples  could  especially  be  very  effective  in  this 
respect . 

Figure  7  shows  the  temperature  dependence  of  the  thermal  dif¬ 
fusivity  of  the  samples  of  zone-refined  and  industrial  silicon. 

For  these  samples  also,  the  relative  difference  in  thermal  dif¬ 
fusivity  decreases  with  increasing  temperature.  Of  interest  is 
the  irreversible  increase  in  the  thermal  diffusivity  of  the  in¬ 
dustrial  silicon  following  heating  and  cooling.  This  effect  could 
arise  from  the  formation  of  precipitates  of  impurties  originally 


Fig.  7. 


Temperature  dependence  of  the  thermal  diffusivity  of  an 
industrial  grade  and  purified  silicon. 


in  solid  solution.  No  evidence  was  obtained  in  support  of  this  hy¬ 
pothesis  . 


The  experimental  results  obtained  in  this  study  indicate  that 
the  heat  conduction  properties  of  SiC-Si  composites  can  be  affect¬ 
ed  by  a  number  of  variables.  The  evidence  presented  indicates  that 
impurities  are  particularly  effective  in  this  respect.  Of  partical 
significance  is  that  such  an  effect  would  allow  some  degree  of  con¬ 
trol  over  the  desired  thermal  conductivity  or  thermal  diffusivity 
of  such  composites.  High  impurity  contents  (in  the  absence  of  other 
adverse  effects,  of  course)  would  be  desirable  to  maximize  thermal 
insulating  properties.  On  the  other  hand,  components  or  structures 
subject  to  thermal  stress  failure  should  be  as  pure  as  possible  in 
order  to  obtain  high  thermal  conductivity. 
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Abstract 


The  variables  which  affect  the  formation  of  micro-cracks  due  to 
thermal  expansion  mismatch  in  brittle  composites  and  their  effect  on  heat 
conduction  properties  are  discussed.  Experimental  data  for  metal  oxide- 
silicon  carbide  composites  are  presented  which  show  that  microcrack  for¬ 
mation  can  lower  the  thermal  diffusivity  and  conductivity  significantly. 
Sub-critical  growth  and  healing  of  the  microcracks  during  thermal  cycling 
are  also  shown  to  have  an  effect  on  heat  transport.  Evidence  is  present¬ 
ed  that  microcrack  formation  in  brittle  composites  requires  the  presence 
of  a  residual  pore  phase. 


Introduction 


The  conduction  of  heat  through  solids  occurs  primarily  by  phonon, 
electron  or  photon  transport.  Fundamentally,  the  heat  flux  by  any  of 
these  transport  mechanisms  is  governed  by  wave-  and  quantum-mechanical 
principles  [Kittel  (1962);  Berman  (1976);  Kingery  et  al.  (1976)].  The 
heat  flux  is  affected  also  by  any  factor  which  contributes  to  phonon, 
electron  or  photon  scattering.  At  the  macro-structural  level  such 
scattering  can  occur  at  surfaces  and  interfaces.  At  the  atomic  level, 
the  presence  of  vacancies,  foreign  atoms,  dislocations,  radiation  damage, 
other  phonons  and/or  electrons,  absorption  centers,  any  other  structural 
or  chemical  imperfections,  and  elastic  or  optical  discontinuities  which 
promote  phonon,  electron  or  photon  scattering  tends  to  decrease  the  heat 
flux  significantly. 

At  the  microstructural  level,  heat  conduction  can  be  affected  by 
grain  boundaries,  intergranular  phases  and  the  existence  of  texture. 
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Inclusions  with  heat  conduction  properties  which  differ  from  those  of  the 
matrix,  coupled  with  the  corresponding  effect  of  their  geometry  and  orien¬ 
tation,  also  can  have  a  profound  effect  on  the  magnitude  of  heat  flux. 

Microcracks  represent  a  particular  type  of  microstruc tural  imperfec¬ 
tion.  For  purposes  of  this  paper,  microcracks  are  defined  as  cracks 
which  are  very  small  compared  to  the  size  of  the  sample  or  components  in 
which  they  are  contained.  The  effect  of  such  microcracks  on  heat  conduc¬ 
tion  properties  will  be  the  subject  of  the  present  paper,  with  special 
emphasis  on  microcracked  brittle  composites.  Such  materials  are  of  cur¬ 
rent  practical  interest  because  of  their  excellent  thermal  shock  resis¬ 
tance,  improved  fracture  toughness,  and  machinability  coupled  with  major 
decreases  in  elastic  moduli.  Most  importantly,  the  presence  of  the 
microcracks  changes  the  highly  unstable  failure  mode  typical  for  brittle 
materials  to  the  far  more  desirable  stable  mode  of  crack  propagation. 

These  improvements  in  properties  are  critical  for  structural  materials 
intended  for  use  at  high  temperatures. 

Microcracks  represent  barriers  to  phonon  and  electron  transport  and 
should  have  a  significant  effect  on  the  heat  conduction  properties  at 
levels  of  temperature  at  which  radiation  across  the  cracks  is  insignifi¬ 
cant.  The  specific  purpose  of  this  paper  is  to  discuss  the  variables 
which  affect  microcracking  and  their  corresponding  effect  on  thermal 
transport  properties  of  brittle  composites.  Experimental  data  are  pre¬ 
sented  in  support  of  the  findings  obtained. 


Analysis 

Effect  of  Cracks  on  Thermal  Conductivity 

Equations  for  the  effect  of  cracks  on  thermal  conductivity  were  de¬ 
rived  by  Hasselman  (1978)  from  the  solutions  of  Fricke  (1924)  and  Powers 
(1961)  for  the  thermal  conductivity  of  composites  consisting  of  a  con¬ 
tinuous  matrix  with  dilute  dispersions  in  the  shape  of  oblate  ellipsoids 
of  revolution  as  follows: 

1.  Randomly  oriented  cracks: 

K  =  K  (1  +  SNU3/^)-1  (1) 

o 

where  K  and  K0  are  the  thermal  conductivity  of  the  microcracked  and  crack- 
free  solids,  respectively,  N  is  the  number  of  cracks  per  unit  volume  and 
l  is  the  crack  radius. 

2.  Cracks  oriented  perpendicularly  to  the  direction  of  heat  flow: 

K  =  K  (1  +  8Nf3/3) 
o 

3.  Cracks  parallel  to  the  heat  flow: 


(2) 
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Comparison  of  eqs.  1,  2  and  3  shows  that  cracks  oriented  perpendicu¬ 
lar  to  the  heat  flow  are  most  effective  in  decreasing  the  thermal  conduc¬ 
tivity,  followed  by  randomly  oriented  cracks  and  with  no  effect  for  cracks 
oriented  parallel  to  the  direction  of  heat  flow. 


Variables  Which  Affect  Microcracking  in  Brittle  Composites 

A  mismatch  in  the  coefficients  of  thermal  expansion  of  the  individual 
components  is  the  principal  cause  of  microcrack  formation  in  brittle  com¬ 
posites.  A  thermal  expansion  mismatch  results  in  the  generation  of  inter¬ 
nal  stresses  whenever  the  temperature  of  the  composite  is  changed  from 
the  value  at  which  such  internal  stresses  are  zero.  This  latter  value 
normally  corresponds  to  the  manufacturing  temperature.  In  brittle  com¬ 
posites,  such  internal  stresses  can  reach  high  magnitude  because  internal 
stress  relief  by  plastic  (i.o.  dislocation)  flow  cannot  occur.  Stress 
relief  could  occur  by  diffusions!  processes,  but  only  at  high  tempera¬ 
tures.  For  practical  purposes,  brittle  composites  exhibit  elastic  be¬ 
havior  with  an  internal  stress  magnitude  directly  proportional  to  the 
range  of  temperature  over  which  t lie  composite  is  cooled. 

The  internal  stresses  can  lead  to  the  generation  of  microcracks  when¬ 
ever  their  magnitude  exceeds  the  fracture  stress  of  either  component  or 
the  cohesive  strength  of  the  interface  between  the  components.  It  is 
generally  recognized  that  microcracks  are  initiated  from  pre-existing 
flaws  in  the  form  of  small  pores  or  other  crack-like  defects  within  any 
of  the  components  or  at  the  interface.  For  this  reason,  an  analysis  of 
microcrack  formation  should  be  based  on  fracture-mechanical  principles. 

In  view  of  the  generally  non-uniform  distribution  of  internal  stresses 
which  results  from  a  thermal  expansion  mismatch,  numerical  techniques 
based  on  finite  element  methods  are  required  (Evans  and  Clark,  1980) .  An 
exception  to  this  is  a  spherical  inclusion  contained  within  an  infinite 
matrix.  For  such  an  inclusion,  the  state  of  internal  stress  due  to  ther¬ 
mal  expansion  mismatch  is  uniform.  By  permitting  microcrack  formation 
within  the  inclusion  only  and  treating  the  relevant  properties  affected 
by  the  microcracks  as  pseudo-continuum  properties,  analytical  solutions, 
as  shown  by  Niihara  (1980),  are  easily  obtained.  From  such  solutions, 
the  relevant  varibles  which  affect  microcracking  are  easily  assessed,  as 
follows: 

For  a  spherical  inclusion,  the  total  elastic  energy  (Wg)  of  the  re¬ 
sidual  stress  field  in  both  the  inclusion  and  matrix  due  to  the  mismatch 
in  the  coefficients  of  thermal  expansion  can  be  obtained  from  the  solu¬ 
tion  of  Davidge  and  Green  (1968) . 

Wv  =  2tt(A«AT)2R3/[(1-2v  )/E  +  (Hv  )/2E  ]  (4) 

£.  p  p  m  m 

where  Aa  is  the  mismatch  in  the  coefficients  of  thermal  expansion  of  the 
inclusion  and  the  matrix,  AT  is  the  temperature  difference  over  which  the 
composite  is  cooled  from  the  value  at  which  the  internal  stresses  are 
zero  (assuming  linearly  elastic  behavior  throughout),  R  is  the  radius 
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of  the  inclusion,  v  is  Poisson’s  ratio,  E  is  Young's  modulus  and  p  and  m 
refer  to  the  inclusion  and  matrix,  respectively. 

For  randomly  oriented  penny-shaped  microcracks  of  equal  size,  Young's 
modulus  of  the  inclusion  can  be  expressed  (Salganik,  1973) : 

E  =  E  [l+16(l-v  2)N£3/9(l-2v  )  ]_1  (5) 

p  o  p  p 

where  N  is  the  number  of  microcracks  per  unit  volume  of  inclusion  and  E^ 
is  Young's  modulus  of  the  crack-free  inclusion. 

The  creation  of  the  microcracks  within  the  inclusion,  requires  an 
expenditure  of  fracture  energy  (Wp) : 

W  =  4irR3Ny  it22/3  (6) 

F  i 

where  Yf  is  the  total  energy  dissipated  in  creating  unit  area  of  fracture 
surface. 

Following  the  approach  of  Griffith  (1924),  crack  instability  occurs 

when: 

d/dX.(WE  -  Wp)  -  0  (7) 


which  yields  the  critical  value  of  temperature  (ATc)  required  for  micro¬ 
crack  instability: 


AT 

c 


7ryf(l-2vp) 

2  2  2 
(AarE  (1-v  ) 

o  p 


h 


E  (1+v  )  16(l-v  2)N£3 

_o - ! - +1+  - E - 

2E  (l-2v  )  9(l-2v  ) 

m  p  p  _ 


(8) 


Alternatively,  eq.  1  can  be  rewritten  in  terms  of  the  minimum  size 
of  precursor  microcrack  required  for  microcrack  formation  for  a  given 
value  of  AT. 

Cracks  can  also  propagate  in  a  subcritical  mode  for  Ys  <  Y  <  Yf» 
where  Ys  is  the  thermodynamic  surface  free  energy.  The  minimum  value  of 
AT  required  for  subcritical  crack  propagation  (ATSC)  can  be  obtained  from 
eq .  8  by  substitution  of  Ys  for  Yf-  For  values  of  AT  <  ATSC  crack  heal¬ 
ing  will  occur.  Of  course,  such  subcritical  crack  growth  and  healing  will 
take  place  only  under  environmental  conditions  and  levels  of  temperature 
at  which  the  appropriate  reaction  rates  of  diffusional  processes  are 
high  enough. 

Eq.  8  shows  that  the  microcrack  is  unstable  between  a  lower  and  up¬ 
per  value  of  crack  length.  This  implies  that  once  a  microcrack  is  formed 
by  propagation  of  a  precursor  flaw,  it  will  arrest  at  the  higher  value  of 
crack  length.  It  should  be  recognized,  however,  that  for  small  precursor 
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microcracks,  at  instability  at  AT  =  ATC,  the  strain  energy  release  rate 
exceeds  the  fracture  surface  energy.  This  causes  the  crack  to  propagate 
in  an  unstable  mode.  The  crack  will  have  maximum  kinetic  energy  when  it 
achieves  the  upper  value  of  crack  length  for  stable  crack  propagation. 
This  will  result  in  an  "overshoot",  with  crack  arrest  following  only  at  a 
value  of  crack  length  ( Jlf )  when  all  the  released  potential  energy  is 
transformated  into  fracture  surface  energy  which  can  be  derived  to  be: 

lc  =  [ 3E  (1+v  )/16  E  (1-v  )N*.  (9) 

f  L  o  m  m  p  o 

Because  of  the  "overshoot",  a  crack  with  length  If  will  be  subcriti- 
cal  with  respect  to  the  value  of  ATC  at  which  it  formed  from  the  pre¬ 
cursor  microcrack.  Because  of  their  subcritical  nature,  such  microcracks 
may  exhibit  crack-healing,  unless  the  temperature  difference  was  raised 
to  a  value  well  in  excess  of  ATC. 

Equations  8  and  9  for  ATC,  ATSC  and  if  can  be  illustrated  graphical¬ 
ly.  Unfortunately,  because  of  the  form  of  these  equations  this  cannot  be 
done  in  a  convenient  non-dimensional  form,  but  requires  the  choice  of  a 
specific  combination  of  inclusion  and  matrix. 

For  a  MgO  inclusion  within  a  SiC  matrix,  Fig.  1  shows  ATC,  ATSC  and 
if  for  an  arbitrary  value  of  crack  density  N  =  108cm-3  and  Yf  and  Ysc 
equal  to  2  and  0.5  J.m-2,  respectively.  The  other  pertinent  property 
data  are  listed  in  Table  1  given  in  the  following  chapter  devoted  to  ex¬ 
perimental  data.  In  Fig.  1,  the  upper  and  lower  solid  curves  correspond 
to  ATC  and  ATSC  for  critical  and  subcritical  crack  propagation,  respec¬ 
tively.  The  dotted  curve  represents  the  crack  length  (fcf)  which  results 
from  unstable  crack  propagation  given  by  eq.  9.  The  regions  marked  A,  B 
and  C  represent  regimes  of  values  of  AT  and  crack  length  in  which 


Fig.  1.  Stability  diagram  of  microcracks  in  MgO 
inclusion  within  SiC  matrix. 


374 


L.  D.  BENTSEN,  D.  P.  H.  HASSELMAN,  N .  CLAUSSEN 


the  microcracks  exhibit  critical  crack  propagation,  subcritical  crack 
propagation,  and  crack  healing,  respectively.  Included  in  Fig.  1  are 
paths  of  propagation  of  microcracks  with  three  different  initial  sizes. 
The  microcrack  with  initial  size  propagates  entirely  in  a  subcritical 
mode.  For  the  initial  size  the  crack  propagates  in  an  unstable  mode 
followed  by  subcritical  crack  growth.  For  a  very  short  initial  crack 
length,  the  crack  propagates  in  an  unstable  mode  followed  by  crack 
healing.  More  complex  paths,  which  combine  both  subcritical  and  critical 
crack  propagation  can  occur  as  well,  depending  on  the  rate  of  temperature 
change  and  the  presence  or  absence  of  a  mechanism  for  subcritical  crack 
growth. 

It  is  important  to  conclude  from  the  above  derivation  and  discussion 
that  a  unique  value  of  crack  length  is  associated  with  each  value  of  AT. 
For  this  reason  any  changes  in  AT  will  require  a  corresponding  change  in 
crack  length.  For  this  reason,  it  is  expected  that  microcracked  compo¬ 
sites  under  non-isothermal  conditions  will  exhibit  time-dependent  behavior 
of  all  properties  affected  by  the  presence  of  the  microcracks. 

As  a  final  remark,  it  should  be  noted,  as  mentioned  earlier,  that  a 
precursor  such  as  a  pore  at  a  triple-point  or  grain  boundary  is  required 
for  microcrack  formation.  For  this  reason,  it  is  anticipated  that  the 
amount  of  residual  pore  phase  due  to  lack  of  complete  densification  may 
play  a  vital  role  in  affecting  the  microcrack  density  and  resulting  heat 
conduction  properties.  Experimental  evidence  for  this  effect  will  be 
presented . 

In  summary  then  the  above  discussion  indicates  that  the  extent  of 
microcracking  in  brittle  composites  is  governed  by  the  thermal  expansion 
mismatch  and  the  elastic  properties  of  the  components,  the  amount  of  pore 
phase,  the  level  of  temperature,  and  any  time-dependent  subcritical  crack 
growth  and  healing. 


The  Thermal  Conductivity  of  Microcracked  Composites 

Many  expressions  are  available  from  the  literature  for  the  thermal 
conductivity  of  composites  for  many  types  of  distributions  of  the  indi¬ 
vidual  components  within  the  composites. 

For  a  matrix  with  dilute  concentrations  of  spherical  dispersions 
the  thermal  conductivity  (Kc)  can  be  expressed,  (Raleigh,  1892;  Maxwell, 
1904): 


K  =  K 
c  m 


2K  +  K  -  2  V  (K  -K  ) 
m  p  p  m  p 


2K  +  K  +  V  (K  -K  ) 

m  p  p  m  p 


(10) 


where  K  is  the  thermal  conductivity,  Vp  is  the  volume  fraction  of  dis¬ 
persed  phase  and  the  subscripts  m  and  p  refer  to  the  matrix  and  dispersed 
phase,  respectively. 
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The  effect  of  microcracking  in  either  the  matrix  or  dispersed  phase 

or  both  on  the  thermal  conductivity  of  the  composite  can  be  accounted  for 

by  the  appropriate  substitution  of  eq.  1  in  the  case  of  randomly  oriented 

microcracks  for  K  and  K  in  eq.  10. 

m  p  ^ 

In  general,  as  decided  by  the  direction  of  the  thermal  expansion 
mismatch,  microcrack  formation  will  occur  in  either  the  matrix  or  the 
dispersed  phase.  For  this  reason,  the  composite  thermal  conductivity  as 
expressed  by  eq .  10  will  be  governed  by  changes  in  Kp  or  K^,  only.  An 
expression  for  the  thermal  conductivity  of  a  dispersed  phase  wi  t'n  randomly 
oriented  microcracks  can  be  obtained  readily  by  substitution  of  the  value 
of  Nb^  obtained  from  eqs.  8  or  9  into  eq.  1,  which  in  turn,  upon  substi¬ 
tution  into  eq.  10  yields  the  value  for  the  thermal  conductivity  of  the 
composite.  This  should  result  in  a  reasonably  reliable  estimate. 

Microcracking  in  the  matrix  occurs  in  the  immediate  vicinity  of  the 
inclusions.  For  this  reason,  the  simultaneous  effect  of  the  inclusions 
on  the  thermal  conductivity  of  the  composite  and  the  associated  micro¬ 
cracking  in  the  matrix  are  expected  to  be  coupled.  Therefore,  a  simple 
substitution  of  the  value  of  the  thermal  conductivity  of  the  microcracked 
matrix  phase  for  Km  in  eq.  10,  possibly  may  not  yield  quantitatively  re¬ 
liable  results.  An  analysis  of  such  a  coupling  effect  could  be  of  inter¬ 
est  to  the  reader  of  this  paper. 

In  general,  the  above  analysis  and  discussion  shows  that  the  thermal 
conductivity  of  microcracked  brittle  composites  is  expected  to  be  a  func¬ 
tion  of  the  volume  fraction,  thermal  expansion,  and  elastic  behavior  of 
the  individual  components  as  well  as  their  crack  propagation  characteris¬ 
tics  including  critical  and  time-dependent  subcritical  propagation  as  well 
as  crack  healing.  The  purpose  of  this  paper  is  to  present  experimental 
data  which  illustrate  a  number  of  these  effects. 


Materials  and  Experimental  Procedures 

The  composite  materials  for  this  study  consisted  of  continuous 
matrices  of  aluminum  oxide,  beryllium  oxide  and  magnesium  oxide  contain¬ 
ing  a  range  of  volume  fractions  of  a  dispersed  phase  of  silicon  carbide. 
The  SiC  particle  size  was  of  the  order  of  a  few  microns.  These  compo¬ 
sites  were  prepared  by  vacuum-hot-pressing  mixtures  of  the  appropriate 
powders  at  temperatures  and  over  time  periods  required  to  achieve  nearly 
full  densif ication.  Samples  for  the  BeO-SiC  and  MgO-SiC  composites  were 
obtained  from  a  commercial  source.*  The  Al203~SiC  composites  were  made 
in  the  laboratory  of  one  of  the  present  writers. §  Table  1  lists  the 
approximate  values  for  the  coefficients  of  thermal  expansion  and  elastic 
properties  of  the  four  individual  components.  The  direction  of  the  mis¬ 
match  in  the  coefficients  of  expansion  of  the  matrix  and  dispersed  SiC 


*Ceradyne  Corporation,  Santa  Ana,  California. 
§Max-Planck  Institute,  Stuttgart,  FRG. 
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phase  is  such  that  on  cooling  from  the  hot-pressing  temperature,  micro¬ 
crack  formation  is  expected  to  occur  within  the  oxide  matrix.  Scanning 
electron  micrography  revealed  that  such  microcracking  occurred  primarily 
by  inter-granular  fracture. 


Table  1  Approximate  Property  Values  for  Components 
of  Silicon  Carbide-Metal  Oxide  Composites 


Property 

SiC 

A12°3 

BeO 

MgO 

Coefficient  of  thermal  expansion 
OC_1(xl06) 

5 

7 

10 

13 

Young's  modulus  (lO^MPa) 

4.5 

4.1 

3.0 

3.0 

Poisson's  ratio 

0.18 

0.26 

0.34 

0.30 

The  heat  conduction  behavior  of  the  composites  was  established  by 
measuring  their  thermal  diffusivity  by  the  laser-flash  method  [Parker,  et 
al.,  196l]  using  specimens  cut  to  the  appropriate  geometry  and  dimensions. 
Details  of  the  equipment  were  described  earlier  [siebeneck,  et  al.,  1976]. 
Since  microcracking  has  a  negligible  effect  on  the  heat  capacity  per  unit 
volume,  the  data  for  the  thermal  diffusivity  closely  reflect  the  corres¬ 
ponding  behavior  of  the  thermal  conductivity . 

Experimental  Results  and  Discussion 

Figure  2  shows  the  experimental  data  for  the  thermal  diffusivity  at 
room  temperature  for  the  Al203~SiC  composites.  All  samples  for  these 
data  were  in  the  as-received  condition  and  were  subjected  to  cooling 
from  the  hot-pressing  temperature  only,  without  any  additional  heating 
or  cooling.  Included  in  Fig.  2  is  the  dependence  of  the  thermal  diffu¬ 
sivity  on  silicon  carbide  content  for  the  crack-free  composites  calculat¬ 
ed  from  the  theory  for  the  thermal  conductivity  of  composites  given  by 
eq.  10  and  the  corresponding  changes  in  density  and  specific  heat.  It 
can  be  noted  that  the  experimental  data  fall  far  below  those  calculated 
for  the  crack-free  material.  The  relative  deviation  increases  with  in¬ 
creasing  SiC  content,  as  expected,  since  the  microcrack  density  is  ex¬ 
pected  to  be  proportional  to  the  number  of  SiC  particles  per  unit  volume. 
The  data  also  indicate  that  for  this  composite  system,  microcracking  re¬ 
duces  the  thermal  diffusivity  (and  therefore  the  thermal  conductivity)  by 
ar  much  as  a  factor  of  two. 
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SiC  CONTENT  (wt  %) 


Fig.  2.  Thermal  diffusivity  of 
A^O^-SiC  composites. 


SiC  CONTENT  (wt  %) 

Fig.  3.  Thermal  diffusivity  of 
BeO-SiC  composites. 


Figure  3  shows  the  data  for  the  thermal  diffusivity  of  the  BeO-SiC 
samples.  The  specimens  for  these  measurements  as  for  the  AI2O3-S1C  data, 
were  in  the  as-received  condition.  These  data  indicate  that  the  thermal 
diffusivity  does  not  decrease  linearly  with  SiC  content,  as  indicated  by 
an  increase  in  the  thermal  diffusivity  at  50  wt.%  SiC.  A  similar  but 
much  greater  increase  was  observed  also  in  a  previous  study  [Singh,  et 
al ,  198l]  which  suggested  a  near  complete  suppression  of  micro-crack  for¬ 
mation  near  50  wt.%  SiC.  It  is  of  interest  to  point  out  that  experimen¬ 
tal  data  for  the  fracture  toughness  and  strength  of  BeO-SiC  composites, 
obtained  by  others  [Rubin,  1979],  show  a  similar  suppression  of  micro¬ 
cracking  near  50%  SiC.  Near  20  and  80%  SiC,  toughness  and  strength 
exhibit  maxima  due  to  microcrack  toughening,  whereas  no  such  toughening 
appears  to  exist  near  50%  SiC.  These  independent  observations  indicate 
that  they  are  not  spurious  and  confirm  that  the  decrease  in  the  thermal 
diffusivity  observed  in  the  present  study  indeed  are  caused  by  micro- 
cracking. 

The  observations  on  the  MgO-SiC  composites  also  were  informative  in 
this  respect.  Figure  4  shows  the  thermal  diffusivity  as  a  function  of 
SiC  content  for  a  set  of  samples,  designated  Set  I,  investigated  earlier. 
Figure  5  shows  the  thermal  diffusivity  of  MgO-SiC  composites,  designated 
Set  II,  obtained  more  recently.  For  this  latter  set  the  data  shown  re¬ 
present  the  average  for  5  samples  at  a  given  SiC  content. 
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Fig.  4.  Thermal  diffusivity  of  set 
I  of  MgO-SiC  composites. 


Fig.  5.  Thermal  diffusivity  of  set 
II  of  MgO-SiC  composites. 


Comparison  of  the  data  in  Fig.  4  labelled  "before  thermal  cycle"  (i. 
e.  immediately  after  cooling  from  the  processing  temperature)  with  those 
in  Fig.  5,  indicates  that  for  the  MgO-SiC  samples  considerable  differences 
in  thermal  diffusivity  can  exist  even  at  the  same  SiC  content.  Comparing 
the  data  for  all  three  types  of  composites  shows  that  the  largest  relative- 
reductions  are  observed  for  the  MgO-SiC  series.  This  is  expected,  since 
for  this  composite  system,  the  mismatch  in  coefficients  of  therma]  expan¬ 
sion  is  higher  than  for  the  other  two  series.  Comparison  of  all  four  sets 
of  data  shows  that  the  dependence  of  the  thermal  diffusivity  on  silicon 
carbide  content  differs  not  only  quantitatively  but  also  qualitatively 
from  composite  series  to  series.  The  suppression  of  microcracking  at  the 
intermediate  values  of  SiC  content  was  attributed  previously  to  differ¬ 
ences  in  the  nature  of  the  distribution  of  the  components  within  the  com¬ 
posite  at  low,  high  and  intermediate  contents  of  silicon  carbide  [Singh, 
198l].  Such  differences  not  only  give  rise  to  differences  in  the  magni¬ 
tude  of  the  internal  stress,  but  also  affect  the  probability  of  microcrack 
formation . 

Such  effects  undoubtedly  are  still  operative.  However,  with  the  ac¬ 
cumulation  of  additional  data  for  the  other  composite  series  an  additional 
explanation  can  be  given,  based  on  differences  in  residual  pore  content 
from  sample  to  sample,  which  also  accounts  for  the  qualitative  differ¬ 
ences  between  the  different  composite  and  sample  series. 


Evidence  for  the  val idity  of  this  hypothesis  is  presented  in  Figs. 
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Fig.  6.  Thermal  diffusivity  of 

micro-cracked  composites 
relative  to  crack-free 
values  as  a  function  of 
pore  content,  a.  AI2O3- 
SiC;  b.  BeO-SiC  and  c. 
MgO-SiC. 


c  0  12  3  4  5 

%  POROSITY 

6a,  b  and  c  which  show  the  value  of  the  thermal  diffusivity  of  the  AI2O3 
BeO  and  MgO-based  composites,  respectively,  relative  to  the  crack-free 
value  plotted  as  a  function  of  pore  content.  It  can  be  noted  that  for 
all  three  series  the  relative  thermal  diffusivity,  albeit  with  some 
scatter,  decreases  very  rapidly  with  increasing  pore  content.  In  inter¬ 
preting  this  effect,  it  should  be  noted  that  this  inverse  correlation 
between  thermal  diffusivity  and  pore  content  cannot  be  regarded  as  a....-  "" 
direct  cause  and  effect.  Any  decrease  in  density  due  to  the  presence 
of  the  microcracks  is  less  than  a  fraction  of  one  percent.  For  this  rea¬ 
son  the  values  of  pore  content  shown  in  Fig.  6  reflect  the  lack  of  com¬ 
plete  densif ication  during  the  hot-pressing  operation.  The  resulting 
residual  pore  phase  can  be  regarded  as  a  special  case  of  spherical  in¬ 
clusions  with  zero  thermal  conductivity.  For  such  a  dispersed  phase, 
eq.  10  shows  that  at  low  volume  fractions,  each  percent  of  pore  content 
reduces  the  thermal  conductivity  by  about  1.5%.  The  specific  heat  per 
unit  volume  also  is  reduced  by  a  fractional  amount  equal  to  the  fractional 
pore  volume.  These  two  effects  combined  result  in  a  net  decrease  in  the 


380 


L.  D.  BENTSEN,  D.  P.  H.  HASSEI.MAN,  N.  CLAUSSEN 


thermal  diffusivity  of  0.5%  for  each  percent  of  pore  phase.  This  value  is 
less  than  the  observed  decreases  shown  in  Figs.  6a,  b  and  c  by  well  in 
excess  of  an  order  of  magnitude. 

It  is  suggested  that  the  data  shown  in  Fig.  6  constitute  evidence  for 
the  existence  for  pore-assisted  microcrack  formation.  As  discussed  ear¬ 
lier,  for  microcrack  precursors  in  the  form  of  pores  at  triple  points  or 
grain-boundaries  which  are  too  small,  microcracks  will  not  form  regardless 
of  the  magnitude  of  internal  stress.  For  such  small  pores,  the  corres¬ 
ponding  composites  will  exhibit  near  theoretical  density  and  values  of 
thermal  diffusivity  near  the  value  for  the  crack-free  material.  For  a 
given  number  of  pores  per  unit  volume,  an  increase  in  pore  content  is  ac¬ 
companied  by  an  increase  in  pore  size.  This  in  turn,  would  promote  micro¬ 
crack  formation  and  a  corresponding  decrease  in  thermal  diffusivity,  as 
observed.  It  is  expected  that  for  a  given  number  of  pores  per  unit 
volume,  some  minimum  ("threshold")  value  of  pore  content  would  be  required 
before  microcracks  can  be  generated.  Unfortunately,  the  number  of  data 
points  shown  in  Fig.  6  is  insufficient  to  unequivocally  demonstrate  the 
existence  of  such  "threshold"  pore  content.  The  finding  of  such  a  "thres¬ 
hold"  value  also  is  expected  to  be  complicated  by  the  usual  existence  of  a 
pore-size  distribution.  It  should  be  noted  also  that  the  effect  of  pore 
phase  on  the  microcracking  and  the  thermal  diffusivity  for  these  compo¬ 
sites  is  coupled  with  the  content  of  silicon  carbide,  which  itself  may 
lead  to  an  increase  in  pore  content  due  to  the  increasing  degree  of  diffi¬ 
culty  of  densification.  The  role  of  the  pore  phase  on  microcracking  can 
be  ascertained  in  more  quantitative  detail  for  specimens  of  constant  SIC 
content.  These  specimens  are  being  obtained  with  results  to  be  presented 
in  a  future  report. 

Regardless  of  the  quantitative  details,  the  data  in  Fig.  6  indicate 
that  by  careful  control  over  pore  content  achieved  by  corresponding  con¬ 
trol  over  processing  conditions,  the  properties  of  composites  affected  by 
the  presence  of  microcracks  may  be  tailored  to  specific  values  required 
for  specific  engineering  applications. 

Heat  treatments,  in  addition  to  the  cooling  from  the  hot-pressing 
temperature,  also  have  a  significant  effect  on  the  thermal  diffusivity. 
Experimental  evidence  for  this  effect  is  shown  in  Fig.  4  which  shows  data 
for  the  thermal  diffusivity  before  and  after  a  thermal  cycle  to  1400°C 
over  a  time  period  of  approximately  8  hours.  Comparison  of  the  data  shows 
that  this  thermal  cycle  resulted  in  a  general  decrease  in  the  thermal  dif¬ 
fusivity.  Previously,  samples  of  MgO  +  20  wt.%  SiC  held  at  1300°C  for  0, 

3  and  10  hours  showed  decreases  In  the  thermal  diffusivity  on  return  to 
room  temperature  of  8.1,  11.3  and  15.4%,  respectively.  These  observations 
are  thought  to  constitute  evidence  of  the  subcritical  growth  of  microcrack 
precursors  which  were  not  of  sufficient  size  to  result  in  microcracks 
during  the  initial  cooling  from  the  hot-pressing  temperature.  The  addi¬ 
tional  decrease  in  the  thermal  diffusivity  following  thermal  cycles  is  the 
result  of  additional  microcrack  formation.  The  above  also  clearly  indi¬ 
cated,  as  discussed  earlier,  that  microcracked  materials  are  non-equili¬ 
brium  materials.  In  fact.  It  can  be  said  that  microcracks  can  be  super¬ 
cooled  . 
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Fig.  7.  Effect  of  thermal  cycling  Fig.  8.  Effect  of  thermal  cycling 
on  thermal  diffusivity  of  on  thermal  diffusivity  of 

MgO  +  20  wt.%  SiC.  MgO  +  30  wt.%  SiC. 


Evidence  for  crack  healing  was  also  found.  Figure  7  shows  the  ex¬ 
perimental  data  for  the  thermal  diffusivity  of  a  specimen  of  MgO  +  2071 
SiC  during  heating  and  cooling  and  holding  for  3  hours  at  1300°C.  It  may 
be  noted  that  at  the  higher  temperatures  on  cooling,  the  values  for  the 
thermal  diffusivity  are  in  excess  of  those  observed  during  heat-up. 

It  is  thought  that  this  effect  arises  from  the  healing  of  microcracks 
which  formed  during  the  initial  cooling.  On  return  to  room  temperature 
however  the  thermal  diffusivity  after  the  thermal  cycle  and  hold  is  nov 
less  than  before.  This,  as  discussed  earlier,  is  evidence  for  addition 
microcrack  formation  due  to  growth  of  precursors  initially  too  small.  it 
appears  then  that  the  changes  in  the  thermal  diffusivity  on  thermal  cycl¬ 
ing  can  be  affected  by  the  simultaneous  healing  of  microcracks  and  prowth 
of  their  precursors. 

Figure  8  shows  additional  evidence  for  this  effect  occurring  in  a 
sample  of  MgO  +  30%  SiC  cycled  to  1400°C  over  a  period  of  some  8  hrs. 
These  data  show  clear  evidence  of  pop-in  of  microcracks  once  the  compo¬ 
site  is  cooled  over  a  temperature  range  such  that  the  stresses  become 
sufficiently  large  that  either  healed  microcracks  can  regenerate  or  that 
precursor  cracks  have  reached  sufficient  size  to  form  new  microcracks. 


Regardless  of  the  underlying  details,  the  data  presented  in  this 
paper  clearly  indicate  that  the  properties  of  microcracked  materials  are 
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not  "hand-book  values",  but  are  strong  functions  of  past  processing  and 
thermal  history.  The  existence  of  such  variations  must  be  recognized  in 
the  applications  of  such  materials  for  engineering  purposes. 
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ABSTRACT 

It  is  suggested  that  the  thermal  conductivity  of  very  fine  fibers  can  be 
evaluated  indirectly  with  the  aid  of  composite  theory  using  the  experimental 
data  for  the  heat  transport  properties  of  an  appropriate  composite  which 
contains  the  fibers.  The  feasibility  of  this  approach  was  investigated 
by  determining  the  thermal  conductivity  and  diffusivity  of  fibers  of  amorphous 
silicon  carbide  from  25°C  to  1000°C  contained  within  a  lithium  aluminosilicate 
glass-ceramic  using  the  laser-flash  technique  for  measurement  of  the  thermal 
diffusivity  of  the  composite.  Due  to  the  amorphous  nature  of  the  fibers, 
their  values  for  the  thermal  conductivity  and  diffusivity  were  found  to  be  far 
less  than  the  corresponding  data  for  crystalline  silicon  carbide.  The  positive 
temperature  dependence  of  the  thermal  conductivity  coupled  with  the  independent 
observation  of  an  increase  in  thermal  conductivity  with  specimen  thickness 
suggests  that  radiative  heat  transfer  makes  a  significant  contribution  to  the 
total  heat  transferred.  A  number  of  advantages  and  limitations  of  the  composite 
method  for  the  evaluation  of  thermal  transport  properties  of  fibers  are  dis¬ 


cussed. 


1. 


INTRODUCTION 


The  technical  and  scientific  need  for  quantitative  data  for  the  heat 
conduction  properties  of  solids  has  led  to  the  development  of  nearly  in¬ 
numerable  methods  for  the  measurement  of  these  properties  [l],  For  a.  given 
material  the  choice  of  the  most  appropriate  method  depends  on  a  number  of 
variables  such  as  the  temperature  range  over  which  the  data  are  required  and 
the  configuration  of  the  material  to  be  measured.  In  this  respect,  micron¬ 
sized  fibers  present  special  difficulties. 

Materials  in  fiber  form  are  used  extensively  in  the  development  of 
mechanically  reinforced  composites.  One  such  composite  developed  during 
the  last  few  years  consists  of  a  glass  or  glass-ceramic  reinforced  with  fibers 
of  silicon  carbide  [2].  The  mechanical  properties  of  these  composites 
indicate  excellent  service  performance  at  elevated  temperature.  For  a  micro- 
mechanical  anlaysis  of  these  composites  under  conditions  of  transient  heat 
flow,  experimental  data  for  the  thermal  conductivity  and  diffusivity  were 
required  as  a  function  of  temperature.  The  choice  of  the  method  for  the 
measurement  of  these  properties  presented  a  number  of  difficulties.  Firstly, 
their  method  of  preparation  precluded  the  synthesis  of  larger  samples  with 
properties  identical  to  those  of  the  fibers.  For  this  reason,  only  the 
fibers  themselves  had  to  be  used.  Secondly,  the  silicon  carbide  fibers  were 
not  electrical  conductors.  This  eliminated  the  resistive  heating  method  [3] 
for  measuring  the  thermal  conductivity.  Although,  in  principle,  a  number  of 
techniques  could  be  devised,  it  was  decided  to  measure  the  thermal  conductivity 
and  diffusivity  of  these  fibers  by  an  indirect  method.  This  method  consisted 
of  measuring  the  thermal  diffusivity  of  the  composite  as  well  as  the  matrix 
(without  fibers),  followed  by  calculation  of  the  required  data  by  means  of 
the  theory  for  the  thermal  conductivity  of  composites.  It  is  the  purpose  of 


this  paper  to  report  the  data  obtained. 

2 .  EXPERIMENTAL 
2.1  Materials 

The  silicon  carbide  fibers  made  from  an  organrmetallic  polymer 
by  the  method  of  Yajima  [4]  were  obtained  from  a  commercial  source. 

Chemical  analysis  showed  the  fibers  to  consist  of  approximately  65%  SiC, 

25%  Si02  and  about  10%  C.  The  fibers  were  in  the  form  of  tows  of  yarn 
containing  ^500  fibers  per  tow.  The  average  fiber  diameter  was  approxi- 
mately  10  pm  and  the  density  was  approximately  2.55  g/cm  .  X-ray  analysis 
showed  the  fibers  to  be  amorphous  with  a  crystallite  size  of  25-30  X 
determined  from  peak  broadening. 

3 

The  matrix  material  with  density  of  2.52  g/cm  consisted  of  a  lithium 
aluminosilicate  glass-ceramic  essentially  identical  to  that  of  a  commercial 
glass-ceramic^  except  that  the  ^3  wt.%  Ti02  nucleating  agent  was  replaced 
by  ^2  wt.%  Zr02  due  to  reactivity  between  the  Ti02  and  the  SiC  fibers. 

The  composite  was  manufactured  by  passing  the  SiC  yarn  through  a  slurry 
of  powder  of  the  uncrystallized  glass  and  isopropyl  alcohol.  Following 
drying,  the  coated  yam  was  cut  to  fit  the  graphite  die  for  hot-pressing 
in  vacuum  for  times  of  5  min.  to  1  hr.  at  1400  to  1500°C  and  a  pressure  of 
14  MPa.  Following  hot-pressing,  crystallization  of  the  glass  was  carried  out 
by  heat  treating  for  1  to  2  hr  at  a  temperature  ranging  from  880  to  1100°C. 
Composite  samples  were  made  in  which  the  fibers  were  all  parallel  or  with 
alternate  layers  oriented  at  90°,  Figure  1  shows  a  photomicrograph  of  such 
a  0/90  cross-ply  composite.  A  sample  of  the  matrix  phase  without  fibers 
was  also  prepared  as  well. 
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Nippon  Carbon  Company,  Japan 
^C-9608,  Corning  Glass  Works 
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2.2  Measurement  of  the  Thermal  Properties 


The  thermal  diffusivity  of  the  composite  and  matrix  samples  was  measured 
by  the  laser  flash  method  [5]  using  a  Nd-glass  laser.  The  specimens  used  for 
these  measurements  were  cut  from  the  hot-pressed  blanks,  and  were  in  the  form 
of  square  plates  9-12  mm  on  the  side  by  about  2  mm  thick.  Direct  transmission 
of  the  laser  beam  was  prevented  by  coating  the  specimen  surfaces  with  carbon. 
The  transient  temperature  response  of  the  specimen  rear  surface  was  monitored 
by  optical  means.  For  measurements  above  room  temperature  the  specimens  were 
held  in  an  appropriate  holder  contained  within  a  carbon  resistance  furnace  with 
nitrogen  atmosphere.  In  the  evaluation  of  the  thermal  diffusivity  from  the 
transient  temperature  response,  the  corrections  for  heat  loss  were  taken 
into  account  using  the  analysis  given  by  Heckman  [6].  Changes  in  specimen 
thickness  with  temperature  due  to  thermal  expansion  were  taken  into  account 
as  well. 

The  specific  heat  of  the  fibers  and  the  glass-ceramic  matrix  was 

* 

determined  by  differential  scanning  calorimetry  using  appropriate  equipment. 


2.3  Evaluation  of  the  Thermal  Conductivity  and  Diffusivity  of  the 
Fibers  From  Composite  Theory 


For  heat  flow  parallel  to  uniaxially  aligned  fibers,  the  thermal  con¬ 
ductivity  (K£)  of  a  composite  is: 


K 


c 


K  V 


m  m 


+  K  V 
P  P 


(1) 


where  K  is  the  thermal  conductivity,  V  is  the  volume  fraction  and  the  subscripts 

c,  m  and  p,  refer  to  the  composite,  matrix  and  fibers,  respectively. 

_ 

Dupont  990  Thermal  Analyzer. 
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For  heat  flow  perpendicular  to  the  fiber  direction,  the  thermal  conducti¬ 


vity  as  derived  by  Bruggeman  [7]  can  be  written: 
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+  K 
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(2) 


From  the  measured  value  of  the  thermal  diffusivity,  the  corresponding 
value  of  the  thermal  conductivity  (K)  can  be  calculated  from: 


K  =  <pc  (3) 

where  k  is  the  thermal  diffusivity,  p  is  density  and  c  is  the  specific 
heat.  The  specific  heat  of  the  composite  can  be  calculated  from  the  measured 
values  for  the  specific  heat  of  the  fibers  and  the  matrix  by  means  of  the 
rule  of  mixtures.  Substitution  of  the  values  for  the  thermal  conductivity 
of  the  matrix  and  the  composite  into  Eqs.  1  or  2  permits  calculations  of  the 
thermal  conductivity  of  the  fibers.  The  thermal  diffusivity  may  then  be 
determined  using  Eq.  3. 


3.  EXPERIMENTAL  RESULTS 


Figure  2  shows  the  experimental  data  for  the  specific  heat  as  a  function 
of  temperature  for  the  glass-ceramic  matrix  and  the  silicon  carbide  fibers.  The 
data  for  the  glass-ceramic  agree  very  well  with  those  obtained  by  the 
manufacture  of  a  glass-ceramic  of  similar  composition  [8] .  The  specific 
heat  of  the  silicon  carbide  fibers  significantly  exceeds  the  corresponding 
values  for  crystalline  silicon  carbide  [9].  This  difference  most  likely  can 
be  attributed  to  the  non-crystalline  nature  of  the  fibers,  since  large  relative 
differences  in  the  specific  heat  of  amorphous  and  crystalline  polymers  are 
also  observed  [10]. 


Table  1  lists  the  experimental  data  for  the  thermal  diffusivity  at 
room  temperature  for  a  number  of  composite  samples  of  different  thickness  with 
a  range  of  volume  fractions  and  both  fiber  orientations.  Also  included  in  Table 
1  are  the  values  for  the  room  temperature  thermal  conductivity  and  thermal 
diffusivity  of  the  fibers  calculated  from  the  data  for  the  specific  heat, 
density  and  thermal  diffusivity. 

Figure  3  shows  the  experimental  data  for  the  temperature  dependence  of 
the  thermal  diffusivity  of  the  glass-ceramic  matrix  and  two  composite  samples 
containing  49  vol.%  SiC  with  heat  flow  parallel  and  perpendicular  to  the  fibers. 

The  data  for  the  LAS  matrix  show  reasonable  agreement  with  data  obtained  earlier 
[ll]  for  a  sample  of  a  similar  glass-ceramic  (C9608)  supplied  by  the  manufacturer. 

Figure  4  shows  the  values  for  the  thermal  conductivity  and  thermal  diffusivity 
of  the  silicon  carbide  fibers  parallel  and  perpendicular  to  the  fiber  axis  cal¬ 
culated  from  the  experimental  data  given  in  Figs.  2  and  3. 

4.  DISCUSSION 

The  values  of  the  thermal  conductivity  and  thermal  diffusivity  of  the 
silicon  carbide  fibers  are  far  below  the  corresponding  values  of  crystalline 
silicon  carbide  [12,13].  The  low  values  for  the  present  fibers  are  most  likely 
the  result  of  their  amorphous  nature  which  limits  the  contribution  of  the 
phonon  conductivity. 

The  relatively  small  negative  temperature  dependence  of  the  thermal 
diffusivity  and  the  positive  temperature  dependence  of  the  thermal  conductivity 
of  the  fibers  is  an  indication  that  a  significant  fraction  of  the  total 
heat  conducted  occurs  by  radiative  heat  transfer.  Additional  independent 
support  for  the  latter  conclusion  is  also  provided  by  the  observations  shown  in 
Table  1,  that  the  thermal  diffusivity  of  the  composite  samples  with  the  SiC  fi¬ 
bers  increases  with  increasing  specimen  thickness  in  accordance  with  theoretical 
expectations  [14] .  Assuming  that  at  room  temperature  the  radiation  contribution 


is  negligible,  the  data  shown  in  Fig.  4  suggest  that  at  1000°C  the  radiation 
contribution  is  in  excess  of  half  of  the  total  heat  transferred.  This 
suggests  that  although  no  supporting  experimental  data  appear  to  be  available, 
these  SiC  fibers  possibly  are  excellent  transmitters  of  infra-red  radiation 
of  relatively  long  wavelength. 

The  thermal  diffusivity  and  conductivity  perpendicular  to  the  fiber  axis 
is  somewhat  less  than  parallel  to  the  fiber  axis.  Possibly,  this  could  be 
attributed  in  part  to  a  structural  anisotropy  within  the  fiber  resulting  from 
its  original  organic  structure.  No  such  preferred  orientation,  however, 
was  detected  by  Yajima  et  al  [4].  Possibly  also,  an  interfacial  contact 
resistance  due  to  a  lack  of  perfect  adhesion  or  the  factor  2-3  thermal  expansion 
mismatch  between  the  glass-ceramic  matrix  and  fibers  could  be  responsible 
for  this  effect.  This  would  manifest  itself  as  an  apparent  decrease 
of  the  thermal  conductivity  perpendicular  to  the  fiber  axis,  because  Eq.  3 
does  not  take  such  an  interfacial  resistance  into  account.  Such  an  inter¬ 
facial  resistance  would  have  a  much  smaller  (if  any)  effect  for  heat  flow 
parallel  to  the  fiber  axes.  For  this  reason  the  evaluation  of  fiber  con¬ 
ductivity  by  the  composite  method  preferably  should  rely  on  experimental 
data  using  uniaxially  aligned  composite  samples  with  heat  flow  parallel  to 
the  fiber  length. 

Some  general  remarks  are  in  order  on  the  general  feasibility  of  the 
"composite  method"  for  obtaining  heat  conduction  data  for  fibers.  It  is 
critical  to  note  that  this  method  relies  on  the  theory  for  the  thermal 
conductivity  of  composites.  Generally,  such  theory  is  exact  only  for  dilute 
concentrations  of  the  second  phase.  At  higher  volume  fractions,  at  which  the 
local  temperature  fields  around  the  disperions  will  interact,  composite 
theory  must  be  considered  to  be  an  approximation  only.  Any  deviations  between 


composite  values  calculated  from  such  approximations  and  actual  values  increase 
with  the  relative  difference  in  the  values  of  the  thermal  conductivity  of  the 
individual  phases.  For  this  reason,  increased  accuracy  of  the  value  for  the 
thermal  conductivity  of  fibers  obtained  by  the  "composite  method"  can  be 
achieved  by  choosing  a  matrix  material  with  a  value  of  the  thermal  conductivity 
as  close  as  possible  to  that  of  the  fibers.  This  condition  is  met  for  the 
composites  of  this  study.  For  the  same  reason,  the  choice  of  the  composite 
equation  used  to  calculate  the  thermal  conductivity  of  the  SiC  fibers  is 
not  critical. 

Furthermore,  for  the  "composite  method"  to  yield  reliable  data  for  the 
heat  transport  properties  of  the  fibers,  it  is  essential  that  the  scale  of 
the  microstructure  and  the  size  of  the  sample  of  the  composite  is  such  that, 
in  its  transient  response,  the  sample,  in  effect,  behaves  as  a  continuum. 
Criteria  for  this  condition  were  examined  by  Kerrisk  [15,16],  Lee  and 
Taylor  [l7],  and  Nomura  and  Chon  [l8].  For  the  samples  used  in  the  present 
study,  these  criteria  were  met. 

In  summary,  it  is  proposed  that  the  thermal  conductivity  of  fibers  can 
be  obtained  from  the  corresponding  data  for  composites  by  means  of  the  theory 
for  the  thermal  conductivity  of  composites.  The  data  presented  indicate? 

4 

the  feasibility  of  this  method. 
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Fig.  2.  Specific  heat  of  lithium-alumino-silicate  glass- 
ceramic  and  of  amorphous  silicon  carbide  fibers. 
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Experimental  data  for  the  thermal  diffusivity  as 
a  function  of  temperature  for  lithium-alumino¬ 
silicate  glass-ceramic  matrix  with  0,  48  and  49 
vol.  %  SiC  fibers  for  specimen  thicknesses  of 
1.32,  1.29  and  1.29  mm,  parallel  and  perpendicular 
to  the  fiber  axis,  respectively. 
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Fig.  A.  Temperature  dependence  of  the  thermal  con¬ 
ductivity  and  diffusivity  of  amorphous 
silicon  carbide  fibers  calculated  from  the 
composite  data  shown  in  Fig.  3. 
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CHAPTER  X 


THERMAL  STRESS  IN  MATERIALS  HEATED  INTERNALLY  3Y  RADIATION  ABSORPTION 


J.  R.  Thomas,  Jr.,  J.  P.  Singh,  and  D.  P.  H-  Kasseiman 


Department  of  Mechanical  Engineering 
and  Materials  Engineering 

Virginia  Polytechnic  Institute  and  Statu-  '’niversity 
Blacksburg,  Virginia  .^"il 


Analysis  and  results  are  reported  for  thermal  stress  in  materials  L eared  internal!.-  by 
radiation  absorption.  Two  modes  of  heating  are  considered:  radiation  suddenly  applied 
from  one  side  and  held  constant  thereafter,  and  radiation  incident  from  one  side  but 
having  a  periodic  variation  in  time. 

In  both  cases  it  is  found  that  the  radiation  absorption  coefficient  plays  a  large- 
role  in  determining  the  thermal  stress,  and  conclusions  are  drawn  concerning  desir¬ 
able  optical  thicknesses.  With  periodic  heating,  the  stress  can  cycle  between  ten¬ 
sile  and  compressive  values.  This  phenomenon  could  be  especially  significant  for 
materials  prone  to  eve  lie  fatigue  failure. 


1.  INTRODUCTION 

Structural  materials  for  fusion  reactors  will 
he  subjected  to  rather  unusual  conditions  of 
transient  heating,  depending  on  their  location 
and  exposure.  In  manv  cases  the  heating  will 
result  from  an  incident  surface  heat  flux  or 
from  the  absorption  of  nuclear  radiation,  which 
generates  .ear  internally.  The  heat  source 
will  almost  certain!  he  of  a  pulsed  nature,  so 
thermal  .cling  can  be  expected.  Manv  of  the 
mater*. lis  “'eing  considered  tor  these  applica¬ 
tions**-  include  brittle  materials,  such  as 
: ar  temperature  cerami  s,  which  are  suseep- 

-  e  ata.'.t  ropti  ic  failure  when  subjected  to. 
•  rmMent  thermal  tresses.  Previous  studies 
t  ermal  'tn  ssd!>  m  -».»cn  situations  have 
.  .m;».ed  .  'ns  i.jerat  :**n  to  seating  m  the  form  ot 
in  .n.  .  tent  -u.ir  :ac«.  c.eat  t  lux.  The  present 
lut  i.’rs  ano  coworxers  nave  recent l'*  been  studv- 
ii  rrermai  srresses  in  materials  generated  bv 
ntemal  absorption  of  radiation  .  These 
nvcsr  i  gat  ions  have  vwlded  results  or  impor¬ 
tance  to  several  app 1 1 cat  tons ,  such  as  par- 
transparent  materials  used  in  solar 
receiver.*,  and  aerospace  structures  or  compon- 
:*.t>  >  ir.  as  r .idomes  ind  !R-wi;idows. 

in  the  present  Paper,  we  wish  to  highlight  some 
•'!  the  results  of  our  previous  work  which 
appear  to  be  of  significance  in  fusion  mater¬ 
ial.-.  considerations,  and  to  r-.porr  some  more 
recent  results  for  pulsed  seating  wnich  should 
nave  more  direct  application  to  fusion  reactor 
materials.  The  nroau  objective  ot  tne  work  has 
.»ecn  to  understand  tne  role  or  the  various 
pn\sical  properties  winch  atfYct  the  nature  and 
the  magnitude  of  thermal  stress,  in  order  to 
encourage  intelligent  design  and  materials 
selection.  For  this  reason  we  nave  considered 
onlv  nlate  geometr-’  so  that  solutions  tor  the 
transient  temperature  and  thermal  stress  dis¬ 
tributions  mav  he  obtained  explicitly  bv 


analytical  means.  These  analytical  sol¬ 
utions  often  reveal  the  importance  of 
the  various  material  properties  with  greater 
clarity  chan  would  numerical  solutions  for 
complicated  geometries. 

SUDDENLY  APPLIED  RADIATION 


In  order  to  provide  background  for  understand¬ 
ing  our  result:*  for  a  periodic  heat  source,  we 
would  like  first  to  review  some  previous  re¬ 
sults.  We  have  considered  thermal  stresses 
arising  from  heating  bv  normally  incident 
radiation  arriving  on  either  side  of  a  flat 
plate  of  thickness  2a,  which  is  cooled  on  one 
side  by  convection.  The  radiation  begins 
suddenly  ac  time  t  =  h  and  is  he la  constant 
thereafter.  The  radiation  is  attenuated  expon¬ 
entially  as  it  passes  through,  the  plate,  and 
leads  to  a  volumetric  heat  gource  which  is 
mathematical ly  described  as 


u<x) 


•  lx  -  .O 


for  radiation  incident  from  the  leit.  Here 
uix)  is  the  energy  generation  rate  »  r  in :  t 
volume  at  location  x,  measured  t  rom  f.«-  cuter- 
line,  and  .  is  t:ic  linear  attenuation  c.'efri- 
cient.  For  an  absorbing  material  with  surface 

reflectance  -  , 

r 


wire  re  ;  ^  is  the  incident  radiation  flux.  For 
convenience  we  set  ■  1  -  . 

r 

To  determine  tt.e  resulting  thermal  \e 

solved  the  heat  conduction  equation 


•“T 

— 7  ( x. t  ) 


vX.t 


where  T(x,t)  is  che  temperature  at  location  x 
at  time  t,  k  the  thermal  conductivity  and  <  the 
thermal  d if f us ivitv .  The  general  solution  of 
Eq .  was  required  to  satisfy  boundary  condi¬ 

tions  appropriate  to  Che  cooling  mechanism. 

For  convective  cooling  on  Che  right,  the  condi¬ 
tions  are 

<-a,t)  *  ')  (3) 

■x 

and 

-k  ^  (a,t)  -  h  (T(a,c)  -  T  J  ,  (4) 

where  h  is  Che  convective  film  coefficient  and 
T  is  the  fluid  temperature.  Also  required  is 
an  initial  condition,  which  we  chose  for  sim¬ 
plicity  as 

T(>f,Ch  -  Tx.  (5) 

Once  T(x,t)  was  determined,  an  expression  for 
the  thermal  stresj^dist ribucion  was  developed 
from  the  equation 

\.iZ  -  fz  i-T(x.t)  +  ^  j“  T(-,c)d? 
ix  -a 

+  ,  T(C,t)id5] »  (6) 

2aJ  -a 

where  i  is  the  coefficient  of  thermal  expan¬ 
sion,  E  is  Young's  modulus,  and  v  is  Poisson's 
ratio.  The  decailed  expression  is  given  in 
He:.  7  and  will  not  he  repeated  here. 
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Figure  1.  Tine  dependence  of  Maximum  tensile 
stress  in  j  partially  absorbing  flat  plate  with 
..a  »  3,  s  ■  u.  3  W/cm^C ,  and  a  *  l  cm. 


[n  Figure  I  we  show  the  maximum  tensile  thermal 
stress,  tor  t!ie  case  described  jbove,  is  a 


function  of  time  for  various  values  of  h.  Note 
that  as  steadv  state  is  approached,  the  maximum 
stress  becomes  independent  of  h,  whereas  for 
intermediate  times,  trie  stress  ih  roughly 
proportional  to  h.  The  stress  distribution  for 
t  -  is  shown  in  Figure  J  for  various  values 
of  optical  half-thickness  -a.  Here  it  is  seen 
that  the  stress  is  tensile  at  either  face,  and 
compressive  near  the  center.  The  largest 
stress  occurs  at  the  hotter  face,  a  -.omewhat 
unexpected  result  which  :s  discussed  more 
thoroughly  in  Ref.  7.  It  is  also  apparent  that 
the  stress  is  largest  for  a  J.O. 


on. 


-i0  -o»  o«~  ~  02  o  02  oi  **  i  "c 

NON -DIMENSIONAL  DISTANCE,  t/a 

Figure  J.  Spatial  distribution  of  steady  state 
it  »  **)  thermal  stresses  for  various  values  of 
ua  in  a  partially  absorbing  flat  plate  with  k  * 
0.3  W/cm°C  and  a  *  1  cm. 

For  verv  small  values  of  h,  the  plate  :s  eft  ac¬ 
tively  insulated.  Since  no  heat  is  removed  the 
temperature  increases  mdef  initeiv,  -ut  oon 
achieves  a  t ime- invar i ant  profile.  In  Figure  \ 
we  show  the  stress  distribution  for  tn;>  case 
for  various  values  of  ...i.  Note  that  tne  >tros- 
profile  is  reversed,  as  compared  to  that  M-.ovn 
in  Figure  1  for  finite  Note  also  *:iat  the 
largest  stresses  now  occur  for  a  li»  instead 
of  ..j  *  J.i). 

The  implications  of  these  result**  tor  tt;« 
selection  of  materials  mav  oe  summarised  as 
follows:  for  non-zero  film  coefficient  h , 

intermediate  values  of  optical  thickness,  i.e., 
dwa  -  4.0,  should  be  lvoided;  ..a  "  or  a  *  * 

yield  groatlv  reduced  thermal  -tr. "-s.  fee 
fact  tiiat  the  maximum  tensile  -tress  occur->  u 
the  iiotter  surface  could  V  of  considerable 
importance  for  environmental  conditions  u-iur, 
promote  fatigue  bv  stress  o'rr  *sicn  or  -.imiiar 
median  isms . 

ihe.se  results  were  obtained  tor  a  -udd*n  M 
applied  radiation  f.eld  wnren  remain-  onstant 
after  time  t  a  We  now  treat  in  mor<  ietail 
the  case  of  an  oscillating  -ourc.-  of  rao  s  it  ion 
which  should  ^e  mor<-  r.  - r**-»  ntat  i  :  •  .  :o:i 
reactor  apnl  iv:at  ion-. 
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NON  -OIMENSONAL  INSTANCE  .  «/o 

Figure  3.  Spatial  distribution  of  steady  state 
(t  «  *)  thermal  stresses  for  various  values  of 
ua,  h  -  0,  k  •  0.3  W/cm°C,  and  a  *  1  era. 


3.  OSCILLATING  RADIATION  SOURCE 
3.1  Analysis 

We  consider  a  slab  as  before,  but  subjected  to 
radiation  from  the  left  which  oscillates  in 
time  with  frequency  w.  We  thus  seek  the  solu¬ 
tion  of  Eq.  (2)  subject  to  the  same  boundary 
and  initial  conditions  (3-5),  but  with  a  heat 
source 

u (x)  *  uQe'u(x  +  a,(l  +  1  sin(.t) ] .  (7) 

Here  '  ^  1  is  the  amplitude  of  the  oscillation 
in  the  heat  generation  rate,  and  all  other 
svmbols  are  as  previously  defin^j.  We  find  the 
solution  by  standard  techniques  to  be 

T(x,t)  -  T  +  :  B  (C)  cos  (•  x),  (8) 

*  ,  n  n 

n«l 


where 


(9a) 


(9b) 


N  *  — —  {sin (2'  a) cos (2-*  a)  +  1'  a], 
n  2  •  n  n  n 

n 

and  where  the  eigenvalues  are  the  solution  of 
the  transcendental  equation 

,  h 

•.  tan(2*  a)  *  -  . 
n  n  k 

In  Eq.  (9b),  .■  and  c  are  the  density  and  heat 
capacity  of  the  material,  respectively.  The 
temperature  distribution  (8),  when  substituted 
into  Eq.  (6),  yields  an  explicit  expression  for 
the  thermal  stress  distribution.  The  result 
is  identical  to  that  given  in  Ref.  7  with 
B  (t)  given  by  Eqs.  (9)  above,  and  thus  will 
not  be  repeated  here. 

3.2  Discussion 

The  temperature  distribution  given  by  Eqs.  (8) 
and  (9)  is  very  similar  to  that  found  for 
suddenly  applied  radiation,  except  that  a 
sinusoidal  variation  is  superimposed  on  the 
time  variation.  The  amplitude  of  tills  oscil¬ 
lation  is  greatest  in  the  front  face  of  the 
plate  (x  *  -  a)  of  course,  and  there  is  a  phase 
shift  between  the  temperature  oscillations  at 
the  front  and  back  face.  This  difference  in 
amplitude  can  be  seen  in  Figure  4,  where  the 
envelope  of  the  temperature  oscillation  is 
plotted  in  nondimens ional  form  at  long  times  - 
i.e.  near  steady-state  conditions.  This  phase 
shift  is  reflected  in  the  thermal  stress  dis¬ 
tribution  which  is  shown  as  a  function  of  time 
in  nondimens ional  form  in  Figure  5.  This  plot 
was  also  constructed  to  show  the  long-time 
behavior  when  near  steady-state  conditions 
prevail.  Although  the  amplitude  of  the  stress 
oscillation  is  nearly  the  same  at  the  front 
face  as  in  the  center,  at  the  center  the 
oscillation  is  between  tensile  and  compressive 
values,  whereas  the  stress  in  the  front  face 
remains  tensile  throughout  the  oscillation.  Of 
potentially  greatest  importance  is  the  large 
swing  in  the  stress  on  the  r*ar  face,  which 
also  oscillates  between  tension  and  compression. 
Materials  prone  to  cyclic  fatigue  failure  are 
thus  more  likely  to  fail  at  the  back  face  of 
the  plate  which  is  in  contact  with  the  cooling 
medium. 

It  is  recognized  that  the  operating  c*-’cle  of  a 
fusion  reactor  will  likely  consist  of  flat- 
topped  pulses  separated  by  somewhat  longer 
downtimes,  and  is  thus  not  accuratel  ■  modeled 
by  a  sinusoidal  variation.  However  the  purpose 
of  the  present  study  was  to  investigate  Che 
effects  of  periodically  applied  radiation  in 
the  simplest  possible  mathematical  format. 

Future  studies  will  consider  more  realistic 
radiation-source  operating  cycles. 
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NON-DIMENSIONAL  OlSTANCE.  i/O 

Figure  4.  Envelope  of  the  non-dimensional 
temperature  oscillation;  ua  =  3.0,  h  *  0.006 
W/ca^C. 
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NON  DIMENSIONAL  TIME.  *t/o* 


Figure-  3.  Time  depence  of  the  thermal  stress; 
•  a  »  ).<>,  h  «  0.00b  W/cm**C. 


*.  SUMMARY  AND  CONCLUSIONS 


Che  right,  the.-  maximum  steady-state  stress 
occurs  on  the  side  j'aciag  the  radiation  source. 
In  addition  to  the  usual  dependence  on  E,  ,  t, 
thickness  a,  and  thermal  difjusivity  •  ,  the 
magnitude  of  the  stresses  is  rather  stronglv 
dependent  on  the  absorption  coefficient,  with 
intermediate  values  of  the  optical  thickness 
producing  greater  stresses.  Stresses  resulting 
from  an  oscillating  radiation  source  also  have 
an  oscillator^  character,  with  the  maximum 
swing  at  the  rear  face  vnere  the  plate  la 
cooled.  Both  near  the  center  and  at  the  rear 
face  the  stresses  oscillate  between  tension  and 
compression,  which  could  have  important  impli¬ 
cations  for  materials  in  which  fatigue  failure 
is  a  problem. 
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